Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 480,236 C→G 15.5% E25Q (GAG→CAG)  hha ← modulator of gene expression, with H‑NS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913480,2360CG15.5% 282.8 / 41.2 103E25Q (GAG→CAG) hhamodulator of gene expression, with H‑NS
Reads supporting (aligned to +/- strand):  ref base C (49/38);  new base G (7/9);  total (56/47)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 105 total aligned reads displayed.

TATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATT  >  NC_000913/480137‑480333
                                                                                                   |                                                                                                 
tATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTc                                                                                                   <  2:329307/100‑1 (MQ=255)
    cATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGc                                                                                               >  1:394828/1‑100 (MQ=255)
     aTGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCg                                                                                              <  2:1570117/100‑1 (MQ=255)
     aTGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCg                                                                                              <  2:1364739/100‑1 (MQ=255)
       ggTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCgtg                                                                                            <  1:2335404/100‑1 (MQ=255)
       ggTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCgtg                                                                                            >  1:2287188/1‑100 (MQ=255)
       ggTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCgtg                                                                                            >  2:1462643/1‑100 (MQ=255)
       ggTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCgtg                                                                                            <  2:455724/100‑1 (MQ=255)
          cAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCa                                                                                         >  2:646903/1‑100 (MQ=255)
           aaTTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACATGCTCCAGCGTGTCaa                                                                                        >  1:2517493/1‑100 (MQ=255)
           aaTTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCaa                                                                                        >  1:1883073/1‑100 (MQ=255)
           aaTTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCaa                                                                                        <  1:2536364/100‑1 (MQ=255)
           aaTTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCaa                                                                                        >  2:1614306/1‑100 (MQ=255)
           aaTTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCaa                                                                                        >  1:1615510/1‑100 (MQ=255)
            aTTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTGGAGCgtgttatc                                                                                       <  1:1901339/100‑5 (MQ=255)
            aTTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTGGAGCgtgttatc                                                                                       <  1:2078411/100‑5 (MQ=255)
             ttCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGTCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTGGAGCgtgttatcg                                                                                      <  2:213054/100‑6 (MQ=255)
             ttCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTGGAGCgtgttatcg                                                                                      <  1:1901210/100‑6 (MQ=255)
             ttCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTGGAGCgtgttatcg                                                                                      <  2:1387697/100‑6 (MQ=255)
              tCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTGGAGCgtgttatcga                                                                                     <  1:1592247/100‑7 (MQ=255)
              tCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTGGAGCgtgttatcga                                                                                     <  2:683350/100‑7 (MQ=255)
                ggCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTc                                                                                   >  2:2509493/1‑100 (MQ=255)
                  cGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTg                                                                                 <  1:2222401/100‑1 (MQ=255)
                  cGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTg                                                                                 <  1:2222341/100‑1 (MQ=255)
                   gAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTgg                                                                                <  1:1385411/100‑1 (MQ=255)
                    aGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGc                                                                               <  1:2199367/100‑1 (MQ=255)
                     ggCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCa                                                                              >  1:852000/1‑100 (MQ=255)
                     ggCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCa                                                                              <  2:1188584/100‑1 (MQ=255)
                       cGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAac                                                                            >  2:1165208/1‑100 (MQ=255)
                       cGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAac                                                                            >  1:698642/1‑100 (MQ=255)
                        ggTGATCTGCGGCTGAGTACACTACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATCACCCGCTCCAGCGTGTCAATTGTCTGGCAacg                                                                           <  1:2041619/100‑1 (MQ=255)
                           gATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAacgacg                                                                        <  1:1785593/100‑1 (MQ=255)
                           gATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAacgacg                                                                        >  1:2361425/1‑100 (MQ=255)
                           gATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATGTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAacgacg                                                                        <  2:1926606/100‑1 (MQ=255)
                             tCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTa                                                                      >  1:1443007/1‑100 (MQ=255)
                               tGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTaaa                                                                    >  1:1506712/1‑100 (MQ=255)
                                  ggCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTCAACGc                                                                 >  2:215107/1‑100 (MQ=255)
                                   gCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCa                                                                >  2:797209/1‑100 (MQ=255)
                                      gAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTa                                                             >  1:2479655/1‑100 (MQ=255)
                                       aGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTaa                                                            <  1:725139/100‑1 (MQ=255)
                                        gTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACCCGCTCCAGCGTGTCAATTGTCTGTCAACGACGTCACCGCATTaa                                                            >  2:1704121/1‑99 (MQ=255)
                                        gTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTaaa                                                           >  1:2421958/1‑100 (MQ=255)
                                        gTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTaaa                                                           >  2:1902951/1‑100 (MQ=255)
                                        gTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTaaa                                                           <  1:1158308/100‑1 (MQ=255)
                                            aaTACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAaataat                                                       >  2:2513170/1‑100 (MQ=255)
                                             aTACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATc                                                      <  2:1828175/100‑1 (MQ=255)
                                              tACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGGCGTAAACGCATTAAATAATCg                                                     <  2:619199/100‑1 (MQ=255)
                                              tACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCg                                                     <  2:2503760/100‑1 (MQ=255)
                                                ccGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGt                                                   >  1:772380/1‑100 (MQ=255)
                                                ccGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGt                                                   >  2:701432/1‑100 (MQ=255)
                                                ccGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGt                                                   >  1:970715/1‑100 (MQ=255)
                                                ccGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGt                                                   >  1:970601/1‑100 (MQ=255)
                                                 cGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGtt                                                  >  1:1744316/1‑100 (MQ=255)
                                                  gCCAGTTCATTAGCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGttt                                                 >  1:1685396/1‑100 (MQ=255)
                                                     aGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCg                                              <  1:1474200/100‑1 (MQ=255)
                                                       ttCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGtt                                            <  1:345719/100‑1 (MQ=255)
                                                        tCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCTTGTCAATTGTCTGGCAACGACTTAAACGCATTAAATAATCGGTTTTCGTTa                                           >  2:1083244/1‑100 (MQ=255)
                                                        tCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGtt                                            <  2:1974557/99‑1 (MQ=255)
                                                          aTTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTACACGCATTAAATAATCGGTTTTCGTTaaa                                         >  2:1134573/1‑100 (MQ=255)
                                                          aTTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTaaa                                         >  1:2023582/1‑100 (MQ=255)
                                                           ttATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAg                                        <  1:1614317/100‑1 (MQ=255)
                                                           ttATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAg                                        <  1:147152/100‑1 (MQ=255)
                                                             aTCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGt                                      <  2:537253/100‑1 (MQ=255)
                                                                tGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGtttt                                   >  2:2216172/1‑100 (MQ=255)
                                                                tGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGtttt                                   >  1:1960226/1‑100 (MQ=255)
                                                                tGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGtttt                                   >  1:1979297/1‑100 (MQ=255)
                                                                 gATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGttttt                                  >  1:1071060/1‑100 (MQ=255)
                                                                  aTAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTc                                 <  1:797216/100‑1 (MQ=255)
                                                                   tAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCg                                >  1:220014/1‑100 (MQ=255)
                                                                   tAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCg                                <  2:1979329/100‑1 (MQ=255)
                                                                    aaTTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCgg                               >  1:809666/1‑100 (MQ=255)
                                                                     aTTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGa                              <  1:647664/100‑1 (MQ=255)
                                                                      ttCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGAc                             <  2:1766952/100‑1 (MQ=255)
                                                                         gTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATa                          <  1:1385097/100‑1 (MQ=255)
                                                                           atttatttTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGCAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACa                            <  2:290365/96‑1 (MQ=255)
                                                                             ttatttTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTc                      <  2:772373/100‑1 (MQ=255)
                                                                              tatttTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCt                     >  1:2273511/1‑100 (MQ=255)
                                                                                tttttCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTCAAGGTTTTTCGGACATACTTCTAc                   >  2:1094102/1‑100 (MQ=255)
                                                                                tttttCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTAc                   <  1:1793741/100‑1 (MQ=255)
                                                                                tttttCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTAc                   >  2:574507/1‑100 (MQ=255)
                                                                                  tttCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCt                 >  2:2291697/1‑100 (MQ=255)
                                                                                  cagacaattgacacGCTGGAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCt                 >  1:1951988/11‑100 (MQ=38)
                                                                                     ctcGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATg              >  1:1489366/1‑100 (MQ=255)
                                                                                      tcGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATgg             >  2:727055/1‑100 (MQ=255)
                                                                                      caattgacacGCTGGAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGTTTTTTCGGACATATTTCTACCTATg              >  2:2551525/7‑99 (MQ=255)
                                                                                       cGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGt            >  1:1284203/1‑100 (MQ=255)
                                                                                       aattgacacGCTGGAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATgg             >  1:904218/6‑99 (MQ=255)
                                                                                        attgacacGCTGGAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGtt           >  1:60702/5‑100 (MQ=255)
                                                                                        attgacacGCTGGAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGtt           >  1:1481200/5‑100 (MQ=255)
                                                                                        attgacacGCTGGAGCGTGTCAATTGTCTGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTg          >  2:1763648/5‑100 (MQ=255)
                                                                                          tgacacGCTGGAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGt         >  2:494479/3‑100 (MQ=255)
                                                                                          tAACACGCTCCAGCGTGTCAATTGTCTGTCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGt         <  1:1471315/100‑1 (MQ=255)
                                                                                          tAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGt         <  2:751080/100‑1 (MQ=255)
                                                                                          tAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGt         <  1:1531751/100‑1 (MQ=255)
                                                                                            acacGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGtaa       >  1:592509/1‑100 (MQ=255)
                                                                                            acacGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGtaa       >  1:750639/1‑100 (MQ=255)
                                                                                             cacGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGtaat      <  2:733202/100‑1 (MQ=255)
                                                                                             cacGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGtaat      >  2:868073/1‑100 (MQ=255)
                                                                                                 cTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCTGTTTTTTTTCTACCTATGGTGGTAATAAtt  >  2:1180013/1‑100 (MQ=255)
                                                                                                 cTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAAtt  <  1:225075/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
TATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATT  >  NC_000913/480137‑480333

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