Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,353,340 C→G 7.2% P141R (CCG→CGG)  yrpB → 2‑nitropropane dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,353,3400CG 7.2% 263.5 / 10.9 83P141R (CCG→CGG) yrpB2‑nitropropane dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (38/39);  new base G (2/4);  total (40/43)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.81e-01

CAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCATCGCGGCATGT  >  NC_002947/4353195‑4353489
                                                                                                                                                  |                                                                                                                                                     
cAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAg                                                                                                                                                   <  2:404445/150‑1 (MQ=255)
 aGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTGACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAgg                                                                                                                                                  >  1:934244/1‑150 (MQ=255)
 aGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAgg                                                                                                                                                  >  1:525752/1‑150 (MQ=255)
 aGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGATTCCATTTCGGCCTGCCGCAgg                                                                                                                                                  >  1:778576/1‑150 (MQ=255)
   gCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGcc                                                                                                                                                <  2:525752/150‑1 (MQ=255)
      cTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGaa                                                                                                                                             <  1:1582566/150‑1 (MQ=255)
       ttGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAAc                                                                                                                                            <  2:129565/150‑1 (MQ=255)
        tGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCGGCAGGCCGaaat                                                                                                                                           <  2:1579846/150‑3 (MQ=255)
        tGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCAACCAACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAGTTCGGCCTGCGGCAGGCCGaaat                                                                                                                                           <  1:120500/150‑3 (MQ=255)
           gCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCAGTTCGGCCTGCGGCAGGCCGaaatgga                                                                                                                                        <  2:704527/150‑6 (MQ=255)
            ccTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAActtctg                                                                                                                                       >  1:1860662/1‑150 (MQ=255)
            ccTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAActgctg                                                                                                                                       >  2:1514174/1‑150 (MQ=255)
              ttattaCAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCa                                                                                                                                     >  2:897103/1‑150 (MQ=255)
                  taCAGCGAAGTGGGCGCCGCTTTCAGGCCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCgg                                                                                                                                 <  2:934244/150‑1 (MQ=255)
                  taCAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCgg                                                                                                                                 <  1:882167/150‑1 (MQ=255)
                   aCAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCggg                                                                                                                                >  1:840748/1‑150 (MQ=255)
                   aCAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCggg                                                                                                                                <  2:1066001/150‑1 (MQ=255)
                    cAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCgggg                                                                                                                               >  2:1470119/1‑149 (MQ=255)
                     aGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTa                                                                                                                              >  2:1546652/1‑150 (MQ=255)
                      gCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTaa                                                                                                                             <  1:1546652/150‑1 (MQ=255)
                          aGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGc                                                                                                                         <  2:342056/150‑1 (MQ=255)
                               gcgcCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCg                                                                                                                    >  1:1940453/1‑150 (MQ=255)
                                  ccGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTg                                                                                                                 <  2:346588/150‑1 (MQ=255)
                                      tttCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCaa                                                                                                             <  2:1023444/150‑1 (MQ=255)
                                       ttCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAg                                                                                                            <  2:1489497/150‑1 (MQ=255)
                                            ggCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAgg                                                                                                                          >  2:1908529/1‑131 (MQ=255)
                                            ggCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAgg                                                                                                                          <  1:1908529/131‑1 (MQ=255)
                                            ggCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACt                                                                                                       <  1:1280534/150‑1 (MQ=255)
                                                 ccACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCa                                                                                                                                                    <  1:112741/100‑1 (MQ=255)
                                                 ccACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCa                                                                                                                                                    >  2:112741/1‑100 (MQ=255)
                                                 ccACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCa                                                                                                  >  2:1219618/1‑150 (MQ=255)
                                                         gtgtCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGt                                                                                                                  >  1:1696289/1‑126 (MQ=255)
                                                         gtgtCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGt                                                                                                                  <  2:1696289/126‑1 (MQ=255)
                                                          tgtCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGa                                                                                         <  1:269417/150‑1 (MQ=255)
                                                            tCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGAcc                                                                                       >  2:1734600/1‑150 (MQ=255)
                                                               aaCCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTa                                                                                                         <  1:44769/129‑1 (MQ=255)
                                                               aaCCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTa                                                                                                         >  2:44769/1‑129 (MQ=255)
                                                                aCCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGAGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTgg                                                                                   >  2:316324/1‑150 (MQ=255)
                                                                  cGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCgg                                                                                                                                 >  2:892972/1‑102 (MQ=255)
                                                                  cGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCgg                                                                                                                                 <  1:892972/102‑1 (MQ=255)
                                                                     gcgcCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGCGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCCGCGGTGCCAAGGTACTGTCCCGCGCCACGACCGGGgaagaa                                                                              >  1:768530/1‑150 (MQ=255)
                                                                     gcgcCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGgaagaa                                                                              <  2:840748/150‑1 (MQ=255)
                                                                       gcCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTa                                                                                                         >  2:955168/1‑121 (MQ=255)
                                                                       gcCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTa                                                                                                         <  1:955168/121‑1 (MQ=255)
                                                                       gcCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGc                                                                            >  1:372469/1‑150 (MQ=255)
                                                                              gacgaGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCctgg                                                                     <  1:1601495/150‑1 (MQ=255)
                                                                               acgaGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCctggc                                                                    >  1:1356750/1‑150 (MQ=255)
                                                                                 gaGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCctggctg                                                                  >  1:1921568/1‑150 (MQ=255)
                                                                                   gCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGa                                                                <  1:1221863/150‑1 (MQ=255)
                                                                                          tgcctgcTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAAt                                                         <  1:1219618/150‑1 (MQ=255)
                                                                                            cctgcTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATgg                                                       >  1:1461162/1‑150 (MQ=255)
                                                                                                   gTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTg                                                                                    <  2:286852/114‑1 (MQ=255)
                                                                                                   gTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTg                                                                                    >  1:286852/1‑114 (MQ=255)
                                                                                                    tCGAAGCCTTGCGCCCGGAGGTGGTGAGTTTCCATTTCGGCCTCACCATAAACCAGCTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgcgatg                                               <  1:1088671/150‑1 (MQ=255)
                                                                                                      gAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTg                                                                                                                 >  2:1609354/1‑82 (MQ=255)
                                                                                                      gAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTg                                                                                                                 <  1:1609354/82‑1 (MQ=255)
                                                                                                       aaGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTTCAGCGCCAAGTCCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgtgatgcga                                            <  1:728168/150‑1 (MQ=255)
                                                                                                       aaGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGgaag                                                                                <  1:108258/114‑1 (MQ=255)
                                                                                                       aaGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGgaag                                                                                >  2:108258/1‑114 (MQ=255)
                                                                                                                gcCCGGCGATGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCCGCGGGAAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCCCCACGACCCTGGCAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCa                                   >  2:431547/1‑150 (MQ=255)
                                                                                                                    ggAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATggg                               <  2:1548272/150‑1 (MQ=255)
                                                                                                                    ggAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATggg                               <  2:1461162/150‑1 (MQ=255)
                                                                                                                      aGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGt                             <  1:1777274/150‑1 (MQ=255)
                                                                                                                      aGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGt                             <  1:897103/150‑1 (MQ=255)
                                                                                                                           gtgAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGatc                                          <  1:1594617/132‑1 (MQ=255)
                                                                                                                           gtgAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGatc                                          >  2:1594617/1‑132 (MQ=255)
                                                                                                                                cTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTg                                                                                    >  2:1921000/1‑85 (MQ=255)
                                                                                                                                cTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTg                                                                                    <  1:1921000/85‑1 (MQ=255)
                                                                                                                                cTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGgaaga                                                                               >  1:196031/1‑90 (MQ=255)
                                                                                                                                cTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGgaaga                                                                               <  2:196031/90‑1 (MQ=255)
                                                                                                                                 ttCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGt                  >  1:77620/1‑150 (MQ=255)
                                                                                                                                 ttCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGt                  >  1:1432897/1‑150 (MQ=255)
                                                                                                                                  tCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTg                 >  2:1113829/1‑150 (MQ=255)
                                                                                                                                  cccATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTg                 <  1:1297044/149‑1 (MQ=255)
                                                                                                                                   gcagttCGGCCTGCGGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGaagaag                                                                             >  2:411828/5‑89 (MQ=255)
                                                                                                                                   gcagttCGGCCTGCGGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGaagaag                                                                             <  1:411828/85‑1 (MQ=255)
                                                                                                                                     aTTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATgg                                >  2:487409/1‑132 (MQ=255)
                                                                                                                                     aTTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATgg                                <  1:487409/132‑1 (MQ=255)
                                                                                                                                       ttCGGCCTGCGGCAGGCCGAAATGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCAt            >  1:1210488/1‑150 (MQ=255)
                                                                                                                                         cGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCATcg          >  2:868920/1‑150 (MQ=255)
                                                                                                                                           gCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgcgatgc                                              >  2:493979/1‑112 (MQ=255)
                                                                                                                                           gCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCtgcgatgc                                              <  1:493979/112‑1 (MQ=255)
                                                                                                                                             cTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCATCGCGGc      >  1:204683/1‑150 (MQ=255)
                                                                                                                                                  cGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCATCGCGGCATGt  >  2:1819109/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CAGGCACTTGAGCCTTATTACAGCGAAGTGGGCGCCGATTTCACGGCGCCCACGCCGGTGTCCAACCGTGCGCCTTTCGACGAGCAGAGCTGCCTGCTGGTCGAAGCGTTGCGCCCGGAGGTGGTGAGCTTCCATTTCGGCCTGCCGCAGGCCGAACTGCTGCAGCGGGTAAAGGCCAGCGGTGCCAAGGTACTGTCCAGCGCCACGACCGTGGAAGAAGCGGCCTGGCTGGAGCGCAATGGCTGCGATGCGATCATCGCCATGGGGTATGAAGCCGGTGGCCATCGCGGCATGT  >  NC_002947/4353195‑4353489

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: