Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,540,645 G→C 6.4% A649G (GCC→GGC) ‡ narZ ← nitrate reductase Z subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,540,6450GC 6.4% 311.0 / 9.5 93A649G (GCC→GGC) ‡narZnitrate reductase Z subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (43/44);  new base C (2/4);  total (45/48)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

GTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTGGCGCGGTGGTCGG  >  NC_000913/1540499‑1540794
                                                                                                                                                  |                                                                                                                                                     
gttCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTCGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCCGTTTACGCCCCAGCTGCGGCGCAGAAGCTAGCCAGCCCATACGTTCGgcgc                                                                                                                                                    >  2:978816/1‑150 (MQ=255)
    gCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGaaa                                                                                                                                                <  1:379816/150‑2 (MQ=255)
      aGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGt                                                                                                                                              >  2:413666/1‑150 (MQ=255)
         ccATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGCGCCGAACGtatg                                                                                                                                           <  1:751224/150‑4 (MQ=255)
         ccATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGCGCCGAACGtatg                                                                                                                                           <  1:1000913/150‑4 (MQ=255)
          cATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGCGCCGAACGtatgg                                                                                                                                          <  2:424486/150‑5 (MQ=255)
           aTACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGaaa                                                                                                                                         <  1:13511/150‑1 (MQ=255)
           aTACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGCGCCGAACGtatggg                                                                                                                                         <  1:674412/150‑6 (MQ=255)
            tACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAaat                                                                                                                                        <  1:775388/150‑1 (MQ=255)
             aCGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAaatc                                                                                                                                       <  2:305676/150‑1 (MQ=255)
                 aaaTCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAaatcaatc                                                                                                                                   <  1:1120292/150‑1 (MQ=255)
                  aaTCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGATGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAatcaatca                                                                                                                                  <  2:1535122/150‑1 (MQ=255)
                  aaTCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAatcaatca                                                                                                                                  >  1:591054/1‑150 (MQ=255)
                    tCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGa                                                                                                                                >  1:1496316/1‑150 (MQ=255)
                      gCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATg                                                                                                                              >  2:196603/1‑150 (MQ=255)
                      gCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATg                                                                                                                              >  1:690016/1‑150 (MQ=255)
                      gCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATg                                                                                                                              >  1:1564315/1‑150 (MQ=255)
                               cAGCGCCTGGGCCGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGta                                                                                                                     >  2:1005567/1‑150 (MQ=255)
                                 gcgcCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGtact                                                                                                                   <  1:1250931/150‑1 (MQ=255)
                                 gcgcCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCCTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGtact                                                                                                                   <  2:302349/150‑1 (MQ=255)
                                      tGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTg                                                                                                              >  2:408867/1‑150 (MQ=255)
                                       gggCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTgg                                                                                                             <  1:443009/150‑1 (MQ=255)
                                          cGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCAt                                                                                                          >  1:896357/1‑150 (MQ=255)
                                           ggTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCCGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATc                                                                                                         >  1:760267/1‑150 (MQ=255)
                                             tAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCgg                                                                                                       <  1:1005567/150‑1 (MQ=255)
                                                aTTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGa                                                                                                    >  2:478997/1‑150 (MQ=255)
                                                aTTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGa                                                                                                    >  1:546737/1‑150 (MQ=255)
                                                 ttCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAg                                                                                                   <  1:1462137/150‑1 (MQ=255)
                                                   cTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCg                                                                                                 <  2:1141946/150‑1 (MQ=255)
                                                   cTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCg                                                                                                 <  1:812738/150‑1 (MQ=255)
                                                     gtgGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGt                                                                                               <  2:874527/150‑1 (MQ=255)
                                                      tgGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTg                                                                                              <  2:389808/150‑1 (MQ=255)
                                                      tgGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTg                                                                                              <  1:168366/150‑1 (MQ=255)
                                                        ggggATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGaa                                                                                            <  2:562219/150‑1 (MQ=255)
                                                        ggggATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGaa                                                                                            >  2:1530848/1‑150 (MQ=255)
                                                         gggATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGaaa                                                                                           >  1:294495/1‑150 (MQ=255)
                                                         ggcaTAATCCGGCCTTGGCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTTCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTTTACTTACTGGCATCGGCGAGCGGTGaaa                                                                                           <  2:1135171/147‑1 (MQ=255)
                                                          ggATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAg                                                                                          >  2:1359427/1‑150 (MQ=255)
                                                          ggATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAg                                                                                          <  1:193698/150‑1 (MQ=255)
                                                            aTAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTa                                                                                        <  1:104186/150‑1 (MQ=255)
                                                             tAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTaa                                                                                       >  1:788428/1‑150 (MQ=255)
                                                             tAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTCGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTaa                                                                                       >  1:138615/1‑150 (MQ=255)
                                                                 ccGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTcc                                                                                   >  2:814500/1‑150 (MQ=255)
                                                                   ggCCTTGTCGTCTTCAGCTTTAATCGCGAGCGGGTTACGCGCCAGCTGCGGCTCAGAAGGTAGCCAGCCCATAGGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTg                                                                                 <  2:770739/150‑1 (MQ=255)
                                                                   ggCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTg                                                                                 <  1:408867/150‑1 (MQ=255)
                                                                    gCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTCGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGa                                                                                <  2:138615/150‑1 (MQ=255)
                                                                     ccTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAg                                                                               <  2:1222660/150‑1 (MQ=255)
                                                                      cTTGTCGGCTTCAGCTTTAATCCCGAGCGGGGTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGc                                                                              <  2:121039/150‑1 (MQ=255)
                                                                          tCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCagag                                                                          <  1:315479/150‑1 (MQ=255)
                                                                           cGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCagagc                                                                         >  2:1367538/1‑149 (MQ=255)
                                                                              cTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACtt                                                                      >  1:1304007/1‑150 (MQ=255)
                                                                              cTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACtt                                                                      >  1:840005/1‑150 (MQ=255)
                                                                               ttCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACttt                                                                     <  2:1564315/150‑1 (MQ=255)
                                                                                  aGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACttttta                                                                  >  2:1622311/1‑149 (MQ=255)
                                                                                  aGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTc                                                                  >  2:1304344/1‑150 (MQ=255)
                                                                                   gCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCa                                                                 <  1:814500/150‑1 (MQ=255)
                                                                                   gCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCa                                                                 >  1:372335/1‑150 (MQ=255)
                                                                                     tttAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATa                                                               <  2:690016/150‑1 (MQ=255)
                                                                                     tttAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATa                                                               >  2:73422/1‑150 (MQ=255)
                                                                                      ttcaTCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAg                                                              <  2:726072/147‑1 (MQ=255)
                                                                                        aaTCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAgcg                                                            <  2:111369/150‑1 (MQ=255)
                                                                                        aaTCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAgcg                                                            <  2:8765/150‑1 (MQ=255)
                                                                                           cccGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGATCAGAGACTTTTTCATAGCGCCa                                                         >  2:338044/1‑150 (MQ=255)
                                                                                           cccGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCa                                                         >  1:617673/1‑150 (MQ=255)
                                                                                            ccGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCAt                                                        >  2:348329/1‑150 (MQ=255)
                                                                                              gAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTg                                                      >  1:796852/1‑150 (MQ=255)
                                                                                                gCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGc                                                    <  2:546737/150‑1 (MQ=255)
                                                                                                  gggTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTg                                                  >  1:850602/1‑150 (MQ=255)
                                                                                                   ggTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTgg                                                 >  1:752900/1‑150 (MQ=255)
                                                                                                    gTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGa                                                >  1:118756/1‑150 (MQ=255)
                                                                                                          ccccAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGAtt                                          >  1:1230540/1‑150 (MQ=255)
                                                                                                          ccccAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGAtt                                          <  1:1320617/150‑1 (MQ=255)
                                                                                                            ccAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGt                                        <  1:1217377/150‑1 (MQ=255)
                                                                                                             cAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTa                                       >  1:1370077/1‑150 (MQ=255)
                                                                                                             cAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTa                                       >  2:1260135/1‑150 (MQ=255)
                                                                                                              cgCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTCCTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAg                                      <  2:1583574/149‑1 (MQ=255)
                                                                                                              aGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAg                                      <  1:648971/150‑1 (MQ=255)
                                                                                                               gCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGaggatcgg                                     >  1:333529/1‑143 (MQ=255)
                                                                                                                  gcggcgCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGaaaa                                  >  2:706253/1‑150 (MQ=255)
                                                                                                                        cAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAgg                            >  2:1627712/1‑150 (MQ=255)
                                                                                                                        cAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAgg                            >  1:514153/1‑150 (MQ=255)
                                                                                                                        cAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAgg                            >  2:1169131/1‑150 (MQ=255)
                                                                                                                          gAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTg                          <  2:788428/150‑1 (MQ=255)
                                                                                                                          gAAGGTAGCCAGCCCATACGATCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAACGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTg                          >  2:807934/1‑150 (MQ=255)
                                                                                                                            tggTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTATCTGCTGAGCAGAGTCTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGTGGTGCt                        <  2:37586/149‑1 (MQ=255)
                                                                                                                                 gccagccCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCa                   <  1:978816/150‑1 (MQ=255)
                                                                                                                                   cagccCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCAtt                 <  2:1577496/150‑1 (MQ=255)
                                                                                                                                    agccCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCAttt                <  2:840005/150‑1 (MQ=255)
                                                                                                                                     atttcaaCGTTCGCGCCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTg               >  1:1590279/7‑150 (MQ=255)
                                                                                                                                        tcaaCGTTCGCGCCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTGgcg            >  1:53923/4‑150 (MQ=255)
                                                                                                                                          tACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTGGcgcg          >  2:357846/1‑150 (MQ=255)
                                                                                                                                             gTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTGGCGCggtg       >  1:494179/1‑150 (MQ=255)
                                                                                                                                              ttCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGTCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTGGCGCggtgg      <  1:1304344/150‑1 (MQ=255)
                                                                                                                                              ttCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTGGCTCggtgg      <  2:796852/150‑1 (MQ=255)
                                                                                                                                                 ggCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTGGCGCGGTGGTCg   >  1:1359494/1‑150 (MQ=255)
                                                                                                                                                  gcgcGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTGGCGCGGTGGTCgg  <  1:1122053/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGTTACGCCCCAGCTGCGGCGCAGAAGGTAGCCAGCCCATACGTTCGGCGCGAACGTTGAAATCAATCAGATGACCGCTGTACTTACTGGCATCGGCGAGCGGTGAAAGTAACTCCTGAGCAGAGACTTTTTCATAGCGCCATTGGCTGGAATGATTGTAGAAAAACGAGGTGCTGTTCATTTGGCGCGGTGGTCGG  >  NC_000913/1540499‑1540794

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