Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,633,785 C→A 18.4% intergenic (‑40/‑12) yjjX ← / → ytjC non‑canonical purine NTP phosphatase, ITPase/XTPase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,633,7850CA18.4% 598.5 / 114.2 216intergenic (‑40/‑12)yjjX/ytjCnon‑canonical purine NTP phosphatase, ITPase/XTPase/phosphatase
Reads supporting (aligned to +/- strand):  ref base C (70/106);  new base A (19/21);  total (89/127)
Fisher's exact test for biased strand distribution p-value = 3.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 224 total aligned reads displayed.

GGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAACCCGTGCCAAAGAGCTT  >  NC_000913/4633638‑4633934
                                                                                                                                                   |                                                                                                                                                     
ggCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACgg                                                                                                                                                     <  1:1006764/150‑1 (MQ=255)
ggCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACgg                                                                                                                                                     <  1:2955676/150‑1 (MQ=255)
ggCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACgg                                                                                                                                                     <  1:2723549/150‑1 (MQ=255)
   aaCGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGaaa                                                                                                                                                  >  2:174245/1‑150 (MQ=255)
     cGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGaaaaa                                                                                                                                                >  1:700673/1‑150 (MQ=255)
      ggATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGaaaaaa                                                                                                                                               >  2:88150/1‑150 (MQ=255)
      ggATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGaaaaaa                                                                                                                                               >  2:2772751/1‑150 (MQ=255)
       gATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGaaaaaaa                                                                                                                                              <  1:268182/150‑1 (MQ=255)
        aTGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAg                                                                                                                                             <  2:1875036/150‑1 (MQ=255)
        aTGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAg                                                                                                                                             <  2:1445195/150‑1 (MQ=255)
         tGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGc                                                                                                                                            <  2:2064642/150‑1 (MQ=255)
         tGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGc                                                                                                                                            <  1:610316/150‑1 (MQ=255)
         tGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGc                                                                                                                                            <  1:193402/150‑1 (MQ=255)
          gCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCa                                                                                                                                           <  2:2399124/150‑1 (MQ=255)
          gCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCa                                                                                                                                           <  2:2361187/150‑1 (MQ=255)
            aaTATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATg                                                                                                                                         >  1:138502/1‑150 (MQ=255)
             atatGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGt                                                                                                                                        >  1:2156353/1‑150 (MQ=255)
             atatGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGAACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGt                                                                                                                                        >  1:2708706/1‑150 (MQ=255)
              tatGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGtt                                                                                                                                       >  1:1094758/1‑150 (MQ=255)
                  gtaGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAg                                                                                                                                   <  2:2555406/148‑1 (MQ=255)
                   taGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAgg                                                                                                                                  >  1:1314503/2‑150 (MQ=255)
                   cAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAgg                                                                                                                                  >  1:992907/1‑150 (MQ=255)
                       aTCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGtata                                                                                                                              >  2:915528/1‑150 (MQ=255)
                        tCCTTCGCCGAAGATCTCGTGAAAGGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACATACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATAc                                                                                                                             <  2:1840862/150‑1 (MQ=255)
                          ctacGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCt                                                                                                                           >  1:9131/4‑150 (MQ=255)
                          cTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCt                                                                                                                           >  2:297128/1‑150 (MQ=255)
                              gCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGACAAAAAGCATGTTACAGGTATACCTAGTc                                                                                                                       >  2:2742110/1‑150 (MQ=255)
                              gCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATAc                                                                                                                             >  1:2871465/1‑144 (MQ=255)
                               ccGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTcc                                                                                                                      >  1:2437483/1‑150 (MQ=255)
                               ccGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTcc                                                                                                                      >  1:378714/1‑150 (MQ=255)
                                cGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCg                                                                                                                     <  1:1428824/150‑1 (MQ=255)
                                 gAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGc                                                                                                                    >  2:1197829/1‑150 (MQ=255)
                                 gAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGc                                                                                                                    <  2:67715/150‑1 (MQ=255)
                                 gAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGc                                                                                                                    >  2:2197897/1‑150 (MQ=255)
                                       ctcGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGt                                                                                                              <  1:1531861/150‑1 (MQ=255)
                                       ctcGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGt                                                                                                              <  1:1727663/150‑1 (MQ=255)
                                        tcGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTg                                                                                                             >  2:965644/1‑150 (MQ=255)
                                          gTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGaa                                                                                                           >  1:2126409/1‑150 (MQ=255)
                                               aTGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCa                                                                                                      <  2:2928361/150‑1 (MQ=255)
                                               aTGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCa                                                                                                      >  2:2081634/1‑150 (MQ=255)
                                                tGCCTGCAGAATGGCCTGAATTTTAGCGGGTTTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAg                                                                                                     <  2:612564/150‑1 (MQ=255)
                                                 gCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGt                                                                                                    <  1:726549/150‑1 (MQ=255)
                                                 gCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGt                                                                                                    <  2:1794728/150‑1 (MQ=255)
                                                    tGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGa                                                                                                 >  2:2877132/1‑150 (MQ=255)
                                                         aaTGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGcc                                                                                            <  2:1296746/150‑1 (MQ=255)
                                                           tGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGa                                                                                          <  1:2066753/150‑1 (MQ=255)
                                                            ggCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAg                                                                                         <  1:1572458/150‑1 (MQ=255)
                                                             gCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGAGGAACGCCGAGc                                                                                        >  2:195437/1‑150 (MQ=255)
                                                               cTGAATTTTAGCGGGATTGGTGGACGCACAGACAACTAGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGa                                                                                      >  1:463369/1‑150 (MQ=255)
                                                                   aTTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTa                                                                                  <  2:138232/150‑1 (MQ=255)
                                                                      ttAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTc                                                                               >  2:1331619/1‑150 (MQ=255)
                                                                         gCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAgg                                                                            >  2:1461738/1‑150 (MQ=255)
                                                                         gCGGGATTGGTGGTCGCACAGACAACTAGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAgg                                                                            <  2:2717076/150‑1 (MQ=255)
                                                                          cGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAggg                                                                           >  1:2612807/1‑150 (MQ=255)
                                                                           gggATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGc                                                                          <  2:1969842/150‑1 (MQ=255)
                                                                           gggATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGc                                                                          <  2:2165610/150‑1 (MQ=255)
                                                                           gggATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGc                                                                          >  2:761977/1‑150 (MQ=255)
                                                                             gATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCa                                                                        <  1:3036316/150‑1 (MQ=255)
                                                                               ttGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGt                                                                      <  2:834582/150‑1 (MQ=255)
                                                                                  gtggtCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTg                                                                   >  2:1363541/1‑150 (MQ=255)
                                                                                   tggtCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGa                                                                  <  2:1400169/150‑1 (MQ=255)
                                                                                   tggtCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGa                                                                  <  1:2823954/150‑1 (MQ=255)
                                                                                     gtCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACa                                                                <  1:532407/150‑1 (MQ=255)
                                                                                        gCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGcc                                                             <  1:973636/150‑1 (MQ=255)
                                                                                          acaGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCg                                                           <  2:2259202/150‑1 (MQ=255)
                                                                                          acaGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCg                                                           <  1:1153296/150‑1 (MQ=255)
                                                                                          acaGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCg                                                           <  2:333428/150‑1 (MQ=255)
                                                                                          acaGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGGATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCg                                                           <  1:815907/150‑1 (MQ=255)
                                                                                             gACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTg                                                        >  2:3055600/1‑150 (MQ=255)
                                                                                             gACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTg                                                        >  1:2917401/1‑150 (MQ=255)
                                                                                              aCAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGa                                                       >  1:1825632/1‑149 (MQ=255)
                                                                                                 aCTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCg                                                    <  1:2772751/150‑1 (MQ=255)
                                                                                                         cATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGt                                            <  2:331131/150‑1 (MQ=255)
                                                                                                         cATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGt                                            <  2:778683/150‑1 (MQ=255)
                                                                                                           tAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGa                                          <  2:3033853/150‑1 (MQ=255)
                                                                                                              tCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCa                                       >  2:2651836/1‑150 (MQ=255)
                                                                                                               cAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCaa                                      <  2:822720/150‑1 (MQ=255)
                                                                                                                 gCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACa                                    >  2:2004399/1‑150 (MQ=255)
                                                                                                                 gCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACa                                    <  1:567713/150‑1 (MQ=255)
                                                                                                                  cATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAg                                   <  1:880134/150‑1 (MQ=255)
                                                                                                                        tCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATg                             >  1:1692359/1‑150 (MQ=255)
                                                                                                                        tCAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATg                             >  1:1347215/1‑150 (MQ=255)
                                                                                                                         cAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGc                            <  2:2437483/150‑1 (MQ=255)
                                                                                                                         cAGAAAATTAACGTTACAGCAGTATAAGGAAAAAAAGCTTGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGc                            <  1:2926422/150‑1 (MQ=255)
                                                                                                                          aGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCa                           >  1:2413236/1‑150 (MQ=255)
                                                                                                                           gAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAg                          <  1:1371893/150‑1 (MQ=255)
                                                                                                                           gAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAg                          <  1:115846/150‑1 (MQ=255)
                                                                                                                             aaaTTAACGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGt                        <  2:57451/150‑1 (MQ=255)
                                                                                                                              aaTTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTg                       <  1:1378045/150‑1 (MQ=255)
                                                                                                                                  atcttttcagcacTATAAGGTATAAAAGCATGTTACAGGTATAGCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAAGTGGCGATGCAGGTGGCaa                   <  2:2015823/143‑1 (MQ=255)
                                                                                                                                  aaCGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAACCGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCaa                   <  2:1112542/150‑1 (MQ=255)
                                                                                                                                   aCGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAAc                  <  2:2156353/150‑1 (MQ=255)
                                                                                                                                   aCGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTc                                                                               <  2:1472285/89‑1 (MQ=255)
                                                                                                                                   aCGTTACAGCAGTATAAGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTc                                                                               >  1:1472285/1‑89 (MQ=255)
                                                                                                                                      ttACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAACCCg               >  2:1547273/1‑150 (MQ=255)
                                                                                                                                             agTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAACCCGTGCCaaa        <  1:1516251/150‑1 (MQ=255)
                                                                                                                                                 taCGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAACCCGTGCCAAAGAGc    >  1:1140013/1‑150 (MQ=255)
                                                                                                                                                 taCGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAACCCGTGCCAAAGAGc    >  1:1010106/1‑150 (MQ=255)
                                                                                                                                                   cGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAACCCGTGCCAAAGAGCtt  <  1:1062296/150‑1 (MQ=255)
                                                                                                                                                   cGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAACCCGTGCCAAAGAGCtt  <  1:1062137/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACAGGTATACCTAGTCCGCCACGGTGAAACGCAGTGGAACGCCGAGCGACGTATTCAGGGCCAGTCTGACAGCCCGCTGACCGCCAAAGGTGAGCAACAGGCGATGCAGGTGGCAACCCGTGCCAAAGAGCTT  >  NC_000913/4633638‑4633934

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: