Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 454,853 T→C 20.9% intergenic (+64/‑280) bolA → / → tig stationary‑phase morphogene, transcriptional repressor for mreB; also regulator for dacA, dacC, and ampC/peptidyl‑prolyl cis/trans isomerase (trigger factor)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913454,8530TC20.9% 195.2 / 50.9 102intergenic (+64/‑280)bolA/tigstationary‑phase morphogene, transcriptional repressor for mreB; also regulator for dacA, dacC, and ampC/peptidyl‑prolyl cis/trans isomerase (trigger factor)
Reads supporting (aligned to +/- strand):  ref base T (42/38);  new base C (14/7);  total (57/45)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.84e-01

Only 100 of 102 total aligned reads displayed.

CTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAAT  >  NC_000913/454755‑454950
                                                                                                  |                                                                                                 
ctcctcCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTg                                                                                                  >  2:114295/1‑100 (MQ=255)
  cctccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCg                                                                                                >  1:1195182/1‑100 (MQ=255)
     cccTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGc                                                                                             >  2:513071/1‑100 (MQ=255)
       cTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCg                                                                                           >  2:2109241/1‑100 (MQ=255)
        tGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgt                                                                                          <  2:1988770/100‑1 (MQ=255)
         gTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgtt                                                                                         >  2:2198520/1‑100 (MQ=255)
          tCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttt                                                                                        <  1:2133139/100‑1 (MQ=255)
           cGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgtttt                                                                                       >  1:308356/1‑100 (MQ=255)
           cGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgtttt                                                                                       >  1:504032/1‑100 (MQ=255)
            gTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttg                                                                                      <  2:526108/100‑1 (MQ=255)
            gTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttg                                                                                      <  2:812432/100‑1 (MQ=255)
              ggAGCAGGAAGGATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGATGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                    <  1:1771851/100‑4 (MQ=255)
              ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                    <  2:2014110/100‑4 (MQ=255)
              ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                    <  2:511638/100‑4 (MQ=255)
              ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                    <  1:2641385/100‑4 (MQ=255)
               gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                   <  1:1564291/100‑5 (MQ=39)
                aGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGCTGCTAAAGAATTTATGAAACACGGCCAGCGGGCCgttttgtttt                                                                                  >  1:1002047/1‑100 (MQ=255)
                 gCAGGAAGCATCGTGTAAAAGCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcata                                                                                  <  1:565587/99‑6 (MQ=38)
                 gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGttttgttttg                                                                                 >  2:1159083/1‑100 (MQ=255)
                 gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGttttgttttg                                                                                 <  2:24193/100‑1 (MQ=255)
                  cAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcataaa                                                                                <  2:1193011/100‑8 (MQ=38)
                    ggAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTtttgttttgttt                                                                              >  1:305514/1‑98 (MQ=255)
                    ggAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCt                                                                              <  1:1294395/100‑1 (MQ=255)
                     gAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTg                                                                             <  1:2134964/100‑1 (MQ=255)
                     gAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTg                                                                             <  1:1028699/100‑1 (MQ=255)
                         cATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTACAGTATGTTGCTAAAGATTCTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGAtt                                                                         <  1:1588010/100‑1 (MQ=255)
                          aTCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGAttt                                                                        >  2:791948/1‑100 (MQ=255)
                            cgcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTg                                                                      <  1:2399098/100‑1 (MQ=255)
                             gcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgc                                                                     <  1:749707/100‑1 (MQ=255)
                             gcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgc                                                                     <  1:1937829/100‑1 (MQ=255)
                               gTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTTCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgc                                                                     <  2:1576361/98‑1 (MQ=255)
                               gTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgcgc                                                                   >  2:1243941/1‑100 (MQ=255)
                                tAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCt                                                                  >  2:1726615/1‑100 (MQ=255)
                                 aaaaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCtt                                                                 >  1:604670/1‑100 (MQ=255)
                                    aaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCttttt                                                              >  2:1929164/1‑100 (MQ=255)
                                     aCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTg                                                             <  2:2671893/100‑1 (MQ=255)
                                       gCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGcc                                                           <  2:2783937/100‑1 (MQ=255)
                                        cATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGccc                                                          <  1:11112/100‑1 (MQ=255)
                                         aTTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGccc                                                          <  1:1822306/99‑1 (MQ=255)
                                         aTTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCa                                                         <  2:2708330/100‑1 (MQ=255)
                                          tttGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAg                                                        <  2:1943469/100‑1 (MQ=255)
                                           ttGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGc                                                       <  2:511475/100‑1 (MQ=255)
                                              cAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCAtt                                                    >  2:1588622/1‑100 (MQ=255)
                                              cAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCAtt                                                    <  1:1929156/100‑1 (MQ=255)
                                                  cgatctcgcTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGa                                                <  2:1840183/94‑1 (MQ=255)
                                                   gTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGAc                                               <  2:2872798/100‑1 (MQ=255)
                                                      ggCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaa                                            >  2:1126907/1‑100 (MQ=255)
                                                       gcgcTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGGTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaa                                            >  1:2164510/1‑99 (MQ=255)
                                                        cgcTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCagaagatcg                                          >  1:1840175/1‑97 (MQ=255)
                                                        cgcTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaaaa                                          >  1:2784226/1‑100 (MQ=255)
                                                         gcTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAAt                                         <  2:440258/100‑1 (MQ=255)
                                                          cTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAAtt                                        >  1:2525746/1‑100 (MQ=255)
                                                           ttttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCGGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTg                                       >  1:2248496/1‑100 (MQ=255)
                                                           ttttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTg                                       <  1:1021350/100‑1 (MQ=255)
                                                            tttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGc                                      >  2:1012162/1‑100 (MQ=255)
                                                             ttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGcc                                     <  1:972907/100‑1 (MQ=255)
                                                             ttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGcc                                     <  2:1104616/100‑1 (MQ=255)
                                                             ttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGcc                                     >  1:1029860/1‑100 (MQ=255)
                                                              tCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGTCCAGCATTCAGAAGAAAATTGcac                                    >  1:573229/1‑98 (MQ=255)
                                                               ccAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCg                                   <  2:2070211/100‑1 (MQ=255)
                                                               ccAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCg                                   >  2:1716967/1‑100 (MQ=255)
                                                                cAGGATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCGCAGCATTCAGACGAAAATTGCCCgg                                  >  1:1885649/1‑100 (MQ=255)
                                                                 aGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCggg                                 <  2:2535644/100‑1 (MQ=255)
                                                                 aGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCggg                                 <  1:1177370/100‑1 (MQ=255)
                                                                 aGTAAGTGGCTAAAGAGTTTATGAAGAACGGCCTGCGGGCCGTTGCGTTTTGTCGGGATTGTGCGCTTTCGGCCCAGCATTCAGACGAAAATTGCCCggg                                 <  1:1656653/100‑1 (MQ=255)
                                                                      gTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtg                            >  2:2011697/1‑100 (MQ=255)
                                                                        tGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtgtg                          <  1:255136/100‑1 (MQ=255)
                                                                          cTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaa                        <  1:2144753/100‑1 (MQ=255)
                                                                           tAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaa                       <  2:1672965/100‑1 (MQ=255)
                                                                           tAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaa                       >  2:2738483/1‑100 (MQ=255)
                                                                            aaaGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaaa                      <  2:2156030/100‑1 (MQ=255)
                                                                            aaaGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaaa                      <  1:2384988/100‑1 (MQ=255)
                                                                            aaaGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaaa                      >  2:602667/1‑100 (MQ=255)
                                                                                agacaaaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATa                  >  2:1132700/10‑100 (MQ=38)
                                                                                aTTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATa                  <  1:2011690/100‑1 (MQ=255)
                                                                                  tttATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcg                >  2:2063340/1‑100 (MQ=255)
                                                                                   caaaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGCATTGTGAAAAAATAcgc               >  1:2061687/7‑100 (MQ=38)
                                                                                   caaaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcgc               >  1:2666590/7‑100 (MQ=38)
                                                                                   caaaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCCGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcgc               >  2:1657931/7‑100 (MQ=38)
                                                                                     aaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGa             >  2:812872/5‑100 (MQ=39)
                                                                                     aaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGa             >  2:885030/5‑100 (MQ=39)
                                                                                     aaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGa             >  2:1162389/5‑100 (MQ=39)
                                                                                     aaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGa             >  2:2497660/5‑100 (MQ=39)
                                                                                      aacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc            >  2:911173/4‑100 (MQ=255)
                                                                                       gAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGGATCGTGAAAAAACACGCGACa           >  1:351662/1‑100 (MQ=255)
                                                                                       gAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACa           >  2:2728994/1‑100 (MQ=255)
                                                                                       gAAAAACGGCCTGCGGGCCGTGTTGTTTTGTCTGGGTTTTGCGCTTTTTGCCCAGCATTCAGACGAAAGTTGCCCGGGACTTGTGACAAAATACGCGACa           >  2:1276506/1‑100 (MQ=255)
                                                                                       acaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACa           >  2:891512/3‑100 (MQ=255)
                                                                                        aaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCGCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAAGCGACAg          >  1:628448/1‑100 (MQ=255)
                                                                                         naaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgc         >  1:2294109/2‑100 (MQ=255)
                                                                                         aaaaCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgc         >  2:2024188/1‑100 (MQ=255)
                                                                                         aaaaCGGCCTGCGGGCCGTCATGTTTTGTCTGGATGTTGCGCTCTGTGCACAGCATTCAGACGACAGTTGTGCGGGAGTTGTGAAAAAATACGCGACa           >  1:1971825/1‑98 (MQ=255)
                                                                                         aaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgc         >  2:1189964/1‑100 (MQ=255)
                                                                                          aaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgcg        >  2:1854462/1‑100 (MQ=255)
                                                                                             cGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgcgcgc     >  2:2704723/1‑100 (MQ=255)
                                                                                              ggCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCa    >  1:2588835/1‑100 (MQ=255)
                                                                                              ggCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCa    >  2:2613788/1‑100 (MQ=255)
                                                                                               gCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCaa   <  2:548190/100‑1 (MQ=255)
                                                                                               gCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCaa   >  2:1036454/1‑100 (MQ=255)
                                                                                                ccTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAAt  >  1:812251/1‑100 (MQ=255)
                                                                                                  |                                                                                                 
CTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAAT  >  NC_000913/454755‑454950

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: