Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,671,594 G→A 5.2% A124T (GCA→ACA)  hcaB → 2,3‑dihydroxy‑2,3‑dihydrophenylpropionate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,671,5940GA 5.2% 253.5 / 2.1 78A124T (GCA→ACA) hcaB2,3‑dihydroxy‑2,3‑dihydrophenylpropionate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (34/39);  new base A (2/2);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

CTGACTCGTTCCGGCAAGCTGGATTGTTTTATCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATG  >  NC_000913/2671447‑2671741
                                                                                                                                                   |                                                                                                                                                   
cTGACTCGTTCCGGCAAGCTGGATTGTTTTATCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCTGTTACCTGCTGGgcg                                                                                                                                                     >  3:242857/1‑148 (MQ=255)
cTGACTCGTTCCGGCAAGCTGGATTGTTTTATCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGgcgc                                                                                                                                                    <  4:279482/149‑1 (MQ=255)
   aCTCGTTCCGGCAAGCTGGATTGTTTTATCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCaa                                                                                                                                                  >  8:191547/1‑148 (MQ=255)
   aCTCGTTCCGGCAAGCTGGATTGTTTTATCGGCAATGCAGGCATCTGGGATCACAATGACTCACTGGTTAATACTCCCGCAGAGACGCTCGAACCCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCaa                                                                                                                                                  >  8:143646/1‑148 (MQ=255)
    cTCGTTCCGGCAAGCTGGATTGTTTTATCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCaaaa                                                                                                                                                <  3:29879/149‑1 (MQ=255)
       gTTCCGGCAAGCTGGATTGTTTTATCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGc                                                                                                                                              <  8:187135/148‑1 (MQ=255)
       gTTCCGGCAAGCTGGATTGTTTTATCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGc                                                                                                                                              <  2:219852/148‑1 (MQ=255)
        ttCCGGCAAGCTGGATTGTTTTATCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCt                                                                                                                                            >  5:242060/1‑149 (MQ=255)
          ccGGCAAGCTGGATTGTTTTATCGGCAATGCAGGCATCTGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTgc                                                                                                                                          <  3:99338/148‑1 (MQ=255)
           cGGCAAGCTGGATTGTTTTATCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTgc                                                                                                                                          <  8:163017/148‑1 (MQ=255)
              cAAGCTGGATTGGTTTATCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTc                                                                                                                                      >  8:307090/1‑149 (MQ=255)
                aGCTGGATTGTTTTATCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCACAAAAGCCTGCGCTCCg                                                                                                                                    >  7:221965/1‑149 (MQ=255)
                     gATTGTTTTATCGGCAATGTAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGt                                                                                                                                >  8:933/1‑148 (MQ=255)
                       ttGTTTTATCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTaa                                                                                                                             <  3:240548/149‑1 (MQ=255)
                           tttATCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGc                                                                                                                         >  4:279848/1‑149 (MQ=255)
                              aTCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCa                                                                                                                       <  7:205925/148‑1 (MQ=255)
                               tCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTCCCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGt                                                                                                                     <  5:181010/149‑1 (MQ=255)
                               tCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCt                                                                                                                                            <  2:219857/126‑1 (MQ=255)
                               tCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCt                                                                                                                                            >  1:219857/1‑126 (MQ=255)
                               tCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGt                                                                                                                     <  3:337875/149‑1 (MQ=255)
                                      tgtaGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGgca                                                                                                              <  7:933/146‑1 (MQ=255)
                                        cAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGgcagc                                                                                                            <  2:173597/149‑1 (MQ=255)
                                             aTCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGAt                                                                                                       <  1:20113/149‑1 (MQ=255)
                                              tCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGAtt                                                                                                      >  1:214249/1‑149 (MQ=255)
                                              tCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGAtt                                                                                                      >  7:140645/1‑149 (MQ=255)
                                                  ggATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGAtttt                                                                                                    <  2:214249/147‑1 (MQ=255)
                                                  ggATCACAATGCCTCACTGGGTACTACTCCCGCAGAGACGCTCGAAACCTGCTTCCACGAGCGGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTca                                                                                                  <  4:242857/149‑1 (MQ=255)
                                                    aTCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGgcagca                                                                                                           >  8:200455/1‑138 (MQ=255)
                                                    aTCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGgcagca                                                                                                           <  7:200455/138‑1 (MQ=255)
                                                     tCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTGAATCGCCAGTGAAg                                                                                                                 >  6:216881/1‑131 (MQ=255)
                                                       acaATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTg                                                                                             >  8:201874/1‑149 (MQ=255)
                                                        caATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCaaaa                                                                                                                                                >  1:184930/1‑97 (MQ=255)
                                                        caATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCaaaa                                                                                                                                                <  2:184930/97‑1 (MQ=255)
                                                             ccTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAAt                                                                                       <  3:167966/149‑1 (MQ=255)
                                                             ccTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAAt                                                                                       <  5:275885/149‑1 (MQ=255)
                                                             ccTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCACAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCaaa                                                                                        <  8:221965/148‑1 (MQ=255)
                                                              cTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCACAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTc                                                                                           <  5:247349/144‑1 (MQ=255)
                                                              cTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCACAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTc                                                                                           >  6:247349/1‑144 (MQ=255)
                                                               tCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATgc                                                                                     >  6:276744/1‑149 (MQ=255)
                                                                   tGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTTATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATgccgcc                                                                                 >  6:296710/1‑149 (MQ=255)
                                                                     gagaaTACTCCCGCAGATACGCTCGACACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCGTGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCAGACTGTCAAATgccgcc                                                                                 <  8:318979/144‑1 (MQ=255)
                                                                     gTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCt                                                                                <  6:242060/148‑1 (MQ=255)
                                                                      ttAATACTCCCTCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCt                                                                                <  7:1015/147‑1 (MQ=255)
                                                                      ttAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCt                                                                                >  8:1015/1‑147 (MQ=255)
                                                                          tACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATc                                                                          <  7:307090/149‑1 (MQ=255)
                                                                           aCTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATcc                                                                         >  3:267857/1‑149 (MQ=255)
                                                                             tCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTg                                                                       >  7:153743/1‑149 (MQ=255)
                                                                                      gaCGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTg                                                              <  4:278214/149‑1 (MQ=255)
                                                                                      gaCGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCAGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGtc                                                              >  5:316/1‑148 (MQ=255)
                                                                                        cGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTACCCTGGCGGCGGTGGc                                                            <  4:144137/149‑1 (MQ=255)
                                                                                         gCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATc                                                                          >  7:134869/1‑134 (MQ=255)
                                                                                         gCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATc                                                                          <  8:134869/134‑1 (MQ=255)
                                                                                         gCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTGGcc                                                           >  8:76451/1‑149 (MQ=255)
                                                                                              aaaCCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCg                                                                                                                          <  1:279958/81‑1 (MQ=255)
                                                                                              aaaCCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCg                                                                                                                          >  2:279958/1‑81 (MQ=255)
                                                                                              aaaCCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTGGCCCGCTg                                                      >  3:74540/1‑149 (MQ=255)
                                                                                              aaaCCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCAGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTGGCCCGCTg                                                      <  7:350939/149‑1 (MQ=255)
                                                                                               aaCCGGCTTCCTCGAGCTGTTTAACGTCAATGTTCTCGGTTACATGCTGGGCGCAATAGCCTGCGCTCCGGCGt                                                                                                                                <  6:71753/74‑1 (MQ=255)
                                                                                                       tCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTgg                                                             >  2:232094/1‑133 (MQ=255)
                                                                                                          aCGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGc                                                                                                                         <  6:211514/70‑1 (MQ=255)
                                                                                                          aCGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGc                                                                                                                         >  5:211514/1‑70 (MQ=255)
                                                                                                              gCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTaaa                                       >  2:217917/1‑148 (MQ=255)
                                                                                                               cTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACa                                     >  5:9537/1‑149 (MQ=255)
                                                                                                                 ttttAACTTCAATGTTCTCTTTCACCTGCAGGGCTCAAAAGCCTGCGCTCCGTCGTTAATCCCCAGTGAAGGCAGCATTATTTTCACACTGTCAAATGCCTCCTGGTATCCTGGCGGCGGTGGCCCGATGTACACCGCCAGTAAACATg                                   <  6:290872/148‑1 (MQ=255)
                                                                                                                  tttAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATg                                   <  4:267857/148‑1 (MQ=255)
                                                                                                                        gTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAAc                            >  4:179126/1‑149 (MQ=255)
                                                                                                                            aTGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGa                        >  8:204089/1‑149 (MQ=255)
                                                                                                                             tGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTTATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGAc                       <  5:296710/149‑1 (MQ=255)
                                                                                                                             tGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGcc                                 >  6:105358/1‑139 (MQ=255)
                                                                                                                             tGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGcc                                 <  5:105358/139‑1 (MQ=255)
                                                                                                                             tGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGAc                       <  7:76451/149‑1 (MQ=255)
                                                                                                                                tctcGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTa                    >  2:156916/1‑149 (MQ=255)
                                                                                                                                    ggTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTGGCCCTCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTc                 <  8:36604/148‑1 (MQ=255)
                                                                                                                                    ggTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTc                 <  3:279848/148‑1 (MQ=255)
                                                                                                                                         ccTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAAc           >  4:263200/1‑149 (MQ=255)
                                                                                                                                            gctgGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAAc                            >  7:6931/1‑129 (MQ=255)
                                                                                                                                            gctgGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAAc                            <  8:6931/129‑1 (MQ=255)
                                                                                                                                             ctgGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGc       <  6:181749/149‑1 (MQ=255)
                                                                                                                                              tgGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCt      >  2:363377/1‑149 (MQ=255)
                                                                                                                                                ggCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTa    <  5:276744/149‑1 (MQ=255)
                                                                                                                                                   gcAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATg  >  6:85060/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CTGACTCGTTCCGGCAAGCTGGATTGTTTTATCGGCAATGCAGGCATCTGGGATCACAATGCCTCACTGGTTAATACTCCCGCAGAGACGCTCGAAACCGGCTTCCACGAGCTGTTTAACGTCAATGTTCTCGGTTACCTGCTGGGCGCAAAAGCCTGCGCTCCGGCGTTAATCGCCAGTGAAGGCAGCATGATTTTCACACTGTCAAATGCCGCCTGGTATCCTGGCGGCGGTGGCCCGCTGTACACCGCCAGTAAACATGCCGCAACCGGACTTATTCGCCAACTGGCTTATG  >  NC_000913/2671447‑2671741

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: