Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,033,421 C→A 6.4% L287M (CTG→ATG)  pgi‑I → glucose 6‑phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,033,4210CA 6.4% 211.9 / 4.5 78L287M (CTG→ATG) pgi‑Iglucose 6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base C (35/38);  new base A (2/3);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

GTACCGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGACGCAAACAGC  >  NC_002947/2033282‑2033570
                                                                                                                                           |                                                                                                                                                      
gTACCGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACGGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCgg                                                                                                                                              <  1:92888/150‑1 (MQ=21)
 tACCGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGc                                                                                                                                             <  2:1622584/150‑1 (MQ=34)
 cccggCCACTTCATTGCCTTACCCAGCAACAAGGCCGTGGTTGTAGCTTGCGGTATCCGGGAAGAAAACATCTCCCCGATGTGGGCCCGGGGGGGGGGCCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGc                                                                                                                                             <  1:1495758/146‑1 (MQ=11)
    cGCCACTTCATTGCCAAAGCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACt                                                                                                                                          >  1:559115/1‑150 (MQ=255)
      ccACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGc                                                                                                                                        >  1:1017352/1‑150 (MQ=35)
         cTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCTTGCCGAACGCCCTGGCCATCGGCACtagcaa                                                                                                                                     >  2:504024/1‑145 (MQ=255)
         cTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATGGCCAGGGCGATCGGc                                                                                                                                             >  2:1449122/1‑142 (MQ=11)
         cTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATGGCCAGGGCGATCGGc                                                                                                                                             <  1:1449122/142‑1 (MQ=11)
         cTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCaa                                                                                                                                     >  1:1115049/1‑150 (MQ=34)
             aTTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTg                                                                                                                                                      >  2:653256/1‑129 (MQ=35)
             aTTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTg                                                                                                                                                      <  1:653256/129‑1 (MQ=35)
              ttGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGGATCCGCGAAGAAAACATCTTCCCGATGGGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCa                                                                                                                                >  1:594696/1‑150 (MQ=25)
                gCCGTATCCAGAAACAAGGCCGCTGCTGGAGCGTTCGGTATCCGGGAAGAAAACATCTTCCCGATGTGGCGTAGTAACCAAGGCCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAg                                                                                                                              <  2:559115/150‑1 (MQ=11)
                cccGTATCCAGCAACAAGGCCGCTGCTGAAGCGTTCGGTAGCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAg                                                                                                                              <  2:751978/149‑1 (MQ=21)
                         aGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAActgctg                                                                                                                     <  2:1152603/150‑1 (MQ=34)
                         aGCAACAAGGCCGCTGCTGTAGCGTTCGGTACCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAActgctg                                                                                                                     >  1:1513386/1‑150 (MQ=25)
                         aGCAACAAGGCAGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAActgctg                                                                                                                     >  1:899348/1‑150 (MQ=25)
                                 ggCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGc                                                                                                            <  2:1017352/150‑1 (MQ=37)
                                  gCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAAGGGTTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGcc                                                                                                           >  1:646924/1‑150 (MQ=255)
                                       tgctgTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCCTGATGACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTacac                                                                                                      >  1:991437/1‑150 (MQ=11)
                                       tgctgTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTacac                                                                                                      >  1:157539/1‑150 (MQ=39)
                                          tgTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCAt                                                                                                   >  1:1480646/1‑150 (MQ=38)
                                          tgTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCATGGCCATCGGCACTg                                                                                                                                         >  2:736738/1‑113 (MQ=255)
                                          tgTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCATGGCCATCGGCACTg                                                                                                                                         <  1:736738/113‑1 (MQ=25)
                                           gTAGCGTTCGGTAGCCGCGAAAAAAACATCTCCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATg                                                                                                  <  2:1097169/150‑1 (MQ=255)
                                              gCGTGCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGGGGGTGGGCGCTACTCGCTGTGGTCGACCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGAc                                                                                               <  1:218914/150‑1 (MQ=255)
                                                       aTCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTc                                                                                      >  2:1193869/1‑150 (MQ=38)
                                                           gcgAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGa                                                                                  >  1:1462120/1‑150 (MQ=39)
                                                            cgAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCCGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGAc                                                                                 <  1:472619/150‑1 (MQ=38)
                                                                aaaaaCATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGGCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCa                                                                             <  2:655737/149‑1 (MQ=38)
                                                                 aaaaCATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCa                                                                                             >  2:988078/1‑133 (MQ=40)
                                                                 aaaaCATCTTCCCGATGTGGGACGGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCTCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCa                                                                                             <  1:988078/133‑1 (MQ=255)
                                                                  aaaCATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCAcc                                                                           >  2:613172/1‑150 (MQ=38)
                                                                   aaCATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCg                                                                          >  2:1478395/1‑150 (MQ=38)
                                                                     cATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGtt                                                                        >  2:1246571/1‑150 (MQ=38)
                                                                     cATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGtt                                                                        >  2:935721/1‑150 (MQ=39)
                                                                             cGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGa                                                                >  2:912064/1‑150 (MQ=38)
                                                                             cGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCATGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGa                                                                >  1:827790/1‑150 (MQ=38)
                                                                              gATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGaa                                                               >  1:1512128/1‑150 (MQ=37)
                                                                                    ggACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGACATCGGAATCCCGATCGCCCTATCCGTCGGCACAGCAAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGcc                                                         >  1:1693557/1‑150 (MQ=11)
                                                                                     gACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCg                                                        <  1:935721/150‑1 (MQ=35)
                                                                                     gACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCg                                                        >  2:62571/1‑150 (MQ=35)
                                                                                     gACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCg                                                        <  2:743899/150‑1 (MQ=35)
                                                                                      aCTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCgg                                                       >  2:661999/1‑150 (MQ=35)
                                                                                            tgggtgggCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCTCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTAc                                                 <  2:1047326/150‑1 (MQ=25)
                                                                                            tgggtgggCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGaa                                                               >  2:1354173/1‑136 (MQ=35)
                                                                                            tgggtgggCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGAGCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGaa                                                               <  1:1354173/136‑1 (MQ=25)
                                                                                                ggggCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGc                                             <  1:1100956/149‑1 (MQ=17)
                                                                                                 gggCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGcc                                            <  1:661999/150‑1 (MQ=32)
                                                                                                      cTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCctgct                                       >  2:1365140/1‑150 (MQ=32)
                                                                                                       tACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCctgctg                                      <  1:1478395/150‑1 (MQ=32)
                                                                                                       tACTCGCTGTGGTCGGCCATCGGCCTGCCGAGCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCctgctg                                      <  1:912064/150‑1 (MQ=17)
                                                                                                        aCTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGGACTTCCAGACCCCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTgg                                     <  2:594696/150‑1 (MQ=17)
                                                                                                          tCGCTGTGGTCGGCCATCGGCTTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGc                                   <  1:504024/150‑1 (MQ=17)
                                                                                                          tCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCg                                                        <  1:281960/129‑1 (MQ=33)
                                                                                                          tCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCg                                                        >  2:281960/1‑129 (MQ=33)
                                                                                                           cGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCg                                                                               <  1:1465618/105‑1 (MQ=32)
                                                                                                           cGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCg                                                                               >  2:1465618/1‑105 (MQ=32)
                                                                                                            gCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCgt                                 >  2:868230/1‑150 (MQ=32)
                                                                                                                    cGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATg                         <  1:868230/150‑1 (MQ=32)
                                                                                                                       ccATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCa                      >  2:1602024/1‑150 (MQ=32)
                                                                                                                          tCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACt                   >  1:1494420/1‑150 (MQ=32)
                                                                                                                              gctctcGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAGCTGCTGTCCGGTGCCTACACCCGGTACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTg               <  1:28146/145‑1 (MQ=11)
                                                                                                                              cgtGCCGATCGCCCTGGCCATCGGCACTGCCCACTTACATGGACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTg               <  2:646924/148‑1 (MQ=11)
                                                                                                                              ccTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTg               <  1:1902494/150‑1 (MQ=32)
                                                                                                                               cTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTgg              <  1:62571/150‑1 (MQ=32)
                                                                                                                               cTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTgg              <  2:583516/150‑1 (MQ=32)
                                                                                                                               cTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTgg              <  1:1193869/150‑1 (MQ=32)
                                                                                                                               cTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTgg              <  2:1512128/150‑1 (MQ=32)
                                                                                                                               cTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTgg              <  2:1159096/150‑1 (MQ=32)
                                                                                                                                 gCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGa            <  2:1462120/150‑1 (MQ=32)
                                                                                                                                 gCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGa            <  2:827790/150‑1 (MQ=32)
                                                                                                                                  ccGATGGCCAGGGCGATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCg                                                        <  2:19119/105‑1 (MQ=11)
                                                                                                                                  ccGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGAc                                                                                               <  1:429818/66‑1 (MQ=37)
                                                                                                                                  ccGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGAc                                                                                               >  2:429818/1‑66 (MQ=37)
                                                                                                                                    gATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCACCTTCt                >  1:75332/1‑143 (MQ=17)
                                                                                                                                    gATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCt                <  2:75332/143‑1 (MQ=32)
                                                                                                                                    gATCGACCTGGATCTCGGCACTGAAAATTTCAAGGACCCGCTGTGCCGGTTCGACCCCATGAACCAGAACTTTCAGACCGCACCGTTCGACAAGCACATGCCGGTGTTACTGGCCCTGCTGGGCGTGTGGTATGGCGACTTCTggagact         >  2:28146/1‑145 (MQ=2)
                                                                                                                                        ggccTGGCCATCGGCACTGCCAACTTCAAGCAATTGCTGTACGGTGCCTAGACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTTCTGGTCGTGTAGTATGGCAACTTCTGGGACGCAAAc     <  2:907481/148‑1 (MQ=11)
                                                                                                                                         cccTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGACGCAAACa    >  1:515243/1‑150 (MQ=32)
                                                                                                                                          ccTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGACGCAAACAg   >  1:796923/1‑150 (MQ=32)
                                                                                                                                           cTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGACGCAAACAGc  >  1:951764/1‑150 (MQ=32)
                                                                                                                                           |                                                                                                                                                      
GTACCGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCGATGTGGGACTGGGTGGGTGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGTGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGACGCAAACAGC  >  NC_002947/2033282‑2033570

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: