Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 480,236 C→G 15.1% 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.1% 481.9 / 69.7 178E25Q (GAG→CAG) hhamodulator of gene expression, with H‑NS
Reads supporting (aligned to +/- strand):  ref base C (83/66);  new base G (10/17);  total (94/84)
Fisher's exact test for biased strand distribution p-value = 9.40e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

Only 100 of 179 total aligned reads displayed.

ATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTC  >  NC_000913/480138‑480334
                                                                                                  |                                                                                                  
aTTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTcc                                                                                                   <  1:1371381/100‑1 (MQ=255)
 ttCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCa                                                                                                  >  1:2155982/1‑100 (MQ=255)
    aTGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATTACACGCTCCAGCg                                                                                               <  1:2670697/100‑1 (MQ=255)
    aTGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGcn                                                                                               <  2:180087/100‑2 (MQ=255)
    aTGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCg                                                                                               <  2:2030535/100‑1 (MQ=255)
      ggTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGAGAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCgtg                                                                                             >  2:2147930/1‑100 (MQ=255)
       gTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATGCGTATTGATTTTTCTCGATAACACGCTCCAGCg                                                                                               >  2:3575987/1‑97 (MQ=255)
       gTCAATTCAGCGAGGCGGTGATCTGCGGCTTAGTAAAATACCGCCAGTTCATGTACTGCTAAGTCGTATTTATTTTTCTCGATAACACGCTCCAGCgtgt                                                                                            >  2:3551472/1‑100 (MQ=255)
        tCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTc                                                                                           <  1:1798575/100‑1 (MQ=255)
         cAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCa                                                                                          <  2:117478/100‑1 (MQ=255)
         cAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCa                                                                                          >  2:2744699/1‑100 (MQ=255)
         cAATTCGGCGAGGCGGTGATATGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATGCGTATTTATGTTACTCGATAACACGCTACAg                                                                                                 >  2:891340/1‑93 (MQ=255)
          aaTTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTGGAGCgtgttat                                                                                         <  1:1516832/100‑4 (MQ=255)
           aTTCGGCGAGGTGGTGATCTGCTGCGAAGTAAAATACCGCCAGTTCATTATCTGAGAATTCGTATTTATTTTTCGCGATAACACGCTCCAGCGTGTCAAt                                                                                        <  2:3381881/100‑1 (MQ=255)
           aTTCGGCGAGGCGGTGATCTGCGGCTTACTAAAATATCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTGGAGCgtgttatc                                                                                        <  2:3353804/100‑5 (MQ=255)
           aTTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTGGAGCgtgttatc                                                                                        <  1:98361/100‑5 (MQ=255)
           aTTCGGCGAGGCGCTGATCTGGGGCTTAGTAAAATACCGACAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTGGAGCgtgttatc                                                                                        <  2:368638/100‑5 (MQ=255)
             tCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAATACGCTCCAGCGTGTCAATTg                                                                                      <  1:3713552/100‑1 (MQ=255)
             tCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTGGAGCgtgttatcga                                                                                      <  2:1208746/100‑7 (MQ=255)
             tCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTGGAGCgtgttatcga                                                                                      <  1:1643887/100‑7 (MQ=255)
             tCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTGGAGCgtgttatcga                                                                                      <  2:1662350/100‑7 (MQ=255)
              cGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTGGAGCgtgttatcgag                                                                                     <  2:2310539/100‑8 (MQ=255)
               ggCGTGGCTGTGAGCTGCGGCTTAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATACCACGCTGGAGCgtgttatcgag                                                                                     <  2:2283484/99‑8 (MQ=255)
               ggCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTc                                                                                    >  1:1435191/1‑100 (MQ=255)
                gCGAGGCGGTGATCTGCGGCTGAGTAAAATTCCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTGTAGCgtgttatcg                                                                                       <  2:1238299/96‑6 (MQ=255)
                 cGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTg                                                                                  >  1:3948022/1‑100 (MQ=255)
                 cGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTg                                                                                  <  2:2264799/100‑1 (MQ=255)
                 cGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTg                                                                                  >  1:1771402/1‑100 (MQ=255)
                   attctgTGATGTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATGATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCt                                                                                   <  2:525471/92‑1 (MQ=255)
                       ggTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAacg                                                                            >  1:241509/1‑100 (MQ=255)
                        gTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAacga                                                                           >  1:973819/1‑100 (MQ=255)
                        gTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAacg                                                                            >  2:3241206/1‑99 (MQ=255)
                          gATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAacgacg                                                                         >  1:1719228/1‑100 (MQ=255)
                          gATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAacgacg                                                                         >  1:3502955/1‑100 (MQ=255)
                          gATATGCGGCTGGGTAAAATACCGCCAGTTCATTATGTGAGAATTCGTATTTATTTTTCTCGATAACACGCTCCGGCGTGTCAATTGTCTGGCAacgacg                                                                         <  1:1935752/100‑1 (MQ=255)
                           aTCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGt                                                                        >  1:189998/1‑100 (MQ=255)
                             cTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATCATTCGTATTTATTTTTCTCGATAACACTCTCCAGCGTGTCCATTTTCTGGCAACGACGt                                                                        >  2:2631694/1‑98 (MQ=255)
                             cTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTaa                                                                      >  1:240477/1‑100 (MQ=255)
                              tGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTaaa                                                                     >  1:959201/1‑100 (MQ=255)
                                cGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACg                                                                   >  1:3171522/1‑100 (MQ=255)
                                cGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACg                                                                   >  2:3928724/1‑100 (MQ=255)
                                    tGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCAtt                                                               <  1:1017993/100‑1 (MQ=255)
                                    tGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGCCAATTGTCTGGCAACGACGTAAACGCAtt                                                               <  1:1533579/100‑1 (MQ=255)
                                      aGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTaa                                                             >  1:1258924/1‑100 (MQ=255)
                                       gTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTATCTGGCAACGACGTAAACGCATTa                                                              >  2:1878989/1‑98 (MQ=255)
                                            aTACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATc                                                       <  1:2141366/100‑1 (MQ=255)
                                             tACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCg                                                      <  2:2873799/100‑1 (MQ=255)
                                              aCCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCgg                                                     >  2:521542/1‑100 (MQ=255)
                                               ccGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGt                                                    >  1:752841/1‑100 (MQ=255)
                                               ccGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGt                                                    <  2:722387/100‑1 (MQ=255)
                                                cGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGtt                                                   >  2:93960/1‑100 (MQ=255)
                                                 gCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCgg                                                     >  2:3779738/1‑97 (MQ=255)
                                                 gCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGttt                                                  >  1:1577147/1‑100 (MQ=255)
                                                 gCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGGCAATTGTCTGGCAACGACGTAAACCCATTAAATAATCGGttt                                                  >  2:3291316/1‑100 (MQ=255)
                                                     gTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGt                                              >  2:1947741/1‑100 (MQ=255)
                                                      ttCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGtt                                             <  2:344577/100‑1 (MQ=255)
                                                      ttCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGtt                                             <  2:189998/100‑1 (MQ=255)
                                                       tCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTa                                            <  2:1579738/100‑1 (MQ=255)
                                                       tCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTa                                            >  2:2968458/1‑100 (MQ=255)
                                                        cATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGtt                                             <  2:644452/98‑1 (MQ=255)
                                                         aTTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTaaa                                          >  1:711458/1‑100 (MQ=255)
                                                           tATCTGATAATTCGGATTTATTTTTCTCGATAACACGCTCCCGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAgg                                        <  2:2427936/100‑1 (MQ=255)
                                                            aTCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGt                                       >  1:1725953/1‑100 (MQ=255)
                                                               tGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGtttt                                    >  1:951552/1‑100 (MQ=255)
                                                                   aaTTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCgg                                <  1:3311609/100‑1 (MQ=255)
                                                                     ttCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGAc                              <  1:3447603/100‑1 (MQ=255)
                                                                     ttCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGAc                              <  1:1046630/100‑1 (MQ=255)
                                                                     ttCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGAc                              <  1:3294692/100‑1 (MQ=255)
                                                                      tCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACa                             <  1:455290/100‑1 (MQ=255)
                                                                          atttatttTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACt                         <  2:1792245/100‑1 (MQ=255)
                                                                          atttatttTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACt                         >  1:889814/1‑100 (MQ=255)
                                                                          atttatttTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACt                         <  2:3129677/100‑1 (MQ=255)
                                                                             tatttTTCTCGATAACACGCTCCAGCGTGTCAATTGTATGGCAACGACATAAACGCAATAAATAATCGGTTTTCGCTAAAAGCTTTTCGGAAATACtt                        >  2:733191/1‑98 (MQ=255)
                                                                               tttttCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTAc                    >  2:12013/1‑100 (MQ=255)
                                                                                ttttCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTAcc                   <  1:1000335/100‑1 (MQ=255)
                                                                                 tttCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCt                  >  1:467540/1‑100 (MQ=255)
                                                                                  agacaattgacacGCTGGAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTa                 >  2:3140633/10‑100 (MQ=255)
                                                                                   tctcGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTAt                <  1:2295736/100‑1 (MQ=255)
                                                                                   gacaattgacacGCTGGAGCGTGTTAATTGTCTGGCAACGACGTAAACGCATTCAATAATCGGTTTTCGTTAAAGATTTTTCGGACATACTTCTACCTAt                >  2:649940/9‑100 (MQ=37)
                                                                                    ctcGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATg               >  2:1079464/1‑100 (MQ=255)
                                                                                    ctcGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATg               >  1:1908004/1‑100 (MQ=255)
                                                                                    ctcGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATg               >  2:552202/1‑100 (MQ=255)
                                                                                    ctcGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATg               >  1:3124439/1‑100 (MQ=255)
                                                                                     tcGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATgg              >  2:1824775/1‑100 (MQ=255)
                                                                                     tcGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATgg              >  1:1197284/1‑100 (MQ=255)
                                                                                      cGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTATAGGTTTTTCGGACATACTTCTACCTATGGt             <  2:2588498/100‑1 (MQ=255)
                                                                                      cGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGt             >  1:2430025/1‑100 (MQ=255)
                                                                                      cGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGt             >  1:1102943/1‑100 (MQ=255)
                                                                                      cGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGt             >  1:3034469/1‑100 (MQ=255)
                                                                                       attgacacGCTGGAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGtt            >  1:3075827/5‑100 (MQ=255)
                                                                                       attgacacGCTGGAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGtt            >  2:2193557/5‑100 (MQ=255)
                                                                                       attgacacGCTGGAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGtt            >  1:2872942/5‑100 (MQ=255)
                                                                                        aTAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTg           <  2:425138/100‑1 (MQ=255)
                                                                                        aTAACACACTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCTGACATACTTCTACCTATGGTTg           <  2:1475346/100‑1 (MQ=255)
                                                                                         tgacacGCTGGAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGt          >  1:3640655/3‑100 (MQ=255)
                                                                                         tAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGt          <  1:2986852/100‑1 (MQ=255)
                                                                                           acacGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGtaa        >  2:3811145/1‑100 (MQ=255)
                                                                                           acacGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGtaa        >  1:625095/1‑100 (MQ=255)
                                                                                           acacGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGtaa        >  2:690360/1‑100 (MQ=255)
                                                                                                 tCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTc  <  2:486097/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
ATTCATGGTCAATTCGGCGAGGCGGTGATCTGCGGCTGAGTAAAATACCGCCAGTTCATTATCTGATAATTCGTATTTATTTTTCTCGATAACACGCTCCAGCGTGTCAATTGTCTGGCAACGACGTAAACGCATTAAATAATCGGTTTTCGTTAAAGGTTTTTCGGACATACTTCTACCTATGGTTGTAATAATTC  >  NC_000913/480138‑480334

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