Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,364,995 T→G 8.9% intergenic (+16/+23) nudI → / ← ais nucleoside triphosphatase/putative LPS core heptose(II)‑phosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,364,9950TG 8.9% 341.1 / 22.4 124intergenic (+16/+23)nudI/aisnucleoside triphosphatase/putative LPS core heptose(II)‑phosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base T (57/56);  new base G (4/7);  total (61/63)
Fisher's exact test for biased strand distribution p-value = 5.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.21e-01

Only 100 of 125 total aligned reads displayed.

GCGTTTCTGCCAACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAAT  >  NC_000913/2364846‑2365142
                                                                                                                                                     |                                                                                                                                                   
gCGTTTCTGCCAACCGAGAAGTGTAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAt                                                                                                                                                     >  2:1340407/1‑150 (MQ=255)
      cTGCCAACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTGCGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTgg                                                                                                                                               >  2:188012/1‑150 (MQ=255)
       tGCCAACCGAGAAGTGAAACTAAACGAAGAGTATCACGACTACGCGTGGGTAATACCTGAAGATCTGGTGCATTAAGATTTGAATGTCGCCACCCGAAAAACGTTACGTCTGAAAGGTCTTCTGTAATAATGATGACAGCCATTACTGGc                                                                                                                                              <  1:1117187/150‑1 (MQ=255)
        gCCAACCGAGAAGTGAAAATAAACGAAGAGTTTCATGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCt                                                                                                                                             <  2:1272167/150‑1 (MQ=255)
        gCCAACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCt                                                                                                                                             <  2:466785/150‑1 (MQ=255)
         ccACCCGAGTAGTGAAATTAATCGAAGAGTTTCAGGTCTACGCGTGGGTAAATCCTGAAGATCTTGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTg                                                                                                                                            <  1:182065/150‑1 (MQ=255)
              cGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATCACAGCCAGTAATGGCTGTCatc                                                                                                                                       <  1:1492785/150‑3 (MQ=255)
              cGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATCACAGCCAGTAATGGCTGTCatc                                                                                                                                       <  1:1197900/150‑3 (MQ=255)
                 gaAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTg                                                                                                                                    >  1:1646214/1‑150 (MQ=255)
                 gaAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTg                                                                                                                                    >  2:1266082/1‑150 (MQ=255)
                   aGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGaa                                                                                                                                  >  1:810381/1‑150 (MQ=255)
                   aGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGaa                                                                                                                                  >  1:1218772/1‑150 (MQ=255)
                    gTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGaaa                                                                                                                                 >  1:1820284/1‑150 (MQ=255)
                     tGAAAATAAACGAAGCGTTTCAGGACTACGCTTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCTACACCCGAAAAACGTTACTTTTGAAAGGTCTTCTGTAACAATGATGACAGACATTACTGTATGTCAGATTTaaaa                                                                                                                                >  2:1423883/1‑150 (MQ=255)
                     tGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGaaaa                                                                                                                                <  1:1643279/150‑1 (MQ=255)
                      gAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAAt                                                                                                                               >  1:1323156/1‑150 (MQ=255)
                         aaTAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTa                                                                                                                            <  2:781728/150‑1 (MQ=255)
                             aaCGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTAATGGCTGTCAGATTGAAAATTTAGTgg                                                                                                                        >  2:1079673/1‑150 (MQ=255)
                               cGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGtt                                                                                                                      >  1:252907/1‑150 (MQ=255)
                                gAATAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTTCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTa                                                                                                                     <  1:1034561/150‑1 (MQ=255)
                                     gTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGaa                                                                                                                >  2:750875/1‑150 (MQ=255)
                                     gTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGaa                                                                                                                >  2:1660997/1‑150 (MQ=255)
                                      tttCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAAt                                                                                                               <  1:1030994/150‑1 (MQ=255)
                                       ttCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAAtt                                                                                                              >  1:1665512/1‑150 (MQ=255)
                                       ttCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAAtt                                                                                                              >  2:318392/1‑150 (MQ=255)
                                        tCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTc                                                                                                             >  1:773087/1‑150 (MQ=255)
                                         cAGGACTACGCGTGGGTAAAACCTGAATATCTAGAGCATTATTATTTGAATGTCGACACCCGAAAAACGTTACGTTTGAAAGGTCTTCTCTAACAATGATGACAGCCATTACTGGATGTCAGATTGAAAATTTAGTGGTTAACTAATTca                                                                                                            >  1:839316/1‑149 (MQ=255)
                                         cAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGAGGTTAACGAATTcc                                                                                                            <  1:402856/150‑1 (MQ=255)
                                          aGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTccc                                                                                                           <  1:1740023/150‑1 (MQ=255)
                                          aGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATATGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACACCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTccc                                                                                                           <  1:1740689/150‑1 (MQ=255)
                                           ggACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTcccc                                                                                                          >  1:1531618/1‑150 (MQ=255)
                                             aCTACGCGTTGGTAAAACCTCAAGATTTGGTGCATTATGATTTGAATTTCGCCACCCGAAAAACGTTACGTTTGAAATGTCTTATGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCAt                                                                                                        <  1:1132044/150‑1 (MQ=255)
                                              gtACGCGTGGGTAAAACCTGAAGATCTGGTGCATTAAGATTAGAATCTCGCCACCCGAAAAACGATACGTTTGAAAGGTCTTCTGAAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATc                                                                                                       >  1:338523/2‑150 (MQ=255)
                                              cTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATc                                                                                                       >  2:1591029/1‑150 (MQ=255)
                                               tACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATcc                                                                                                      <  1:1809840/150‑1 (MQ=255)
                                                aCGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCa                                                                                                     <  2:313963/150‑1 (MQ=255)
                                                aCGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCa                                                                                                     >  1:597344/1‑150 (MQ=255)
                                                 cgcgTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAg                                                                                                    <  1:1722631/150‑1 (MQ=255)
                                                  tagtgtgtAAAACATGAAGATCTGGTGAATTATGATTTGAATGTCGCAACCCGAAAAACGTTACTTTTTAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGa                                                                                                   <  1:324860/148‑1 (MQ=255)
                                                   cgTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGAt                                                                                                  >  2:1604922/1‑150 (MQ=255)
                                                   cgTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGAt                                                                                                  >  1:1177569/1‑150 (MQ=255)
                                                    gggggTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATa                                                                                                 <  2:1619782/148‑1 (MQ=255)
                                                    gTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATa                                                                                                 >  2:949880/1‑150 (MQ=255)
                                                         tAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACtt                                                                                            >  1:1368900/1‑150 (MQ=255)
                                                             aCCTGAAGATCTGGTGCATTATGATATGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACACTGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGcc                                                                                        <  2:196739/150‑1 (MQ=255)
                                                              ccTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCt                                                                                       >  1:511712/1‑150 (MQ=255)
                                                              ccTGAAGATCTGGTGCATGATGATTTGAATGTCGCCACCAGAAAAACGTTACGTTTGAAATGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCACCATCCAGATAAACTTTTCCt                                                                                       <  2:1168798/150‑1 (MQ=255)
                                                                tGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCttt                                                                                     >  2:673799/1‑150 (MQ=255)
                                                                 gAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCtttt                                                                                    >  1:1017304/1‑150 (MQ=255)
                                                                   aGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGTCTTTTTc                                                                                  <  1:1068449/150‑1 (MQ=255)
                                                                    gATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTcc                                                                                 <  1:1178937/150‑1 (MQ=255)
                                                                     aTCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCca                                                                                <  1:1708659/150‑1 (MQ=255)
                                                                      tCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCcac                                                                               <  2:1218772/150‑1 (MQ=255)
                                                                       cTGGTGCATTATGATTTGAATGGCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCcaca                                                                              <  2:1383416/150‑1 (MQ=255)
                                                                        tGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcat                                                                             >  1:739356/1‑150 (MQ=255)
                                                                         ggTGCATTATGATTTGCATGTCGCCACCCGAAAAACGGTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatg                                                                            <  2:1721421/150‑1 (MQ=255)
                                                                         ggTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatg                                                                            <  2:1265653/150‑1 (MQ=255)
                                                                          gTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTGCGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgc                                                                           <  1:188012/150‑1 (MQ=255)
                                                                            gCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAACCTTTGCCTTTTTCCAcatgcat                                                                         >  2:1123128/1‑150 (MQ=255)
                                                                             cATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgcatg                                                                        <  1:52/150‑1 (MQ=255)
                                                                              aTTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATTCATGa                                                                       >  2:45054/1‑150 (MQ=255)
                                                                                tATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACt                                                                     >  2:966617/1‑150 (MQ=255)
                                                                                   gATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTaac                                                                  >  2:733802/1‑149 (MQ=255)
                                                                                    aTTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAc                                                                 >  2:963468/1‑150 (MQ=255)
                                                                                     tttGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGCCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAcc                                                                >  2:1381368/1‑150 (MQ=255)
                                                                                        gAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATc                                                             <  2:1518983/150‑1 (MQ=255)
                                                                                             tCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGAt                                                        >  1:28459/1‑150 (MQ=255)
                                                                                                cctcccGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTATCGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAAt                                                     <  2:617903/147‑1 (MQ=255)
                                                                                                  aCCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCa                                                   <  1:991421/150‑1 (MQ=255)
                                                                                                   cccGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAg                                                  >  1:140595/1‑150 (MQ=255)
                                                                                                   cccGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAg                                                  >  1:1336747/1‑150 (MQ=255)
                                                                                                    ccGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAgg                                                 <  1:1510640/150‑1 (MQ=255)
                                                                                                      gAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGtt                                               >  2:1682735/1‑150 (MQ=255)
                                                                                                         aaaCGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTaa                                            <  1:242246/150‑1 (MQ=255)
                                                                                                            cgttacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATg                                         <  2:597344/150‑1 (MQ=255)
                                                                                                             gttacgttTGAAAGGTCTTCTGTATCAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTTCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAAAGt                                        >  2:549250/1‑150 (MQ=255)
                                                                                                                 cgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTcgcg                                    >  2:367849/1‑150 (MQ=255)
                                                                                                                  gttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGt                                   >  1:484415/1‑150 (MQ=255)
                                                                                                                   tttGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTc                                  >  2:445140/1‑150 (MQ=255)
                                                                                                                   tttGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTc                                  >  2:1448987/1‑150 (MQ=255)
                                                                                                                   tttGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTAAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCACATAAACTTTGCCTTTTTCCACATGCATGACTACTCCCTCCAGATAATCAGTTTTAAATGTCGCGTc                                  >  1:281886/1‑150 (MQ=255)
                                                                                                                    ttGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCa                                 >  2:1133896/1‑150 (MQ=255)
                                                                                                                            tcttctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATc                         <  2:1017304/150‑1 (MQ=255)
                                                                                                                                ctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTa                     >  2:796155/1‑150 (MQ=255)
                                                                                                                                   tAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCa                  <  1:1591029/150‑1 (MQ=255)
                                                                                                                                   tAACAATGATGACAGCCATTACTGGATGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACATTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCa                  <  1:318392/150‑1 (MQ=255)
                                                                                                                                    aaCAATGATGACAGCCATTACTGGCTTTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCaa                 <  2:488018/150‑1 (MQ=255)
                                                                                                                                    aaCAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAACATCCAGATAATCAGGTTTAATTGTCGCGTCAAGCTTATATTGAGaaa                 >  2:1143287/1‑150 (MQ=255)
                                                                                                                                     ttcaatctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAat                >  2:1640130/8‑150 (MQ=255)
                                                                                                                                     aCAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAat                >  1:1203357/1‑150 (MQ=255)
                                                                                                                                      cAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAata               >  1:731833/1‑150 (MQ=255)
                                                                                                                                      cAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTAAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATCCATGACTACACCATCCAGATAATCAGGTATAAATGTAGCGTCACTCTTATCATTAGCAata               >  2:1451344/1‑150 (MQ=255)
                                                                                                                                       aaTGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAatat              <  2:140595/150‑1 (MQ=255)
                                                                                                                                           ctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGACTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGt          >  1:1622468/2‑150 (MQ=255)
                                                                                                                                            tgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTc         <  2:723713/150‑1 (MQ=255)
                                                                                                                                            tgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTc         >  1:1192857/1‑150 (MQ=255)
                                                                                                                                             taCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCa        <  2:918442/149‑1 (MQ=255)
                                                                                                                                                aGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAAATCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGc     <  2:810381/150‑1 (MQ=255)
                                                                                                                                                  ccATTACTGGCAGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCaa   >  2:488071/1‑150 (MQ=255)
                                                                                                                                                   cATTACTGGCTGTCAGATTGAAAATTTAGTTGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAAt  <  2:1605613/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCGTTTCTGCCAACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAAT  >  NC_000913/2364846‑2365142

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