Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,212,477 G→A 5.4% intergenic (‑474/+226) iraM ← / ← ycgX RpoS stabilzer during Mg starvation, anti‑RssB factor/DUF1398 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,212,4770GA 5.4% 419.1 / 7.3 130intergenic (‑474/+226)iraM/ycgXRpoS stabilzer during Mg starvation, anti‑RssB factor/DUF1398 family protein
Reads supporting (aligned to +/- strand):  ref base G (52/70);  new base A (3/4);  total (55/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

Only 100 of 134 total aligned reads displayed.

AATAATGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGA  >  NC_000913/1212330‑1212626
                                                                                                                                                   |                                                                                                                                                     
aataatGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCa                                                                                                                                                     <  2:342412/150‑1 (MQ=255)
 ataatGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCAc                                                                                                                                                    <  2:318393/150‑1 (MQ=255)
  taatGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTTCGTTAAGTATTGAAACGTACGTGTGGCACg                                                                                                                                                   >  2:48247/1‑150 (MQ=255)
    atGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGc                                                                                                                                                 <  1:1590233/150‑1 (MQ=255)
    atGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGc                                                                                                                                                 >  2:1679134/1‑150 (MQ=255)
       aaTATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGCCGTGCCACACGTCCgt                                                                                                                                              <  1:1129662/150‑1 (MQ=255)
         tataAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCgttg                                                                                                                                            <  1:366014/150‑1 (MQ=255)
         aataAGCATGACATTAGAATAAGGTTTCAATTTTTGATTTATATAGTAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACTTGTGCGTTAAGTATTGAAACGGCCGTGCCACACGTCCgttt                                                                                                                                            <  1:739882/149‑2 (MQ=255)
          ataAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGCCGTGCCACACGTCCgtttc                                                                                                                                           <  1:1298165/150‑3 (MQ=255)
          ataAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGCCGTGCCACACGTCCgtttc                                                                                                                                           <  1:1239905/150‑3 (MQ=255)
           taAGCATGACATGAGAATAAGGATTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCgttgtt                                                                                                                                          <  2:409907/150‑1 (MQ=255)
             aGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGtttt                                                                                                                                        <  2:424731/150‑1 (MQ=255)
             aGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGtttt                                                                                                                                        <  2:148749/150‑1 (MQ=255)
              gCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGttttt                                                                                                                                       <  2:1366222/150‑1 (MQ=255)
               catgacatgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTta                                                                                                                                      >  1:426177/1‑150 (MQ=255)
                atgccatgaCAATAAGGTTTCAATTTTTGAGTTATATAGGAAAGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTtat                                                                                                                                     <  1:335272/150‑1 (MQ=255)
                atgacatgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTtat                                                                                                                                     <  1:1407472/150‑1 (MQ=255)
                  gacatgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATaa                                                                                                                                   <  2:1351908/150‑1 (MQ=255)
                   acatgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATaaa                                                                                                                                  >  2:1606241/1‑150 (MQ=255)
                      tgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAatat                                                                                                                               <  1:481254/150‑1 (MQ=255)
                      tgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAatat                                                                                                                               >  2:1104452/1‑150 (MQ=255)
                      tgaGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATAAAACCGGATTATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACTGACGTGTGGCACGGCCGTTGTTTTTTTAAatat                                                                                                                               <  1:1021756/150‑1 (MQ=255)
                         gaATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGTACGTGTGGCACGGCCGTTGTTTTTATAAATATGtt                                                                                                                            >  1:1668566/1‑150 (MQ=255)
                           aTAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTaa                                                                                                                          <  2:706532/150‑1 (MQ=255)
                               ggTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACAGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGt                                                                                                                      >  2:1639973/1‑150 (MQ=255)
                                gTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGtt                                                                                                                     <  2:441508/150‑1 (MQ=255)
                                  ttCAATTTTTTATTTATATATGAATTTTTTAACCTGTTCCTGGCTAAAATACATTTAACCGGATGATGACTAAACCATAATACATGTGCGTTATGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTTACCGTtat                                                                                                                   <  1:517560/150‑1 (MQ=255)
                                    cAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATaa                                                                                                                 >  2:368635/1‑150 (MQ=255)
                                    cAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATaa                                                                                                                 <  2:303573/150‑1 (MQ=255)
                                        ttttGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATa                                                                                                             >  1:650885/1‑150 (MQ=255)
                                           tGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCTGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACg                                                                                                          <  2:232452/150‑1 (MQ=255)
                                           tGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGTATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACg                                                                                                          <  1:987372/150‑1 (MQ=255)
                                           tGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACg                                                                                                          <  1:1189874/150‑1 (MQ=255)
                                               ttATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATc                                                                                                      <  1:704611/150‑1 (MQ=255)
                                                 atataGGAATGATTTAACCTGTTCCTGGCTAAAATAAATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATTTTATAACggttagt                                                                                                    <  1:1639973/150‑7 (MQ=255)
                                                    taGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAg                                                                                                 <  1:1402993/150‑1 (MQ=255)
                                                    taGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAg                                                                                                 <  1:1649937/150‑1 (MQ=255)
                                                    taGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAg                                                                                                 <  2:371274/150‑1 (MQ=255)
                                                     aGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGt                                                                                                >  2:1322632/1‑150 (MQ=255)
                                                      ggAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTc                                                                                               <  2:804645/150‑1 (MQ=255)
                                                      ggAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTc                                                                                               >  2:1053186/1‑150 (MQ=255)
                                                       gAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCa                                                                                              >  1:100799/1‑150 (MQ=255)
                                                        aaTGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCaa                                                                                             >  2:768438/1‑150 (MQ=255)
                                                           gATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTg                                                                                          >  1:645023/1‑150 (MQ=255)
                                                           gATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTg                                                                                          <  2:507430/150‑1 (MQ=255)
                                                           gATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTg                                                                                          >  2:1229663/1‑150 (MQ=255)
                                                                 aCCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATcaca                                                                                    <  1:1571306/150‑1 (MQ=255)
                                                                        ccTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCaaa                                                                             <  2:1722368/150‑1 (MQ=255)
                                                                          tGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAata                                                                           <  1:1436257/150‑1 (MQ=255)
                                                                             cTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATGTGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCa                                                                        <  2:581612/150‑1 (MQ=255)
                                                                             cTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCa                                                                        <  2:150516/150‑1 (MQ=255)
                                                                               aaaaTACATATAACCGTATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAg                                                                      <  1:599392/150‑1 (MQ=255)
                                                                                aaaTACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGt                                                                     <  1:1078970/150‑1 (MQ=255)
                                                                                 aaTACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGtt                                                                    <  1:1672319/150‑1 (MQ=255)
                                                                                  aTACATATAAGCGGATGATGACTAATCCATAATACATGTGCGTTTAGTATTGTATCGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATATAATAACGTATCAAAAGTCTAGTGATCACATTTCAAATATCAAGTTg                                                                   <  2:446399/150‑1 (MQ=255)
                                                                                   tACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTAAAAAATAACGTAAAAAAAGTCAAGTGATCACATTTAAAATATCAAGttta                                                                  >  1:1476626/1‑148 (MQ=255)
                                                                                         taACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTAt                                                            >  1:1640106/1‑150 (MQ=255)
                                                                                            ccGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTa                                                             <  1:502649/146‑1 (MQ=255)
                                                                                            ccGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTa                                                             >  2:502649/1‑146 (MQ=255)
                                                                                              ggATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTc                                                       <  2:8745/150‑1 (MQ=255)
                                                                                                atgatgACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTCGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTg                                                     <  2:91251/150‑1 (MQ=255)
                                                                                                  gatgaCTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGt                                                   <  1:618473/150‑1 (MQ=255)
                                                                                                  gatgaCTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGt                                                   <  1:908151/150‑1 (MQ=255)
                                                                                                   atgaCTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTg                                                  >  1:298584/1‑150 (MQ=255)
                                                                                                   atgaCTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATAAAGTTGATAGTATTAGTATGGTg                                                  >  1:448218/1‑150 (MQ=255)
                                                                                                    tgaCTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGa                                                 >  2:632554/1‑150 (MQ=255)
                                                                                                     gaCTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGat                                                >  2:399476/1‑150 (MQ=255)
                                                                                                         aaaCCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGattatt                                            >  2:1574383/1‑150 (MQ=255)
                                                                                                         aaaCCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACAATTCAAATATCAAGTTGATATGATTAGTATGGTGattatt                                            >  2:186096/1‑150 (MQ=255)
                                                                                                              aaaaTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATgg                                       >  2:667366/1‑150 (MQ=255)
                                                                                                               aaaTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATggg                                      <  1:1053186/150‑1 (MQ=255)
                                                                                                               aaaTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATggg                                      <  1:1104452/150‑1 (MQ=255)
                                                                                                               aaaTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATggg                                      <  2:818459/150‑1 (MQ=255)
                                                                                                               aaaTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATggg                                      <  1:646548/150‑1 (MQ=255)
                                                                                                               aaaTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATggg                                      <  1:552036/150‑1 (MQ=255)
                                                                                                                 aTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTg                                    <  2:98464/150‑1 (MQ=255)
                                                                                                                  tACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGa                                   <  2:431083/150‑1 (MQ=255)
                                                                                                                       gtgCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGACTGGTGATTATATATGTGAGACAATa                              >  2:1643713/1‑150 (MQ=255)
                                                                                                                        tgCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATaa                             >  1:300775/1‑150 (MQ=255)
                                                                                                                        tgCGTTAAGTATTGAAACGGACGTGTGGCAAGGCCGTTGTTATTATATATATGTTAACCGTTATAAAATAAAGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTTATAGTATTAGTCTGGTGATTATATATGGGTTACAATaa                             >  1:662966/1‑150 (MQ=255)
                                                                                                                         gCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATTGGTGACAATaaa                            >  2:1307103/1‑150 (MQ=255)
                                                                                                                         gCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATaaa                            >  1:366207/1‑150 (MQ=255)
                                                                                                                          cGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATaaaa                           >  1:919184/1‑150 (MQ=255)
                                                                                                                           gTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTAAAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATAAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATaaaaa                          >  2:235693/1‑150 (MQ=255)
                                                                                                                              aaGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGAc                       <  2:1648943/150‑1 (MQ=255)
                                                                                                                              aaGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGAc                       <  2:650885/150‑1 (MQ=255)
                                                                                                                               aGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACa                      >  1:1634372/1‑150 (MQ=255)
                                                                                                                                taaaaacaaCGGCCGTGCCACACGTCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAg                     >  2:552636/8‑150 (MQ=255)
                                                                                                                                 aaaaacaaCGGCCGTGCCACACGTCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGt                    >  2:990594/7‑150 (MQ=255)
                                                                                                                                  aTTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTa                   >  2:1079016/1‑150 (MQ=255)
                                                                                                                                  aTTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTa                   >  1:415679/1‑150 (MQ=255)
                                                                                                                                   ttGAAACGTACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAAGCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTAt                  <  2:855136/150‑1 (MQ=255)
                                                                                                                                    tGAAACGGACGTGTGGCACGGCCGTTGGTTTTAAAAATATGTTTACCGGTTTAAAAGTACGTTTCAAAAGTCAAGGGTAGACATTTGAAATTTGAAGTTTAAAGTTTTATTGTTGTGCTTTTTTTAGGGTGTCAATAAAAAGACAGTAtt                 >  2:1564584/1‑150 (MQ=255)
                                                                                                                                     gAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTa                >  2:269055/1‑150 (MQ=255)
                                                                                                                                     gAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTa                >  2:1204515/1‑150 (MQ=255)
                                                                                                                                      aaaCGGACGAGTGGCACGGCCGTTGATTTTATAAATATGTAAACCGTTATAAAATAACGAATCAAAAGTCAAGTGATCACATATCAAATATAAAGATGAATCTATTAGTCTCTTGATTATATATGGGTGACAATAAAAAGAAAGTATTaa               >  1:952/1‑150 (MQ=255)
                                                                                                                                        aCGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAatc             >  2:110173/1‑150 (MQ=255)
                                                                                                                                           gACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAatcatc          >  1:29303/1‑150 (MQ=255)
                                                                                                                                                  ggCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATaga  >  2:1710200/1‑149 (MQ=255)
                                                                                                                                                  ggCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATagag   >  2:1203540/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AATAATGAATATAAGCATGACATGAGAATAAGGTTTCAATTTTTGAGTTATATAGGAATGATTTAACCTGTTCCTGGCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGA  >  NC_000913/1212330‑1212626

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: