Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 250,912 G→C 6.2% A97A (GCC→GCG nfeD ← membrane bound peptidase NefD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947250,9120GC 6.2% 386.1 / 12.6 123A97A (GCC→GCGnfeDmembrane bound peptidase NefD
Reads supporting (aligned to +/- strand):  ref base G (51/64);  new base C (4/4);  total (55/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 129 total aligned reads displayed.

TCATCGTCGCCGCCTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTGCGCCTGGGCCTGGCCGA  >  NC_002947/250763‑251057
                                                                                                                                                     |                                                                                                                                                 
tCATCGTCGCCGCCTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCgg                                                                                                                                                   <  1:2084450/150‑1 (MQ=255)
tCATCGTCGCCGCCTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCgg                                                                                                                                                   <  2:1360870/150‑1 (MQ=255)
  aTCGTCGCCGCCTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGgc                                                                                                                                                 >  1:167302/1‑150 (MQ=255)
       cgccgcCTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGcgccgc                                                                                                                                            <  2:1883567/150‑1 (MQ=255)
       cgccgcCTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCCCGcgccgc                                                                                                                                            <  1:2046187/150‑1 (MQ=255)
         ccgccTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTg                                                                                                                                          <  1:2417080/150‑1 (MQ=255)
          cgccTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTgg                                                                                                                                         <  1:2220997/150‑1 (MQ=255)
           gccTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGAATAATTTAGGTTTTTCCCGACAGCCATCGCTGCAATATGGCTGGCGTAGAGAATTTAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTggg                                                                                                                                        <  1:767863/150‑1 (MQ=255)
            ccTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTgggg                                                                                                                                       >  2:2035265/1‑150 (MQ=255)
              tttGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCACTGGGGGc                                                                                                                                     >  1:1284136/1‑150 (MQ=255)
                   cTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGt                                                                                                                                <  1:2263814/150‑1 (MQ=255)
                    tttATCGTCACTGGGGGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTa                                                                                                                               >  1:2050211/1‑150 (MQ=255)
                    tttATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCAGCTGGGAGcgcggcg                                                                                                                               >  2:391912/1‑144 (MQ=255)
                        tCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTg                                                                                                                           >  2:2134527/1‑150 (MQ=255)
                           tCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCg                                                                                                                        >  2:599443/1‑150 (MQ=255)
                            cACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGa                                                                                                                       >  2:217339/1‑150 (MQ=255)
                            cACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGa                                                                                                                       <  1:674343/150‑1 (MQ=255)
                              cTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACg                                                                                                                     >  1:2195997/1‑150 (MQ=255)
                               tGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGTACTTAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACgg                                                                                                                    >  1:982909/1‑150 (MQ=255)
                               tGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACgg                                                                                                                    <  2:1671774/150‑1 (MQ=255)
                                    gTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGgcgc                                                                                                                                               >  1:1113159/1‑118 (MQ=255)
                                       ccAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTgggg                                                                                                                                       <  2:773730/123‑1 (MQ=255)
                                       ccAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTgggg                                                                                                                                       >  1:773730/1‑123 (MQ=255)
                                        cAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGc                                                                                                           <  1:1877014/150‑1 (MQ=255)
                                          tgggcgccgccacgccggtgcggatcggcggcgcgcAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGgcgc                                                                                                                                               >  2:194942‑M2/37‑112 (MQ=255)
                                          tgggcgccgccacgccggtgcggatcggcggcgcgcAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGgcgc                                                                                                                                               <  1:194942‑M2/76‑1 (MQ=255)
                                          ggCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTg                                                                                                         <  2:1791298/150‑1 (MQ=255)
                                            cgcgcCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTGTAGAATATAGGTGCCGTCGCGGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGga                                                                                                       <  1:2122968/150‑1 (MQ=255)
                                             gcgcCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAAAGTAGCTGGCGACGGGTACGGTGCTGgag                                                                                                      >  1:603048/1‑150 (MQ=255)
                                                 cgccgATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGaa                                                                                                  <  1:1093175/150‑1 (MQ=255)
                                                 cgccgATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGaa                                                                                                  <  2:1184666/150‑1 (MQ=255)
                                                     gATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCggg                                                                                                                                                  <  2:2210102/98‑1 (MQ=255)
                                                     gATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTa                                                                                                                               <  1:311689/117‑1 (MQ=255)
                                                     gATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTa                                                                                                                               >  2:311689/1‑117 (MQ=255)
                                                      aTCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGcc                                                                                             >  2:1292387/1‑150 (MQ=255)
                                                      aTCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGcc                                                                                             <  2:815192/150‑1 (MQ=255)
                                                      aTCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGcc                                                                                             <  2:2195997/150‑1 (MQ=255)
                                                      aTCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGcc                                                                                             <  2:1449903/150‑1 (MQ=255)
                                                      aTCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGACATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGCGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGcc                                                                                             <  2:1114434/150‑1 (MQ=255)
                                                       tCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCt                                                                                            <  1:737656/150‑1 (MQ=255)
                                                       tCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCt                                                                                            <  2:1129313/150‑1 (MQ=255)
                                                        ccGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCtt                                                                                           <  2:1425116/150‑1 (MQ=255)
                                                        ccGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCtt                                                                                           <  2:2050211/150‑1 (MQ=255)
                                                          gCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGa                                                                                         >  1:1571319/1‑150 (MQ=255)
                                                          gCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGa                                                                                         >  2:195587/1‑150 (MQ=255)
                                                           cACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGat                                                                                        <  1:2269304/150‑1 (MQ=255)
                                                            aCCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTa                                                                                                                               <  1:2398492/110‑1 (MQ=255)
                                                            aCCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTa                                                                                                                               >  2:2398492/1‑110 (MQ=255)
                                                                gCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAAtt                                                                                                <  2:2192984/137‑1 (MQ=255)
                                                                gCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAAtt                                                                                                >  1:2192984/1‑137 (MQ=255)
                                                                  gTGGCGGCGCCCAGGTTGGTCCAGGGTACCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTgg                                                                                 <  1:262/150‑1 (MQ=255)
                                                                      cggcgCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCa                                                                             <  1:458137/150‑1 (MQ=255)
                                                                       ggcgCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCAt                                                                            <  1:2312171/150‑1 (MQ=255)
                                                                        gcgcCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTTCTGGAGAGAATTGCCTTGATCATCTGGCGCATg                                                                           <  1:391912/150‑1 (MQ=255)
                                                                         cgcCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGATAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATgg                                                                          <  1:1455149/150‑1 (MQ=255)
                                                                         cgcCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATgg                                                                          <  1:2010207/150‑1 (MQ=255)
                                                                           cccAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCa                                                                        >  1:1549059/1‑150 (MQ=255)
                                                                            ccAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTgggg                                                                                                                                       <  1:1835571/86‑1 (MQ=255)
                                                                            ccAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGcc                                                                                             >  1:2097975/1‑128 (MQ=255)
                                                                            ccAGGGTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCAc                                                                       >  2:1806974/1‑150 (MQ=255)
                                                                             cAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGcc                                                                                             >  2:2402571‑M2/2‑127 (MQ=255)
                                                                             cAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGcc                                                                                             <  1:2402571‑M2/126‑1 (MQ=255)
                                                                              aGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTg                                                                     >  2:1629048/1‑150 (MQ=255)
                                                                               ggTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGt                                                                    <  2:1909515/150‑1 (MQ=255)
                                                                                     tCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGcc                                                                                             <  2:502308/119‑1 (MQ=255)
                                                                                       ccGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAgcccg                                                            >  2:1743800/1‑150 (MQ=255)
                                                                                         gggCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAgcccgcc                                                          >  2:1035191/1‑150 (MQ=255)
                                                                                               cATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGt                                                    >  2:2137101/1‑150 (MQ=255)
                                                                                                 tCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATgg                                                                          <  1:2133046/126‑1 (MQ=255)
                                                                                                 tCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATgg                                                                          >  2:2133046/1‑126 (MQ=255)
                                                                                                   gCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATcc                                                <  1:1629048/150‑1 (MQ=255)
                                                                                                     tGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCggg                                                      <  1:60282/142‑1 (MQ=255)
                                                                                                     tGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCggg                                                      >  2:60282/1‑142 (MQ=255)
                                                                                                       cAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGc                                            <  1:1292387/150‑1 (MQ=255)
                                                                                                             ggctggcGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATaa                                      <  1:1296712/150‑1 (MQ=255)
                                                                                                              gctggcGTAGAGACTATAGGTCAGGGCTGTGGTTTCTCGGTCGCCGCTGGGGGCAACGTAGCTGGCGACGGGGACGGTGCTGGAGAGAATTGCCTTGATCAGCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAac                                     <  1:643314/150‑1 (MQ=255)
                                                                                                              gctggcGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAac                                     >  1:1667360/1‑150 (MQ=255)
                                                                                                               ctggcGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAacc                                    <  1:217339/150‑1 (MQ=255)
                                                                                                                tggcGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCa                                                                        >  1:1222342/1‑113 (MQ=255)
                                                                                                                tggcGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCa                                                                        <  2:1222342/113‑1 (MQ=255)
                                                                                                                tggcGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAacca                                   >  1:696454/1‑150 (MQ=255)
                                                                                                                     tAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGt                              >  1:707038/1‑150 (MQ=255)
                                                                                                                            tataGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAg                                                                >  1:495862/1‑109 (MQ=255)
                                                                                                                            tataGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAg                                                                <  2:495862/109‑1 (MQ=255)
                                                                                                                            tataGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTgcgcgc                       <  2:1816083/150‑1 (MQ=255)
                                                                                                                             ataGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGgcg                                                                               >  2:1838832/1‑93 (MQ=255)
                                                                                                                             ataGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGgcg                                                                               <  1:1838832/93‑1 (MQ=255)
                                                                                                                               aGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGGAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCtt                    >  1:788924/1‑150 (MQ=255)
                                                                                                                               aGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCtt                    <  2:1102862/150‑1 (MQ=255)
                                                                                                                                  tGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATCGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTgcg                 >  2:1148625/1‑150 (MQ=255)
                                                                                                                                      ggCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTGCgcctg             >  1:945233/1‑150 (MQ=255)
                                                                                                                                       gcgcTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTGCgcctgg            <  1:1371685/150‑1 (MQ=255)
                                                                                                                                        cgcTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTGCgcctggg           <  2:603048/150‑1 (MQ=255)
                                                                                                                                         cccagcggcgCCCGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTGCgcctgg            >  1:1135336/5‑148 (MQ=255)
                                                                                                                                          ccagcggcgCCCGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTGCgcctgggcc         >  2:1125577/4‑150 (MQ=255)
                                                                                                                                           cagcggcgCCCGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTGCgcctgg            <  2:422839/144‑1 (MQ=255)
                                                                                                                                           cagcggcgCCCGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTGCgcctgg            >  1:422839/3‑146 (MQ=255)
                                                                                                                                               ggcgCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTGCGCCTGGGCCTGGcc    >  1:2434445/1‑150 (MQ=255)
                                                                                                                                                 cgcgGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCACCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTGCGCCTGGGCCTGGCCGa  <  2:1818439/150‑1 (MQ=255)
                                                                                                                                                     ggCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAac                                     <  2:858400/111‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TCATCGTCGCCGCCTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTGCGCCTGGGCCTGGCCGA  >  NC_002947/250763‑251057

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