Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 454,853 T→C 17.2% 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,8530TC17.2% 223.1 / 40.5 105intergenic (+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 (50/37);  new base C (8/10);  total (58/47)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

Only 100 of 106 total aligned reads displayed.

CTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATA  >  NC_000913/454755‑454951
                                                                                                  |                                                                                                  
ctcctcCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTg                                                                                                   <  2:384716/100‑1 (MQ=255)
 tcctccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGc                                                                                                  <  1:303473/100‑1 (MQ=255)
 tcctccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGGCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTg                                                                                                   >  1:2242498/1‑99 (MQ=255)
    tccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCggg                                                                                               <  2:76418/100‑1 (MQ=255)
    tccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCggg                                                                                               <  2:209490/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:311782/1‑100 (MQ=255)
            gTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttg                                                                                       <  2:1078813/100‑1 (MQ=255)
            gTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttg                                                                                       <  2:1270030/100‑1 (MQ=255)
            gTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgttttt                                                                                       <  1:556639/100‑2 (MQ=255)
            gTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgttttt                                                                                       <  1:117645/100‑2 (MQ=255)
               gAGCAGGAATCATCGCGTAAAAACGCATTTGTAACTGTCGGCGCTTTGCCAGTATGTTGCTAAAGATTTTAGGAAAAACGGCCCGCAGGCCgtttttcat                                                                                    <  1:2758933/100‑5 (MQ=37)
               gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgttt                                                                                    <  2:1891953/100‑1 (MQ=255)
               gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                    <  2:2024270/100‑5 (MQ=39)
               gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                    <  2:525479/100‑5 (MQ=39)
               gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                    <  2:1913556/100‑5 (MQ=39)
               gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                    <  2:784335/100‑5 (MQ=255)
               gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                    <  1:2071483/100‑5 (MQ=39)
                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)
                    ggAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTCTGTTGCTAAAGATTTTATGAAAAACGGCATGCGGGCCGTTTTGTTTTGTCt                                                                               >  1:694988/1‑100 (MQ=255)
                    ggAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCt                                                                               <  1:1854112/100‑1 (MQ=255)
                     gAAGTATCCCTTAAAAACGCATTTGCAACTGTCGGCACCTTTCCAGTAGCTTGCTAAAGATTTTATGAAAGACGGTCCGCAGGTCgtttttc                                                                                      <  1:1034746/92‑3 (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:1033522/100‑1 (MQ=255)
                            cgcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTg                                                                       >  1:2240558/1‑100 (MQ=255)
                            cgcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTg                                                                       <  2:2606927/100‑1 (MQ=255)
                            cgcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTg                                                                       >  1:1234904/1‑100 (MQ=255)
                             gcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgc                                                                      >  2:2123585/1‑100 (MQ=255)
                             gcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgc                                                                      >  2:593368/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:2460457/1‑100 (MQ=255)
                                                aCTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCa                                                   >  1:2386896/1‑100 (MQ=255)
                                                aCTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCa                                                   >  1:521349/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:1894753/100‑1 (MQ=255)
                                                                                   ttATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcgc                >  1:451390/1‑100 (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:59174/1‑100 (MQ=255)
                                                                                      tGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc             >  2:890973/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)
                                                                                                 cTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATa  <  1:1906324/100‑1 (MQ=255)
                                                                                                 cTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATa  <  1:1010930/100‑1 (MQ=255)
                                                                                                  tGCGGGCCGTTTTGTTTTGTCTGGATTTTGCAATTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGGGCaaaa  >  1:83417/1‑97 (MQ=255)
                                                                                                  |                                                                                                  
CTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATA  >  NC_000913/454755‑454951

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: