Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,331,039 A→C intergenic (‑371/‑9) yciN ← / → topA DUF2498 protein YciN/DNA topoisomerase I, omega subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,331,0390AC97.1% 377.4 / 1.2 105intergenic (‑371/‑9)yciN/topADUF2498 protein YciN/DNA topoisomerase I, omega subunit
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base C (48/54);  total (49/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

Only 100 of 105 total aligned reads displayed.

GGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCA  >  NC_000913/1330942‑1331135
                                                                                                 |                                                                                                
gggTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAg                                                                                                >  1:861756/1‑100 (MQ=255)
   ttGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTg                                                                                             <  1:1723632/100‑1 (MQ=255)
      atctCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAat                                                                                          <  2:1213970/97‑1 (MQ=255)
        atCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAatat                                                                                        <  2:1815181/100‑1 (MQ=255)
        atCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAatat                                                                                        >  2:2227747/1‑100 (MQ=255)
         tCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATg                                                                                       <  2:203744/100‑1 (MQ=255)
         tCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATg                                                                                       <  1:1336096/100‑1 (MQ=255)
          ccGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAatagat                                                                                      >  1:1213966/1‑97 (MQ=255)
          ccGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATgg                                                                                      >  1:140509/1‑100 (MQ=255)
           cttccgatctgAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATggg                                                                                     <  2:2403928/90‑1 (MQ=255)
           cGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATggg                                                                                     <  1:721535/100‑1 (MQ=255)
           cGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATggg                                                                                     >  2:287373/1‑100 (MQ=255)
             cAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTa                                                                                   >  1:581899/1‑100 (MQ=255)
                 gagCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGc                                                                               <  1:1934678/100‑1 (MQ=255)
                    cGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGctct                                                                            >  1:1211040/1‑100 (MQ=255)
                    cGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGctct                                                                            >  1:452330/1‑100 (MQ=255)
                    cGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGctct                                                                            >  1:452317/1‑100 (MQ=255)
                    cGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGctct                                                                            >  2:674454/1‑100 (MQ=255)
                     gAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATgggagatcggaag                                                                           >  1:2403910/1‑90 (MQ=255)
                     gAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCtt                                                                           <  1:590107/100‑1 (MQ=255)
                        tCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTc                                                                        <  2:311147/100‑1 (MQ=255)
                        tCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTc                                                                        <  1:822313/100‑1 (MQ=255)
                         ccATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCCACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCa                                                                       >  2:1829168/1‑100 (MQ=255)
                         ccATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCa                                                                       <  2:1211044/100‑1 (MQ=255)
                         ccATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCa                                                                       <  2:452331/100‑1 (MQ=255)
                         ccATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCa                                                                       <  2:452318/100‑1 (MQ=255)
                          cATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCAt                                                                      >  2:125455/1‑100 (MQ=255)
                          cATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCAt                                                                      >  2:1928698/1‑100 (MQ=255)
                            tatCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCg                                                                    <  1:973872/100‑1 (MQ=255)
                            tatCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCg                                                                    >  1:570545/1‑100 (MQ=255)
                            tatCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCg                                                                    >  1:1643327/1‑100 (MQ=255)
                               cGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTg                                                                 <  1:1341684/100‑1 (MQ=255)
                               cGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTg                                                                 <  2:861759/100‑1 (MQ=255)
                               cGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTg                                                                 <  1:1068603/100‑1 (MQ=255)
                                ggTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGa                                                                >  2:1046486/1‑100 (MQ=255)
                                 gTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAg                                                               >  1:405732/1‑100 (MQ=255)
                                 gTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAg                                                               <  1:752884/100‑1 (MQ=255)
                                 gTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAg                                                               <  2:800310/100‑1 (MQ=255)
                                  tAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGt                                                              >  2:1398779/1‑100 (MQ=255)
                                   aaCTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTc                                                             <  1:767433/100‑1 (MQ=255)
                                    aCTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTcc                                                            <  2:2231304/100‑1 (MQ=255)
                                       cGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCg                                                         >  2:2225013/1‑100 (MQ=255)
                                          tGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCa                                                      >  1:1027595/1‑100 (MQ=255)
                                          tGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCa                                                      <  2:530444/100‑1 (MQ=255)
                                            ncAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCaaa                                                    >  1:932479/2‑100 (MQ=255)
                                            ccAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCaaa                                                    >  2:1421158/1‑100 (MQ=255)
                                            ccAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCaaa                                                    >  1:2139892/1‑100 (MQ=255)
                                              aGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAg                                                  <  1:168888/100‑1 (MQ=255)
                                              aGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAg                                                  >  2:168888/1‑100 (MQ=255)
                                              aGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAg                                                  <  1:2227733/100‑1 (MQ=255)
                                              aGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATCGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCaaaan                                                  <  1:1396425/100‑2 (MQ=255)
                                              aGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATCGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAg                                                  <  1:1396614/100‑1 (MQ=255)
                                                  gAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCaa                                              >  1:2335613/1‑100 (MQ=255)
                                                   aaGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTGCCCGGCAAAAGCCaaa                                             >  1:1158255/1‑100 (MQ=255)
                                                     ggTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAAc                                           <  2:1080937/100‑1 (MQ=255)
                                                      gTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACg                                          <  1:1090838/100‑1 (MQ=255)
                                                      gTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACg                                          >  1:182404/1‑100 (MQ=255)
                                                      gTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACg                                          <  2:50563/100‑1 (MQ=255)
                                                        ttATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGAt                                        >  2:2386473/1‑100 (MQ=255)
                                                          ttCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATca                                      <  1:1453337/99‑1 (MQ=255)
                                                          aTCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATca                                      <  2:2425807/100‑1 (MQ=255)
                                                           tCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATcaa                                     >  2:1453344/1‑100 (MQ=255)
                                                           tCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATcaa                                     >  2:518670/1‑100 (MQ=255)
                                                            cAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATcaac                                    >  1:244602/1‑100 (MQ=255)
                                                               cGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATcaaca                                   <  2:109762/98‑1 (MQ=255)
                                                               cGTGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAg                                 <  1:488124/100‑1 (MQ=255)
                                                                 tGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTa                               <  1:1592266/100‑1 (MQ=255)
                                                                 tGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTa                               >  2:787212/1‑100 (MQ=255)
                                                                 tGCGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTa                               <  2:1027598/100‑1 (MQ=255)
                                                                   cGACGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATc                             >  2:1600771/1‑100 (MQ=255)
                                                                     aCGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTg                           >  2:2487305/1‑100 (MQ=255)
                                                                     aCGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTg                           <  2:1738086/100‑1 (MQ=255)
                                                                      cGCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTgg                          <  2:251409/100‑1 (MQ=255)
                                                                       gCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTggg                         <  1:518669/100‑1 (MQ=255)
                                                                       gCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTggg                         >  1:427808/1‑100 (MQ=255)
                                                                       gCATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTggg                         <  1:1928688/100‑1 (MQ=255)
                                                                        cATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGt                        >  1:1202040/1‑100 (MQ=255)
                                                                        cATTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGt                        <  2:753372/100‑1 (MQ=255)
                                                                         aTTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTa                       >  2:1040547/1‑100 (MQ=255)
                                                                         aTTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTa                       >  2:418425/1‑100 (MQ=255)
                                                                         aTTCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTa                       >  2:250281/1‑100 (MQ=255)
                                                                          ttCCTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAg                      <  1:1491728/100‑1 (MQ=255)
                                                                            ccTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTg                    <  1:1046483/100‑1 (MQ=255)
                                                                            ccTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTg                    <  2:2335631/100‑1 (MQ=255)
                                                                             cTGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGa                   >  1:2206719/1‑100 (MQ=255)
                                                                              tGGAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGAc                  <  1:2574288/100‑1 (MQ=255)
                                                                               ggAAGAATCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACt                 >  1:489845/1‑100 (MQ=255)
                                                                                 aagaaTCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTAc               <  2:870740/100‑1 (MQ=255)
                                                                                  agaaTCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACg              <  2:140509/100‑1 (MQ=255)
                                                                                   gaaTCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACgt             >  1:1817897/1‑100 (MQ=255)
                                                                                    aaTCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACgtg            >  2:1663300/1‑100 (MQ=255)
                                                                                     aTCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACgtgg           <  2:427809/100‑1 (MQ=255)
                                                                                     aTCAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACgtg            >  1:1329981/1‑99 (MQ=255)
                                                                                     aTCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACgtgg           >  1:694477/1‑100 (MQ=255)
                                                                                      acAAATTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACgtggt          <  1:1544765/99‑1 (MQ=255)
                                                                                         aaTTAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGaa       <  2:489846/100‑1 (MQ=255)
                                                                                           ttAGGTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAAt     <  2:2206732/100‑1 (MQ=255)
                                                                                              ggTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCa  >  1:2476499/1‑100 (MQ=255)
                                                                                              ggTCAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCa  >  1:1203332/1‑100 (MQ=255)
                                                                                              ggTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCa  <  2:27598/100‑1 (MQ=255)
                                                                                                 |                                                                                                
GGGTTGATATCCGCAGAGAGCGAGTCCATATCGGTAACTCGTTGCCAGTGGAAGGTTTATCAACGTGCGACGCATTCCTGGAAGAATCAAATTAGGTAAGGTGAATATGGGTAAAGCTCTTGTCATCGTTGAGTCCCCGGCAAAAGCCAAAACGATCAACAAGTATCTGGGTAGTGACTACGTGGTGAAATCCA  >  NC_000913/1330942‑1331135

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