Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 954,436 T→C 6.7% L92P (CTG→CCG)  cyoE → protoheme IX farnesyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947954,4360TC 6.7% 665.8 / 24.1 208L92P (CTG→CCG) cyoEprotoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base T (106/88);  new base C (7/7);  total (113/95)
Fisher's exact test for biased strand distribution p-value = 7.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 213 total aligned reads displayed.

TGATCGGTACTTCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCGGCAGCCTGT  >  NC_002947/954289‑954585
                                                                                                                                                   |                                                                                                                                                     
tgATCGGTACTTCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGc                                                                                                                                                     >  2:4356127/1‑150 (MQ=255)
  aTCGGTACTTCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTc                                                                                                                                                   >  1:3078858/1‑150 (MQ=255)
     ggTACTTCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCggg                                                                                                                                                <  1:2969356/150‑1 (MQ=255)
     ggTACTTCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCggg                                                                                                                                                <  2:2697804/150‑1 (MQ=255)
        aCTTCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTg                                                                                                                                             >  2:2192606/1‑150 (MQ=255)
          ttCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGc                                                                                                                                           >  1:3560386/1‑150 (MQ=255)
          ttCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGc                                                                                                                                           >  2:266391/1‑150 (MQ=255)
          ttCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCCGAGCAGGGTGGc                                                                                                                                           <  1:3774348/150‑1 (MQ=255)
           tCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCg                                                                                                                                          <  2:3323367/150‑1 (MQ=255)
           tCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCCGAGCAGGGTGGCg                                                                                                                                          <  2:852935/150‑1 (MQ=255)
             gCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCCGAGCAGGGTGGCgta                                                                                                                                        <  2:1879878/150‑3 (MQ=255)
              cTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGt                                                                                                                                       <  1:3568432/150‑1 (MQ=255)
              cTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATAGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGt                                                                                                                                       <  2:560985/150‑1 (MQ=255)
               tggtggTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGtt                                                                                                                                      >  1:2050840/1‑150 (MQ=255)
               tggtggTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGtt                                                                                                                                      >  1:1155007/1‑150 (MQ=255)
                 gtggtTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTc                                                                                                                                    >  2:2476255/1‑150 (MQ=255)
                  tggtTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCa                                                                                                                                   <  2:4101406/150‑1 (MQ=255)
                  tggtTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCa                                                                                                                                   <  1:106340/150‑1 (MQ=255)
                  tggtTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCa                                                                                                                                   <  1:2409540/150‑1 (MQ=255)
                  tggtTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGAAACCCGCCACCCCGAGCAGGGTGGCgtagatca                                                                                                                                   <  1:3456548/150‑8 (MQ=38)
                  tggtTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCGGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCa                                                                                                                                   <  1:101364/150‑1 (MQ=255)
                       gCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCctg                                                                                                                              >  1:4469728/1‑150 (MQ=255)
                         gTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCctgct                                                                                                                            <  2:1518398/150‑1 (MQ=255)
                                    tgtTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTGATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAgg                                                                                                                 >  2:21197/1‑150 (MQ=255)
                                    tgtTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAgg                                                                                                                 >  2:3947782/1‑150 (MQ=255)
                                    tgtTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAgg                                                                                                                 <  2:3533693/150‑1 (MQ=255)
                                      ttCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGcc                                                                                                               >  2:1553692/1‑150 (MQ=255)
                                         aacaacTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAAc                                                                                                            >  2:1730442/1‑150 (MQ=255)
                                         aacaacTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAAc                                                                                                            >  2:2929608/1‑150 (MQ=255)
                                         aacaacTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGACCAAc                                                                                                            >  2:140805/1‑150 (MQ=255)
                                          acaacTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAAcc                                                                                                           >  2:3814520/1‑150 (MQ=255)
                                            aacTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCg                                                                                                         <  1:870626/150‑1 (MQ=255)
                                               tGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTg                                                                                                      >  1:2666844/1‑150 (MQ=255)
                                                 cATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGATTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTc                                                                                                    >  2:1851964/1‑150 (MQ=255)
                                                   tCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTg                                                                                                  >  2:2653126/1‑150 (MQ=255)
                                                      aCCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGt                                                                                               >  1:1535629/1‑150 (MQ=255)
                                                      aCCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACACCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGt                                                                                               >  1:2137750/1‑150 (MQ=255)
                                                        cGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTc                                                                                             >  2:1458347/1‑150 (MQ=255)
                                                         gTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCt                                                                                            <  2:1032371/150‑1 (MQ=255)
                                                          tGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTg                                                                                           <  1:4152779/150‑1 (MQ=255)
                                                          tGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTg                                                                                           <  2:522648/150‑1 (MQ=255)
                                                           gACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTgc                                                                                          <  1:711205/150‑1 (MQ=255)
                                                              aTCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCa                                                                                       <  1:3672685/150‑1 (MQ=255)
                                                               tCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCAc                                                                                      >  2:3045064/1‑150 (MQ=255)
                                                               tCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCAc                                                                                      >  2:133659/1‑150 (MQ=255)
                                                                   ccACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGAt                                                                                  >  2:3985454/1‑150 (MQ=255)
                                                                   ccACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGAt                                                                                  >  2:2286259/1‑150 (MQ=255)
                                                                    cacaAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGAtt                                                                                 >  2:111529/1‑150 (MQ=255)
                                                                    cacaAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGAtt                                                                                 <  1:1202816/150‑1 (MQ=255)
                                                                     acaAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTg                                                                                <  2:59279/150‑1 (MQ=255)
                                                                     acaAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTg                                                                                >  1:2364572/1‑150 (MQ=255)
                                                                      caAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTgg                                                                               <  2:4292095/150‑1 (MQ=255)
                                                                       aaGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGc                                                                              <  1:677247/150‑1 (MQ=255)
                                                                       aaGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGc                                                                              <  1:2892169/150‑1 (MQ=255)
                                                                       aaGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGc                                                                              <  2:1242061/150‑1 (MQ=255)
                                                                          aTGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCTCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCttt                                                                           <  2:3434997/150‑1 (MQ=255)
                                                                          aTGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCttt                                                                           <  1:2497959/150‑1 (MQ=255)
                                                                           tGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTg                                                                          <  2:3287449/150‑1 (MQ=255)
                                                                              aGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTGTgtcg                                                                       >  2:2481255/1‑150 (MQ=255)
                                                                                cgcACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTgtcgtc                                                                     <  2:3373750/150‑1 (MQ=255)
                                                                                cgcACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTgtcgtc                                                                     <  1:3638436/150‑1 (MQ=255)
                                                                                   aCCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTAc                                                                  <  1:1263873/150‑1 (MQ=255)
                                                                                   aCCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCGTTGTCGTCTAc                                                                  >  2:1387123/1‑150 (MQ=255)
                                                                                        gAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCgg                                                             <  2:2569809/150‑1 (MQ=255)
                                                                                           ccGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGttt                                                          <  2:1075944/150‑1 (MQ=255)
                                                                                           ccGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGttt                                                          >  1:323934/1‑150 (MQ=255)
                                                                                           ccGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGttt                                                          <  2:4335483/150‑1 (MQ=255)
                                                                                           ccGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACGGATTGGCTTTGTCGTCTACGTCGGttt                                                          <  1:4356808/150‑1 (MQ=255)
                                                                                            cGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTc                                                         >  2:3767652/1‑150 (MQ=255)
                                                                                                tCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACa                                                     >  2:649146/1‑150 (MQ=255)
                                                                                                  aTGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGc                                                   >  2:3697445/1‑150 (MQ=255)
                                                                                                   tGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACa                                                     >  1:2173082/1‑147 (MQ=255)
                                                                                                       tCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCtgtg                                              <  1:3330519/150‑1 (MQ=255)
                                                                                                          gggCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCt                                           >  1:4264224/1‑150 (MQ=255)
                                                                                                            gcggcATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGa                                         <  1:4303157/150‑1 (MQ=255)
                                                                                                             cggcATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGaa                                        <  1:2467768/150‑1 (MQ=255)
                                                                                                             cggcATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGaa                                        <  2:45573/150‑1 (MQ=255)
                                                                                                             cggcATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGaa                                        <  2:4445448/150‑1 (MQ=255)
                                                                                                                      gCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGAACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATc                               >  1:165957/1‑150 (MQ=255)
                                                                                                                        tGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCgg                             >  1:3300890/1‑150 (MQ=255)
                                                                                                                         gCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGt                            <  2:1639245/150‑1 (MQ=255)
                                                                                                                            ccTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCa                         <  1:1230525/150‑1 (MQ=255)
                                                                                                                              tCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACg                       <  2:217071/150‑1 (MQ=255)
                                                                                                                              tCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACg                       >  2:1935806/1‑150 (MQ=255)
                                                                                                                                 cgcTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCa                    >  2:2437905/1‑150 (MQ=255)
                                                                                                                                    tgaaacccGCCACCCCGAGCAGGGTGGCGTGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCt                 >  2:2900094/8‑150 (MQ=255)
                                                                                                                                    tGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCt                 >  2:728357/1‑150 (MQ=255)
                                                                                                                                    tGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCt                 >  1:846222/1‑150 (MQ=255)
                                                                                                                                       tCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTgg              >  2:3134917/1‑150 (MQ=255)
                                                                                                                                         tACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTc            >  1:63871/1‑150 (MQ=255)
                                                                                                                                          ccGCCACCCCGAGCAGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCg           >  1:208924/2‑150 (MQ=255)
                                                                                                                                          aCGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCg           >  2:3619394/1‑150 (MQ=255)
                                                                                                                                           cGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCgg          >  1:3489509/1‑150 (MQ=255)
                                                                                                                                           cGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCgg          >  1:3995084/1‑150 (MQ=255)
                                                                                                                                           cGCCACCCCGAGCAGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCgg          >  2:1063297/1‑150 (MQ=255)
                                                                                                                                              cACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCGGCAg       >  2:3198824/1‑150 (MQ=255)
                                                                                                                                               aCCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCGGCAGc      <  2:311092/150‑1 (MQ=255)
                                                                                                                                                 ccTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCGGCAGCCt    >  1:4447405/1‑150 (MQ=255)
                                                                                                                                                  cTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCGGCAGCCTg   <  2:4469728/150‑1 (MQ=255)
                                                                                                                                                   tGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCGGCAcccgac  <  1:1675337/150‑4 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGATCGGTACTTCGCTGGTGGTTGCGTCCGGATGCGTGTTCAACAACTGCATCGACCGTGACATCGACCACAAGATGGAGCGCACCAAGAACCGTGTCATGGTTCAGGGCGGCATGTCGCTGCCCCTCGCGCTGATCTACGCCACCCTGCTCGGGGTGGCGGGTTTCAGCCTGCTGTATGTCCAGGCCAACCCGCTGTCTGCGTTCTGCGCACTGATTGGCTTTGTCGTCTACGTCGGTTTCTACAGCCTGTGGCTGAAGCGTAAATCGGTGCACGGCACCTTGGTCGGCAGCCTGT  >  NC_002947/954289‑954585

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