Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,822,089 T→A 26.8% intergenic (+38/‑17) gmk → / → rpoZ guanylate kinase/RNA polymerase, omega subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,822,0890TA26.8% 296.1 / 129.2 177intergenic (+38/‑17)gmk/rpoZguanylate kinase/RNA polymerase, omega subunit
Reads supporting (aligned to +/- strand):  ref base T (65/64);  new base A (29/18);  total (95/82)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01

Only 100 of 179 total aligned reads displayed.

GCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCCGCGCGTCGCGCTCG  >  NC_000913/3821993‑3822188
                                                                                                |                                                                                                   
gcATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCagtg                                                                                                  <  1:1286528/100‑1 (MQ=255)
   tGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAg                                                                                               <  2:424879/100‑1 (MQ=255)
     aGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCt                                                                                             >  2:3057312/1‑100 (MQ=255)
      gccgccAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCTGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCtt                                                                                            >  1:2836115/1‑100 (MQ=255)
      gccgccAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCtt                                                                                            >  1:1149761/1‑100 (MQ=255)
      gccgccAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCtt                                                                                            >  1:3856513/1‑100 (MQ=255)
      gccgccAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCtt                                                                                            >  2:1080072/1‑100 (MQ=255)
      gccgccAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCtt                                                                                            >  2:2543178/1‑100 (MQ=255)
        cgccAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTGCACCTGTGGAGCtttt                                                                                          <  1:1460764/100‑1 (MQ=255)
         gccAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCttttt                                                                                         >  2:2116757/1‑100 (MQ=255)
         gccAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCttttt                                                                                         >  2:2188939/1‑100 (MQ=255)
         gcaaaaaGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAGCCTGATTTCAGTATCATGCCCAGGCATTTCTTCACCTGTGGAGCtttt                                                                                          <  1:2515570/96‑1 (MQ=255)
          ccAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTa                                                                                        <  2:3498218/100‑1 (MQ=255)
           cAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTaa                                                                                       >  1:1853173/1‑100 (MQ=255)
           cAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTaa                                                                                       <  1:743428/100‑1 (MQ=255)
            aaacGCAGCGTCATGACGCTTTAATCAGCACATTGTTGGCAGGCTGAACCCGATTTCAGTATCATGCCCCGCCTTTTCTTCACCCGCGGAGCTTTTGAAg                                                                                      >  1:1941081/1‑100 (MQ=255)
            aaaaGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAg                                                                                      <  2:3040627/100‑1 (MQ=255)
            aaaaGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAg                                                                                      <  2:824777/100‑1 (MQ=255)
            aaaaGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAg                                                                                      <  2:3977071/100‑1 (MQ=255)
             aaaGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGt                                                                                     <  1:2983990/100‑1 (MQ=255)
                gcagcGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATg                                                                                  >  2:2631247/1‑100 (MQ=255)
                   gcGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCa                                                                               <  1:66630/100‑1 (MQ=255)
                   gcGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCAGTGGAGCTTTTTAAGTATGGCa                                                                               >  1:3474362/1‑100 (MQ=255)
                    cGTCATGACGCGTTAATCAGCAAATTGTTGGAACACTGAACCTGATTTCAGTATCATGCCCAGTCATCTCTTCACCTGTGGAGCTTTTTAAGTATGGCAc                                                                              >  1:2247071/1‑100 (MQ=255)
                      tCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCAcgc                                                                            >  2:3386677/1‑100 (MQ=255)
                       cttGACGCTTTACTTGGCAAATTGTTGGCAGACTGAACCTGTTTTCAGTATTATGCCCACTTATTTCTTCACCTGTGGTGGTTTTTACGTATggtag                                                                              >  1:320491/3‑94 (MQ=255)
                       cATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCAcgcg                                                                           <  2:2361184/100‑1 (MQ=255)
                         tGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCAGTGGAGCTTTTTAAGTATGGCACGCGTa                                                                         >  1:644363/1‑100 (MQ=255)
                           aCGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAAc                                                                       <  2:1731987/100‑1 (MQ=255)
                           aCGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCAGTGGAGCTTTTTAAGTATGGCACGCGTAAc                                                                       <  2:1658098/100‑1 (MQ=255)
                           aCGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCAGTGGAGCTTTTTAAGTATGGCACGCGTAAc                                                                       >  2:154168/1‑100 (MQ=255)
                               tttAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGt                                                                    >  1:142217/1‑99 (MQ=255)
                                ttAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTc                                                                  <  1:1750591/100‑1 (MQ=255)
                                ttAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCAGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTc                                                                  <  1:706495/100‑1 (MQ=255)
                                ttAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCAGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTc                                                                  <  2:3094765/100‑1 (MQ=255)
                                 tAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCa                                                                 <  2:1100243/100‑1 (MQ=255)
                                  aaTCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAg                                                                <  2:1157699/100‑1 (MQ=255)
                                  aaTCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAg                                                                <  2:2836126/100‑1 (MQ=255)
                                   aTCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAgg                                                               <  2:1125451/100‑1 (MQ=255)
                                      aGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTTTTTTTTGTACGCGTAACTGTTCAGGAc                                                             >  1:374874/1‑99 (MQ=255)
                                      aGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACg                                                            >  2:1000268/1‑100 (MQ=255)
                                       gCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGGCCGGTCATTTCTTCACCTGTGGAGCTTTTTAGGTATGGCACGCGTAACTTTTTAGGACGc                                                           >  1:857526/1‑100 (MQ=255)
                                       gCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCAGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGc                                                           >  1:3363732/1‑100 (MQ=255)
                                        cAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGc                                                           <  1:1719466/99‑1 (MQ=255)
                                        cAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCt                                                          >  1:4115723/1‑100 (MQ=255)
                                        cAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCt                                                          >  1:4115778/1‑100 (MQ=255)
                                         aaaTTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGc                                                           >  2:3240801/1‑98 (MQ=255)
                                            ttgttgGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGt                                                        <  1:546083/98‑1 (MQ=255)
                                             tgttgGCAGCCTGAACCTGGTTTCAGTATCATGCTCTGTCTTTTCTTCACCTGTTGCGCTTTCTAATTGTGGCACGCGTAACTGTTCAGGAGGCTGTag                                                      >  1:2141307/1‑99 (MQ=255)
                                               ttgGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCAGTGGAGCTTTTTTTTTTTTTCTCTCGTAACTGTTCAGGACGCTGTAgaga                                                   >  1:1559723/1‑100 (MQ=255)
                                                tgGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGaa                                                  >  2:2974645/1‑100 (MQ=255)
                                                 ggCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGaaa                                                 >  1:3607508/1‑100 (MQ=255)
                                                 ggCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCAGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGaaa                                                 >  1:2513970/1‑100 (MQ=255)
                                                   cAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAAt                                               <  1:610255/100‑1 (MQ=255)
                                                   cAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCAGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAAt                                               >  2:526599/1‑100 (MQ=255)
                                                    aGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAAtt                                              >  1:2178482/1‑100 (MQ=255)
                                                      aCTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTgg                                            <  1:1023356/100‑1 (MQ=255)
                                                      aCTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCAGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTgg                                            <  2:2318835/100‑1 (MQ=255)
                                                       cTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGt                                           <  2:3552701/100‑1 (MQ=255)
                                                       cTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCAGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGt                                           >  1:2429626/1‑100 (MQ=255)
                                                          aaCCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAAc                                        <  2:393158/100‑1 (MQ=255)
                                                           aCCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAAcc                                       >  1:2860113/1‑100 (MQ=255)
                                                           aCCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAAcc                                       <  1:2627277/100‑1 (MQ=255)
                                                            ccTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCg                                      >  2:492392/1‑100 (MQ=255)
                                                            ccTGATTTCAGTATCATGCCCAGTCATTTCTTCACCAGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCg                                      >  2:2953311/1‑100 (MQ=255)
                                                              cgacttCAGTATCATGCCCAGTCATTTCTTCACCAGTAGAGCTTTCTAAGTAGGCCACACGTAACTGTTCCGCAAGCTGCAGAGCAAATAGGt                                           >  1:3185694/5‑93 (MQ=255)
                                                                aTTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGtttt                                  <  1:1083133/100‑1 (MQ=255)
                                                                aTTTCAGTATCATGCCCAGTCATTTCTTCACCAGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAAGTGGTAACCGtttt                                  <  1:451244/100‑1 (MQ=255)
                                                                 tttCAGTATCATGCCCAGTCATTTCTTCACCAGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGtttt                                  <  1:806053/99‑1 (MQ=255)
                                                                 tttCAGTATCATGCCCAGTCATTTCTTCACCAGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGAAGCTGTAGAGAAAATTGGTAACCGTTTTg                                 >  1:576876/1‑100 (MQ=255)
                                                                  ttCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGa                                <  1:3945447/100‑1 (MQ=255)
                                                                   tCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGAc                               <  2:2685044/100‑1 (MQ=255)
                                                                    cAGTATCATGCCCAGTCATTTCTTCACCAGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGAcc                              >  2:134227/1‑100 (MQ=255)
                                                                     aGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGAcc                              <  2:3596023/99‑1 (MQ=255)
                                                                      gTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTg                            >  1:3550468/1‑100 (MQ=255)
                                                                      gTATCATGCCCAGTCATTTCTTCACCAGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTg                            <  1:3797023/100‑1 (MQ=255)
                                                                      gTATCAGGCCCAGTCATTTCTTCACCAGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTg                            <  2:722627/100‑1 (MQ=255)
                                                                       tATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTgg                           <  1:4225134/100‑1 (MQ=255)
                                                                       tATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTgg                           <  1:1530841/100‑1 (MQ=255)
                                                                       tATCATGCCCAGTCATTTCTTCACCAGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTgg                           <  1:2958688/100‑1 (MQ=255)
                                                                        aTCATGCCAAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGt                          >  1:409356/1‑100 (MQ=255)
                                                                          cATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCCGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTAc                        <  1:2326187/100‑1 (MQ=255)
                                                                          cATGCCCAGTCATTTCTTCACCAGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTan                        <  2:843972/100‑2 (MQ=255)
                                                                           aTGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACt                       >  1:2857159/1‑100 (MQ=255)
                                                                            tGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTg                      >  1:3961675/1‑100 (MQ=255)
                                                                            tGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGAGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTg                      >  2:2459649/1‑100 (MQ=255)
                                                                             gCCCAGTCATTTCTTCACCAGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTgg                     <  2:3817713/100‑1 (MQ=255)
                                                                              cccAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGt                    >  2:3761487/1‑100 (MQ=255)
                                                                              cccAGTCATTTCTTCACCAGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGt                    >  2:4065554/1‑100 (MQ=255)
                                                                                  gTCATTTCTTCACCAGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTcgcc                >  2:1344835/1‑100 (MQ=255)
                                                                                    cATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTcgccgc              >  1:3066613/1‑100 (MQ=255)
                                                                                     aTTTCTTCACCAGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCcgcg             >  2:2962853/1‑100 (MQ=255)
                                                                                     aTTTCTTCACCAGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCcgcg             >  2:2139554/1‑100 (MQ=255)
                                                                                       ttagtcACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACGGTTTAGGACACGGTAGAGAAAATTGGGAACCGTTTTGACCTGGTACTGGTcgcagc              >  1:1997147/5‑97 (MQ=255)
                                                                                        tcttcACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCCGCGCGt          >  1:2618450/1‑100 (MQ=255)
                                                                                        tcttcACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCCGCGCGt          <  1:1423811/100‑1 (MQ=255)
                                                                                          ttcACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCCGCGCGTcg        <  1:3075917/100‑1 (MQ=255)
                                                                                          ttcACCAGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCCGCGCGTcg        >  2:3460836/1‑100 (MQ=255)
                                                                                            cACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCCGCGCGTcgcg      <  2:583166/100‑1 (MQ=255)
                                                                                            cACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCCGCGCGTcgcg      >  1:514318/1‑100 (MQ=255)
                                                                                                |                                                                                                   
GCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCCGCGCGTCGCGCTCG  >  NC_000913/3821993‑3822188

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: