Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,645,929 G→C 9.4% S490R (AGC→AGG) ‡ gcd ← quinoprotein glucose dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,645,9290GC 9.4% 409.5 / 26.2 128S490R (AGC→AGG) ‡gcdquinoprotein glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (56/60);  new base C (5/7);  total (61/67)
Fisher's exact test for biased strand distribution p-value = 7.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 136 total aligned reads displayed.

ACTGCGCCTTGTGGGGTGGGGATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTGGTCCGGGGTCTGGTT  >  NC_002947/1645780‑1646078
                                                                                                                                                     |                                                                                                                                                     
aCTGCGCCTTGTGGGGTGGGGATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCggctgg                                                                                                                                                       >  1:4883137/1‑150 (MQ=255)
  tGCGCCTTGTGGGGTGGGGATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCggctggct                                                                                                                                                     >  1:4963500/1‑150 (MQ=255)
          gtggggtgggGATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCAGCCAGCCGAccc                                                                                                                                             <  2:3884921/150‑3 (MQ=255)
          gtggggtgggGATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCAGCCAGCCGAccc                                                                                                                                             <  1:3661175/150‑3 (MQ=255)
            ggggtggggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTcc                                                                                                                                           <  1:3803577/150‑1 (MQ=255)
            ggggtggggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCAGCCAGCCGAccctg                                                                                                                                           <  1:5037813/150‑5 (MQ=255)
             gggtggggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCAGCCAGCCGAccctgg                                                                                                                                          <  1:971388/150‑6 (MQ=255)
              ggtggggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCAt                                                                                                                                         <  2:1980889/150‑1 (MQ=255)
               gtggggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATg                                                                                                                                        >  1:4715641/1‑150 (MQ=255)
               gtggggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATg                                                                                                                                        >  1:4715645/1‑150 (MQ=255)
                  gggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTcc                                                                                                                                     >  1:3838011/1‑150 (MQ=255)
                  gggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTcc                                                                                                                                     <  1:4015709/150‑1 (MQ=255)
                    gATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCca                                                                                                                                   <  2:3598540/150‑1 (MQ=255)
                     aTCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCcac                                                                                                                                  <  2:1806625/150‑1 (MQ=255)
                       ctcGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAg                                                                                                                                <  2:1759837/150‑1 (MQ=255)
                         cgcgGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGt                                                                                                                              <  1:2840612/150‑1 (MQ=255)
                          gcgGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTc                                                                                                                             >  2:3993369/1‑150 (MQ=255)
                          gcgGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTc                                                                                                                             >  1:4927125/1‑150 (MQ=255)
                          gcgGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTc                                                                                                                             >  1:1501962/1‑150 (MQ=255)
                          gcgGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTc                                                                                                                             >  1:4451852/1‑150 (MQ=255)
                          gcgGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCACCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACGGGTc                                                                                                                             >  2:2358727/1‑150 (MQ=255)
                            ggATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCgt                                                                                                                           >  1:2465280/1‑150 (MQ=255)
                                  gCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTggtggg                                                                                                                     >  2:2827634/1‑150 (MQ=255)
                                       cttGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAAc                                                                                                                <  1:3661071/149‑1 (MQ=255)
                                       aTTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAAc                                                                                                                <  2:3215771/150‑1 (MQ=255)
                                       aTTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAAc                                                                                                                <  1:4361246/150‑1 (MQ=255)
                                            cGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATa                                                                                                           <  2:1095309/150‑1 (MQ=255)
                                            cGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATa                                                                                                           <  2:1982982/150‑1 (MQ=255)
                                            cGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATa                                                                                                           <  1:961598/150‑1 (MQ=255)
                                            cGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATa                                                                                                           <  1:1228459/150‑1 (MQ=255)
                                             gTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAg                                                                                                          >  2:785665/1‑150 (MQ=255)
                                                 cGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTcc                                                                                                      <  1:4707040/150‑1 (MQ=255)
                                                  gTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCa                                                                                                     <  2:1871328/150‑1 (MQ=255)
                                                    cgcgGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGc                                                                                                   <  1:2892375/150‑1 (MQ=255)
                                                                  aCGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTg                                                                                     >  2:3376154/1‑150 (MQ=255)
                                                                      aCAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCa                                                                                 <  2:5012259/150‑1 (MQ=255)
                                                                      aCAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCa                                                                                 <  2:4749932/150‑1 (MQ=255)
                                                                       cAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCAGGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAg                                                                                <  1:1901400/150‑1 (MQ=255)
                                                                         ggCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGt                                                                              >  1:2321415/1‑150 (MQ=255)
                                                                          gCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTc                                                                             >  1:5034555/1‑150 (MQ=255)
                                                                            tGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCaa                                                                           <  1:3602870/150‑1 (MQ=255)
                                                                             gCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCATGCCGGTGGCCAGGTCAAg                                                                          <  2:1282285/150‑1 (MQ=255)
                                                                              ccTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAgg                                                                         >  1:1128105/1‑150 (MQ=255)
                                                                              ccTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAgg                                                                         <  2:203521/150‑1 (MQ=255)
                                                                               cttgcttgGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAggg                                                                        >  2:166602/1‑150 (MQ=255)
                                                                                ttgcttgGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGc                                                                       >  2:2375862/1‑150 (MQ=255)
                                                                                ttgcttgGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGc                                                                       >  2:4989604/1‑150 (MQ=255)
                                                                                  gcttgGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGa                                                                     >  1:78427/1‑150 (MQ=255)
                                                                                    ttgGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGAcc                                                                   <  1:1547302/150‑1 (MQ=255)
                                                                                    ttgGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGAcc                                                                   >  1:1972720/1‑150 (MQ=255)
                                                                                       gTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACg                                                                >  1:321689/1‑150 (MQ=255)
                                                                                        tCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGc                                                               <  1:1900197/150‑1 (MQ=255)
                                                                                          ggTACGATGCTCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCg                                                             <  1:131456/150‑1 (MQ=255)
                                                                                          ggTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCg                                                             >  2:1958406/1‑150 (MQ=255)
                                                                                          ggTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCg                                                             <  2:1265682/150‑1 (MQ=255)
                                                                                           gTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCgg                                                            >  1:546819/1‑150 (MQ=255)
                                                                                              cGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGgcgc                                                         <  1:3506211/150‑1 (MQ=255)
                                                                                              cGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGgcgc                                                         <  2:588232/150‑1 (MQ=255)
                                                                                               gatgatCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCt                                                        >  1:2296880/1‑150 (MQ=255)
                                                                                               gatgatCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCt                                                        >  2:1337031/1‑150 (MQ=255)
                                                                                               gatgatCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCt                                                        >  2:1138648/1‑150 (MQ=255)
                                                                                                 tgatCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGt                                                      >  1:4560035/1‑150 (MQ=255)
                                                                                                  gatCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTa                                                     <  1:2433395/150‑1 (MQ=255)
                                                                                                   atCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGc                                                               >  1:3630444/1‑139 (MQ=255)
                                                                                                   atCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTAc                                                    <  2:4848215/150‑1 (MQ=255)
                                                                                                    tCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCgtggt                                                                                                                        >  2:692216/1‑81 (MQ=255)
                                                                                                    tCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCgtggt                                                                                                                        <  1:692216/81‑1 (MQ=255)
                                                                                                    tCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACt                                                   >  1:4798984/1‑150 (MQ=255)
                                                                                                    tCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCAGGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACt                                                   >  2:1410180/1‑150 (MQ=255)
                                                                                                    tCGCCGGTTTCAAACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCAGTACt                                                   >  2:1524891/1‑150 (MQ=255)
                                                                                                        cGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTtctc                                               <  2:281186/150‑1 (MQ=255)
                                                                                                        cGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTtctc                                               <  2:281109/150‑1 (MQ=255)
                                                                                                         ggTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCg                                              <  2:4887128/150‑1 (MQ=255)
                                                                                                             tCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCGGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGgcgc                                          >  2:3257412/1‑150 (MQ=255)
                                                                                                               acacCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCgtggtg                                                                                                                       >  1:3349998/1‑71 (MQ=255)
                                                                                                               acacCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGccc                                        >  1:1491181/1‑150 (MQ=255)
                                                                                                                  ccATCGTCGGTTGTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCggg                                     <  2:4221266/150‑1 (MQ=255)
                                                                                                                   cATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCgggg                                    >  1:4008020/1‑150 (MQ=255)
                                                                                                                       gtcgGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCg                                >  1:3431988/1‑150 (MQ=255)
                                                                                                                            ttttCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCgg                           >  2:2179392/1‑150 (MQ=255)
                                                                                                                            ttttCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCgg                           >  2:2179338/1‑150 (MQ=255)
                                                                                                                            ttttCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCgg                           >  2:1094078/1‑150 (MQ=255)
                                                                                                                                 aGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGccc                      <  2:357846/150‑1 (MQ=255)
                                                                                                                                  ggTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGcccc                     <  2:3838011/150‑1 (MQ=255)
                                                                                                                                    tgtACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCAc                   <  2:4715645/150‑1 (MQ=255)
                                                                                                                                    tgtACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCAc                   <  2:4715641/150‑1 (MQ=255)
                                                                                                                                     gtACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACt                  <  1:3850882/150‑1 (MQ=255)
                                                                                                                                      tACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTg                 <  1:520221/150‑1 (MQ=255)
                                                                                                                                      tACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGGTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTg                 <  1:520410/150‑1 (MQ=255)
                                                                                                                                       atggacgTCGGCAGCCAGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTgg                >  2:1453193/7‑150 (MQ=255)
                                                                                                                                        ccAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTGGt               >  1:2192937/1‑150 (MQ=255)
                                                                                                                                        ccAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTGGt               >  2:1240066/1‑150 (MQ=255)
                                                                                                                                         ggacgTCGGCAGCCAGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTGGTc              >  1:3498913/5‑150 (MQ=255)
                                                                                                                                         ggacgTCGGCAGCCAGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTGGTc              >  1:3217553/5‑150 (MQ=255)
                                                                                                                                          gacgTCGGCAGCCAGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTGGTcc             >  2:3201364/4‑150 (MQ=255)
                                                                                                                                            ggTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTGGTCCgg           <  1:3096717/150‑1 (MQ=255)
                                                                                                                                             gTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTGGTCCggg          <  1:973684/150‑1 (MQ=255)
                                                                                                                                                    ggctgCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTGGTCCGGGGTCTGGt   <  2:4883137/150‑1 (MQ=255)
                                                                                                                                                     gctgCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTGGTCCGGGGTCTGGtt  <  2:4451852/150‑1 (MQ=255)
                                                                                                                                                     gctgCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTGGTCCGGGGTCTGGtt  <  1:357476/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACTGCGCCTTGTGGGGTGGGGATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTGGTCCGGGGTCTGGTT  >  NC_002947/1645780‑1646078

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: