Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 454,856 G→A 16.9% intergenic (+67/‑277) 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,8560GA16.9% 229.5 / 32.3 102intergenic (+67/‑277)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 G (50/33);  new base A (7/10);  total (59/43)
Fisher's exact test for biased strand distribution p-value = 1.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

Only 100 of 104 total aligned reads displayed.

TCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAAC  >  NC_000913/454759‑454953
                                                                                                 |                                                                                                 
tccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCggg                                                                                                 <  2:209490/100‑1 (MQ=255)
tccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCggg                                                                                                 <  2:76418/100‑1 (MQ=255)
  ccTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGcc                                                                                               >  2:2004442/1‑100 (MQ=255)
   cTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCg                                                                                              >  2:1738111/1‑100 (MQ=255)
     gTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgtt                                                                                            >  2:5668/1‑100 (MQ=255)
     gTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgtt                                                                                            >  1:2242219/1‑100 (MQ=255)
     gTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgtt                                                                                            >  1:41137/1‑100 (MQ=255)
        gTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttg                                                                                         <  2:1078813/100‑1 (MQ=255)
        gTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttg                                                                                         <  2:1270030/100‑1 (MQ=255)
        gTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttg                                                                                         >  2:311782/1‑100 (MQ=255)
        gTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgttttt                                                                                         <  1:117645/100‑2 (MQ=255)
        gTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgttttt                                                                                         <  1:556639/100‑2 (MQ=255)
         tGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttc                                                                                        <  2:26117/100‑3 (MQ=255)
           gAGCAGGAATCATCGCGTAAAAACGCATTTGTAACTGTCGGCGCTTTGCCAGTATGTTGCTAAAGATTTTAGGAAAAACGGCCCGCAGGCCgtttttcat                                                                                      <  1:2758933/100‑5 (MQ=37)
           gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgttt                                                                                      <  2:1891953/100‑1 (MQ=255)
           gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                      <  2:1913556/100‑5 (MQ=39)
           gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                      <  2:525479/100‑5 (MQ=39)
           gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                      <  1:2071483/100‑5 (MQ=39)
           gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                      <  2:784335/100‑5 (MQ=255)
            aGCAGGAAGCAGCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTGTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtttt                                                                                     <  2:255039/100‑1 (MQ=255)
             gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtttt                                                                                     >  1:2820932/1‑99 (MQ=255)
             gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGttttgttttg                                                                                    >  1:1027641/1‑100 (MQ=255)
             gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTCTATGAAAAACGGCCTGCGGGCCGttttgttttg                                                                                    >  1:1340325/1‑100 (MQ=255)
             gCAGGAAGCATCGCGTAAAAACACATTTGCAACTGTCGGCGCTGTTCAGGTATGTTGCTAAAGATGTTATGAAAAACGTCCTGCCGGCCgttttgcttt                                                                                     >  1:167372/1‑99 (MQ=255)
              cAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTtttgttttgt                                                                                   >  1:1802780/1‑100 (MQ=255)
              cAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTtttgttttgt                                                                                   <  2:560108/100‑1 (MQ=255)
               aGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTc                                                                                  <  1:1178317/100‑1 (MQ=255)
                ggAAGCATCGCGTAAAATCGCCTTTGCGACTGTCGGCGCGTTTCTAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCt                                                                                 <  2:867412/100‑1 (MQ=255)
                ggAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCt                                                                                 <  1:1854112/100‑1 (MQ=255)
                ggAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCt                                                                                 >  1:1678018/1‑100 (MQ=255)
                 gAAGTATCCCTTAAAAACGCATTTGCAACTGTCGGCACCTTTCCAGTAGCTTGCTAAAGATTTTATGAAAGACGGTCCGCAGGTCgtttttc                                                                                        <  1:1034746/92‑3 (MQ=255)
                  aaGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTg                                                                                <  1:690562/99‑1 (MQ=255)
                   aGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGa                                                                              <  2:1831174/100‑1 (MQ=255)
                   aGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGa                                                                              <  1:2868339/100‑1 (MQ=255)
                       tCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGAtttt                                                                          <  1:2870842/100‑1 (MQ=255)
                        cgcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTg                                                                         >  1:1234904/1‑100 (MQ=255)
                        cgcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTg                                                                         <  2:2606927/100‑1 (MQ=255)
                        cgcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTg                                                                         <  1:1033522/100‑1 (MQ=255)
                         gcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgc                                                                        >  2:593368/1‑100 (MQ=255)
                         gcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgc                                                                        >  2:2123585/1‑100 (MQ=255)
                              aaaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTTCGGGCCGTTTTGTTTTTTCTGGATTTTTCGCttt                                                                   >  1:1176145/1‑100 (MQ=255)
                                aaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCATTTTGTTTTGTCTGGATTTTGCGCttttt                                                                 >  2:2903719/1‑100 (MQ=255)
                                   gCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGcc                                                              >  1:2340477/1‑100 (MQ=255)
                                     aTTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCa                                                            <  2:2145428/100‑1 (MQ=255)
                                      tttGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAg                                                           <  1:1825682/100‑1 (MQ=255)
                                      tttGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAg                                                           <  1:2658367/100‑1 (MQ=255)
                                         gccaCTGTCGGCGCCTTTCCCGTATGTTGCTAAAGAGTTTATGAAAAACGGCCTGCGGGCCGTTTTTTTTTTTTTGGATTTTGCGCTTTTTGccc                                                             >  1:2614857/4‑95 (MQ=255)
                                         gCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCAt                                                        >  2:1930890/1‑100 (MQ=255)
                                            aCTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCTCTGCAtt                                                       >  1:1347748/1‑98 (MQ=255)
                                            aCTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCa                                                     >  1:521349/1‑100 (MQ=255)
                                            aCTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCa                                                     >  1:2460457/1‑100 (MQ=255)
                                            aCTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCa                                                     >  1:2386896/1‑100 (MQ=255)
                                             cTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAg                                                    <  1:1099238/100‑1 (MQ=255)
                                             cTGTAGGCGCTTTTGCATTATCTTTCTATAGATTTTATGAAAATCGGCCTGCGTGCCTTTTTGTTTTTTCTGGATTTTGCGCTTTTTGCCCAGCa                                                         <  1:2872683/95‑1 (MQ=255)
                                                  ggCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaa                                               >  2:1409498/1‑100 (MQ=255)
                                                   gcgcTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaaa                                              >  1:718953/1‑100 (MQ=255)
                                                           ccAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCg                                      <  1:792667/100‑1 (MQ=255)
                                                             aGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCTCggg                                    >  2:1743200/1‑100 (MQ=255)
                                                             aGTATGGTGCTAACGATTTTATGAAAAACGGCCTGCGGGCCGTGTTGTGTTGTCTGGAGTTTGCGCTTTTTGCCCAGCATTCAGACGAAAACTGCCCggg                                    <  1:1232046/100‑1 (MQ=255)
                                                              gTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGa                                   >  2:1481621/1‑100 (MQ=255)
                                                                 tGTTGCTAAAGATTTTATGAAAAACGTCCTGCGGGCCGTTTTGTTTTGTCTGGATATTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAAt                                 >  1:2594741/1‑99 (MQ=255)
                                                                  gTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtg                               <  1:974569/100‑1 (MQ=255)
                                                                   ttGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtgt                              >  2:1935493/1‑100 (MQ=255)
                                                                    tGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtgtg                             >  2:2104913/1‑100 (MQ=255)
                                                                     gCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGa                            <  1:697666/100‑1 (MQ=255)
                                                                       tAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaa                          >  2:168863/1‑100 (MQ=255)
                                                                        aaaGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaaa                         >  2:1606083/1‑100 (MQ=255)
                                                                        aaaGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaaa                         <  1:1259072/100‑1 (MQ=255)
                                                                         aaGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCGCGGGAATTGTGaaaaa                        >  1:1703954/1‑100 (MQ=255)
                                                                            agacaaaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATa                     >  1:1661823/10‑100 (MQ=38)
                                                                            aTTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCTGGGAATTGTGAAAAAATa                     >  1:2682778/1‑100 (MQ=255)
                                                                            aTTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATa                     >  2:2862548/1‑100 (MQ=255)
                                                                             ttttATGACAACCGGCTTTCGTGCCGTTTTGTTTTATCTACATTTTGCGCTTTTTACCCAGCATTCAGACGAAAATTGCCCGTGAATTGTGAAAAAATa                     <  1:2651931/99‑1 (MQ=255)
                                                                              tttATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcg                   >  1:2163500/1‑100 (MQ=255)
                                                                              attATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTGGCGCTTTTGGCCCAGCAGGCAGACGAAAATTGCGCGGTAATTGTGAAAAAATAc                    >  1:2432186/2‑99 (MQ=255)
                                                                               ttATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcgc                  <  2:1850690/100‑1 (MQ=255)
                                                                               ttATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcgc                  >  1:451390/1‑100 (MQ=255)
                                                                               ttATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcgc                  <  2:1894753/100‑1 (MQ=255)
                                                                               gaaaacaaaaCGGCCCGCTGGCCGTTTTGTTTTGTGTGGATTTTGCGCTTTTGGCCCAGCATGTAGACGAAAATTGCCCTGGGTTTGTGAACAAATAggc                  >  1:554027/7‑100 (MQ=255)
                                                                                tATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCCGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcgcg                 >  2:1424487/1‑100 (MQ=255)
                                                                                tATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcgcg                 >  1:2863513/1‑100 (MQ=255)
                                                                                 aaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGa                >  1:337570/5‑100 (MQ=39)
                                                                                 aaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGa                >  2:2913548/5‑100 (MQ=39)
                                                                                 aaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATCTTGCGCTTTTTCCCCAGCCTTCCGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGa                >  1:234667/5‑100 (MQ=37)
                                                                                  tGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc               >  2:890973/1‑100 (MQ=255)
                                                                                  tGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc               >  2:59174/1‑100 (MQ=255)
                                                                                  aacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc               >  2:2283904/4‑100 (MQ=255)
                                                                                   acaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACa              >  2:336685/3‑100 (MQ=255)
                                                                                    aaaaaCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAg             >  2:2096919/1‑100 (MQ=255)
                                                                                     aaaaCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc               >  1:957369/1‑97 (MQ=255)
                                                                                      aaaCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgcg           >  2:218064/1‑100 (MQ=255)
                                                                                         cGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgcgcgc        >  2:1906333/1‑100 (MQ=255)
                                                                                          ggCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTGGCCCAGCATTCAGACAAAAATTGCCCGGGAATTTTGAAAAAATACGCGACAgcgcg         <  1:143577/98‑1 (MQ=255)
                                                                                          ggCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCa       >  2:2685706/1‑100 (MQ=255)
                                                                                            ccTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgcgcn         <  1:1621199/96‑2 (MQ=255)
                                                                                            ccTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAAt     <  1:2636556/100‑1 (MQ=255)
                                                                                             cTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATa    <  1:1906324/100‑1 (MQ=255)
                                                                                             cTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATa    <  1:1010930/100‑1 (MQ=255)
                                                                                              tGCGGGCCGTTTTGTTTTGTCTGGATTTTGCAATTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGGGCaaaa    >  1:83417/1‑97 (MQ=255)
                                                                                               gCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAAc  >  2:1917218/1‑100 (MQ=255)
                                                                                                 |                                                                                                 
TCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAAC  >  NC_000913/454759‑454953

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