Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 250,911 G→C 5.7% R98G (CGC→GGC)  nfeD ← membrane bound peptidase NefD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947250,9110GC 5.7% 390.6 / 12.2 122R98G (CGC→GGC) nfeDmembrane bound peptidase NefD
Reads supporting (aligned to +/- strand):  ref base G (51/64);  new base C (4/3);  total (55/67)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

Only 100 of 128 total aligned reads displayed.

TTCATCGTCGCCGCCTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTGCGCCTGGGCCTGGCCGA  >  NC_002947/250762‑251057
                                                                                                                                                     |                                                                                                                                                  
ttCATCGTCGCCGCCTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGcgcg                                                                                                                                                    <  1:246502/150‑1 (MQ=255)
 tCATCGTCGCCGCCTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCgg                                                                                                                                                   <  2:1360870/150‑1 (MQ=255)
   aTCGTCGCCGCCTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGgc                                                                                                                                                 >  1:167302/1‑150 (MQ=255)
        cgccgcCTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCCCGcgccgc                                                                                                                                            <  1:2046187/150‑1 (MQ=255)
          ccgccTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTg                                                                                                                                          <  1:2417080/150‑1 (MQ=255)
           cgccTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTgg                                                                                                                                         <  1:669338/150‑1 (MQ=255)
            gccTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGAATAATTTAGGTTTTTCCCGACAGCCATCGCTGCAATATGGCTGGCGTAGAGAATTTAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTggg                                                                                                                                        <  1:767863/150‑1 (MQ=255)
              cTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTgg                                                                                                                                         >  2:2137649/1‑147 (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)
                               cTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACg                                                                                                                     >  1:2195997/1‑150 (MQ=255)
                                tGGGCGTGCCAGGCGCGCCGCCGATCCGCACCACCGGGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGAAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACgg                                                                                                                    >  1:1773325/1‑150 (MQ=255)
                                     gTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGgcgc                                                                                                                                               <  2:1113159/118‑1 (MQ=255)
                                     gTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGgcgc                                                                                                                                               >  1:1113159/1‑118 (MQ=255)
                                        ccAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTgggg                                                                                                                                       >  1:773730/1‑123 (MQ=255)
                                        ccAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTgggg                                                                                                                                       <  2:773730/123‑1 (MQ=255)
                                         cAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGc                                                                                                           <  1:1877014/150‑1 (MQ=255)
                                           tgggcgccgccacgccggtgcggatcggcggcgcgcAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGgcgc                                                                                                                                               <  1:194942‑M2/76‑1 (MQ=255)
                                           tgggcgccgccacgccggtgcggatcggcggcgcgcAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGgcgc                                                                                                                                               >  2:194942‑M2/37‑112 (MQ=255)
                                           ggCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTg                                                                                                         <  2:390469/150‑1 (MQ=255)
                                           ggCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTg                                                                                                         <  2:1791298/150‑1 (MQ=255)
                                             cgcgcCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTGTAGAATATAGGTGCCGTCGCGGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGga                                                                                                       <  1:2122968/150‑1 (MQ=255)
                                              gcgcCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAAAGTAGCTGGCGACGGGTACGGTGCTGgag                                                                                                      >  1:603048/1‑150 (MQ=255)
                                                  cgccgATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGCCGGGTACGATGCGGGAGAGaa                                                                                                  >  2:1320334/1‑150 (MQ=255)
                                                  cgccgATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGaa                                                                                                  <  1:1093175/150‑1 (MQ=255)
                                                     cgATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTg                                                                                               >  1:2343427/1‑150 (MQ=255)
                                                      gATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCggg                                                                                                                                                  <  2:2210102/98‑1 (MQ=255)
                                                      gATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCggg                                                                                                                                                  >  1:2210102/1‑98 (MQ=255)
                                                      gATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTa                                                                                                                               <  1:311689/117‑1 (MQ=255)
                                                      gATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGc                                                                                              <  2:568752/150‑1 (MQ=255)
                                                       aTCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGcc                                                                                             <  2:1449903/150‑1 (MQ=255)
                                                       aTCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGcc                                                                                             >  2:1292387/1‑150 (MQ=255)
                                                       aTCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGcc                                                                                             <  2:815192/150‑1 (MQ=255)
                                                       aTCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGcc                                                                                             <  1:345893/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:2050211/150‑1 (MQ=255)
                                                           gCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGa                                                                                         >  2:195587/1‑150 (MQ=255)
                                                           gCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGa                                                                                         >  1:1571319/1‑150 (MQ=255)
                                                            cACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGat                                                                                        <  1:2269304/150‑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)
                                                                    tGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGgc                                                                                <  1:1320334/150‑1 (MQ=255)
                                                                       cggcgCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCa                                                                             <  1:458137/150‑1 (MQ=255)
                                                                        ggcgCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCAt                                                                            <  1:2312171/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                                                                                                                                       >  2:1835571/1‑86 (MQ=255)
                                                                             ccAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGcc                                                                                             <  2:2097975/128‑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)
                                                                                   tGGTCCCGGGCGCCATCGCTGACATATGGCTGGCGTAGAGCATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCa                                                                 <  2:2343427/150‑1 (MQ=255)
                                                                                    ggTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAg                                                                <  1:2055072/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)
                                                                                                  tCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATgg                                                                          >  2:2133046/1‑126 (MQ=255)
                                                                                                  tCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATgg                                                                          <  1:2133046/126‑1 (MQ=255)
                                                                                                    gCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATcc                                                <  1:1629048/150‑1 (MQ=255)
                                                                                                      tGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCggg                                                      <  1:60282/142‑1 (MQ=255)
                                                                                                      tGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCggg                                                      >  2:60282/1‑142 (MQ=255)
                                                                                                       gCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATg                                             <  2:1767266/150‑1 (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                                                                        <  2:1222342/113‑1 (MQ=255)
                                                                                                                 tggcGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCa                                                                        >  1:1222342/1‑113 (MQ=255)
                                                                                                                 tggcGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAacca                                   >  1:696454/1‑150 (MQ=255)
                                                                                                                      tAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGt                              >  2:1921465/1‑150 (MQ=255)
                                                                                                                      tAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGt                              >  1:707038/1‑150 (MQ=255)
                                                                                                                             tataGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAg                                                                >  1:495862/1‑109 (MQ=255)
                                                                                                                              ataGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGgcg                                                                               <  1:1838832/93‑1 (MQ=255)
                                                                                                                              ataGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGgcg                                                                               >  2:1838832/1‑93 (MQ=255)
                                                                                                                                aGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGGAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCtt                    >  1:788924/1‑150 (MQ=255)
                                                                                                                                aGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCtt                    <  2:1102862/150‑1 (MQ=255)
                                                                                                                                aGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCtt                    >  2:668368/1‑150 (MQ=255)
                                                                                                                                  gTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTgc                  >  1:2424166/1‑150 (MQ=255)
                                                                                                                                   tGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATCGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTgcg                 >  2:1148625/1‑150 (MQ=255)
                                                                                                                                       ggCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTGCgcctg             >  1:945233/1‑150 (MQ=255)
                                                                                                                                          cccagcggcgCCCGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTGCgcctgg            >  1:1135336/5‑148 (MQ=255)
                                                                                                                                          cccagcggcgCCCGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTGCgcctgg            <  2:1135336/144‑1 (MQ=255)
                                                                                                                                           cTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTGCgcctgggcc         <  2:1571319/150‑1 (MQ=255)
                                                                                                                                            cagcggcgCCCGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTGCgcctgg            <  2:422839/144‑1 (MQ=255)
                                                                                                                                            cagcggcgCCCGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTGCgcctgg            >  1:422839/3‑146 (MQ=255)
                                                                                                                                                  cgcgGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCACCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTGCGCCTGGGCCTGGCCGa  <  2:1818439/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TTCATCGTCGCCGCCTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTGCGCCTGGGCCTGGCCGA  >  NC_002947/250762‑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: