Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,332,822 G→A 16.9% intergenic (‑87/+40) yhbE ← / ← rpmA EamA family inner membrane putative transporter/50S ribosomal subunit protein L27

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,332,8220GA16.9% 260.1 / 41.1 117intergenic (‑87/+40)yhbE/rpmAEamA family inner membrane putative transporter/50S ribosomal subunit protein L27
Reads supporting (aligned to +/- strand):  ref base G (56/41);  new base A (4/16);  total (60/57)
Fisher's exact test for biased strand distribution p-value = 2.74e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.53e-02

Only 100 of 119 total aligned reads displayed.

TGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCA  >  NC_000913/3332725‑3332920
                                                                                                 |                                                                                                  
tgcttcTTCATACTTTATTCCCTGCAAATCGCATTTTCGTATTAACGTCGCACCATGCGCATTTTTACCGGTTATCTAATGTAAAAAGCCCCGCAACg                                                                                                    <  1:538053/98‑1 (MQ=255)
 gctgctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACgtgt                                                                                                 <  2:1506723/100‑1 (MQ=255)
    gctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGc                                                                                              >  2:2186773/1‑100 (MQ=255)
     ctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGTGCATTTTTACCGGTTATCGAATGTAAAATGCCCCGCAACGTGTTGCg                                                                                             <  2:691711/100‑1 (MQ=255)
     ctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCg                                                                                             <  1:643315/100‑1 (MQ=255)
      ttCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCgg                                                                                            <  2:2765645/100‑1 (MQ=255)
       tCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCgag                                                                                           >  2:643311/1‑98 (MQ=255)
        cATACCGTATTCCCTGCAAAGCGCATTTGCGGATTATCGTTGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCgggg                                                                                          <  1:1513138/100‑1 (MQ=255)
        cATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCgggg                                                                                          <  1:2886419/100‑1 (MQ=255)
        cATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCgggg                                                                                          <  2:811718/100‑1 (MQ=255)
        cATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCgggg                                                                                          <  2:2500928/100‑1 (MQ=255)
        cATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCgggg                                                                                          <  2:2376307/100‑1 (MQ=255)
        cATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCgggg                                                                                          <  2:313945/100‑1 (MQ=255)
         aTACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGc                                                                                         <  2:2025937/100‑1 (MQ=255)
         aTACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGc                                                                                         <  1:1657065/100‑1 (MQ=255)
         aTACCGTATTCCCTGCAAAGCGCATTTGCGGATTAAAGTAGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGc                                                                                         <  2:695091/100‑1 (MQ=255)
          tACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCt                                                                                        <  2:1620149/100‑1 (MQ=255)
          tACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGc                                                                                         <  1:1350601/99‑1 (MQ=255)
           aCCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGCtt                                                                                       <  2:2457995/100‑1 (MQ=255)
           aCCGTATGCCCTGCAAAGCGGATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCtt                                                                                       >  2:1766239/1‑100 (MQ=255)
            ccGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCttt                                                                                      >  1:3196/1‑100 (MQ=255)
               tATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCan                                                                                   <  2:1791391/100‑2 (MQ=255)
               tATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGCttttta                                                                                   <  1:384315/100‑4 (MQ=38)
                aTTCCCTGCAAAGCGCATTTGGGGATTAGCGTCGCACCAGGCGCATTTTTAACAGTTATCGAATGTAAAACGCCCCGCAACGTGTTGAGGGGATGTCATc                                                                                  <  2:2701930/100‑1 (MQ=255)
                   gccTGCAACGCGCATTTGCAGATTAACGTTGCGCCGTGCGCATTTTTACCGGTTATCGAATGTAAAAAGCGCCGGAACGTGTTGCGGGTCTTTGATGCg                                                                                >  2:2238685/2‑99 (MQ=255)
                    ccTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGtt                                                                              <  2:2644693/100‑1 (MQ=255)
                       gCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTAcc                                                                           >  1:932414/1‑100 (MQ=255)
                       gCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTAcc                                                                           >  1:1059265/1‑100 (MQ=255)
                       gCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTAcc                                                                           >  2:239090/1‑100 (MQ=255)
                         aaaGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCgg                                                                         <  2:1238458/100‑1 (MQ=255)
                         aaaGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCgg                                                                         <  1:2315319/100‑1 (MQ=255)
                         aaaGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCgg                                                                         >  1:565030/1‑100 (MQ=255)
                          aaGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCggg                                                                        <  1:1485704/100‑1 (MQ=255)
                             cgcATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGAcg                                                                     <  1:1863360/100‑1 (MQ=255)
                              gcATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGAcgc                                                                    <  2:255879/100‑1 (MQ=255)
                               cATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGAcgcg                                                                   >  1:751542/1‑100 (MQ=255)
                                 cttGCGTATTAGCGCCCCACCAGACGCTTGTGTAGCGGTTATCGAATAGAAAAAGCCCCGCAACGTGCTGCGGGGCTTTCATCCGCTACCGGGCcgg                                                                    <  2:1589400/96‑2 (MQ=255)
                                  ttGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGaaa                                                                >  1:1061886/1‑100 (MQ=255)
                                   tGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGTATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGaaa                                                                >  2:432268/1‑99 (MQ=255)
                                   tGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGaaaa                                                               >  1:2705261/1‑100 (MQ=255)
                                   tGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGaaaa                                                               <  2:2320357/100‑1 (MQ=255)
                                    gCGGATTAACGTCGCACCATGCACATTTTTACCGGTTATCAAATGTAAAAAGCCCCGCAACATTTTGCGGGGATTTCATCCGTTACCGAGACGCGaaaaa                                                              >  2:771590/1‑100 (MQ=255)
                                     cGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAAc                                                             >  1:174954/1‑100 (MQ=255)
                                      ggATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAAc                                                             <  2:2597424/99‑1 (MQ=255)
                                      ggATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAAc                                                             <  2:1504887/99‑1 (MQ=255)
                                       gATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACtt                                                           <  2:932408/100‑1 (MQ=255)
                                       gATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACtt                                                           >  1:1695297/1‑100 (MQ=255)
                                        aTTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTa                                                          >  2:363206/1‑100 (MQ=255)
                                         ttAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTAt                                                         <  1:2815166/100‑1 (MQ=255)
                                           aaCGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTc                                                       <  2:1020568/100‑1 (MQ=255)
                                           aaCGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTc                                                       >  2:283997/1‑100 (MQ=255)
                                            aCGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCa                                                      >  2:2894795/1‑100 (MQ=255)
                                             cGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAg                                                     >  1:635323/1‑100 (MQ=255)
                                              gTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGc                                                    <  2:1182625/100‑1 (MQ=255)
                                               tCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCt                                                   <  2:1119287/100‑1 (MQ=255)
                                               tCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCt                                                   >  1:1380500/1‑100 (MQ=255)
                                                 gCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTc                                                 >  2:1379203/1‑100 (MQ=255)
                                                  cACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCg                                                <  1:1087015/100‑1 (MQ=255)
                                                  cACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCg                                                <  1:1944772/100‑1 (MQ=255)
                                                   aCCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCACAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGa                                               <  2:2834442/100‑1 (MQ=255)
                                                        gcgcATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTg                                          >  2:2719600/1‑100 (MQ=255)
                                                         cgcATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGa                                         >  1:752510/1‑100 (MQ=255)
                                                         cgcATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGa                                         >  1:2002319/1‑100 (MQ=255)
                                                         cgcATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGATACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGa                                         >  1:76828/1‑100 (MQ=255)
                                                           cATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATa                                       >  2:679225/1‑100 (MQ=255)
                                                             tttttACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATaaa                                     >  1:2362024/1‑100 (MQ=255)
                                                               tttACCGGTTATCGAATGTAAAAAGCCCAGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGCTGCTGATAAAt                                    >  2:1304305/1‑99 (MQ=255)
                                                                   ccGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCTGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGATGATAAATTTACg                               >  1:1018196/1‑100 (MQ=255)
                                                                   ccGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACg                               <  2:1380487/100‑1 (MQ=255)
                                                                    cGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACgg                              <  2:2838378/100‑1 (MQ=255)
                                                                    cGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACgg                              >  2:2620078/1‑100 (MQ=255)
                                                                     ggTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGt                             <  1:2592823/100‑1 (MQ=255)
                                                                        tATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGtttt                          >  1:2416604/1‑100 (MQ=255)
                                                                         aTCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGttttt                         <  2:1636464/100‑1 (MQ=255)
                                                                         aTCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACACGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGttttt                         >  1:2432553/1‑100 (MQ=255)
                                                                         aTCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGATACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGtttt                          >  2:722476/1‑99 (MQ=255)
                                                                           cGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCg                       >  1:2209427/1‑100 (MQ=255)
                                                                            gAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGttctt                         >  2:2512575/1‑97 (MQ=255)
                                                                            gAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCgg                      >  1:1763189/1‑100 (MQ=255)
                                                                            gAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTACGATGCTGATAAATTTACGGTTTTTCgg                      >  1:171855/1‑100 (MQ=255)
                                                                              aTGTAAAAAGCCCCGCAACGTGTTGCGGTGCTTTCATCCGTTACCGTGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCggg                     <  1:1484321/99‑1 (MQ=255)
                                                                                  aaaaaGCCCCGCAACGTGTTGAGGGGAATACATCCGTTCCGGGGCCGCGAAAAACTTATTCAGCTGCAGTGTTGATAAATTTACGGTTTTTCGGGcc                   >  2:2595584/1‑97 (MQ=255)
                                                                                   aaaattcccGCAACGGGTTGCTTGGCGTTCAGCTGTTACCGGGACTCGAATGACTTATTCTGCTTTGATGCTCATAAATTTACGGTTTTTCGGGGCttt                <  2:555257/93‑1 (MQ=255)
                                                                                   aaaaGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTa               >  1:2797074/1‑100 (MQ=255)
                                                                                     aaGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAAc             <  1:284000/100‑1 (MQ=255)
                                                                                     aaGCCCCGCAACGTGTTGCGAGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAAc             >  2:126075/1‑100 (MQ=255)
                                                                                     aaGCCCCGCAACACGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTaa              >  2:1546571/1‑99 (MQ=255)
                                                                                      aGCCCCGCAACGTGTTGCGCGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGAAGCTGATAAATTTAAGAAATTTAGAGCCTTAAACt            >  1:905654/1‑100 (MQ=255)
                                                                                      aGCCCCGCAACACGTTGCGGGGCTTTCATCCGTTACCGGGACGCTAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACt            <  1:801473/100‑1 (MQ=39)
                                                                                       gCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACtt           >  1:2460053/1‑100 (MQ=255)
                                                                                       gCCCCGCAACACGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACt            >  1:56020/1‑99 (MQ=255)
                                                                                        ccccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTc          >  1:2243262/1‑100 (MQ=255)
                                                                                         cccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGTGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCg         >  2:217607/1‑100 (MQ=255)
                                                                                         cccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCg         <  2:1747000/100‑1 (MQ=255)
                                                                                         cccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCg         >  2:767918/1‑100 (MQ=255)
                                                                                          ccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGa        <  2:2386299/100‑1 (MQ=255)
                                                                                           cGCAACGTGTTGCGGGTCTTTCATCCGTTACCGTGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGaa       <  1:2816761/100‑1 (MQ=255)
                                                                                               aCGTGTTGCGGGGATTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTc   >  2:793430/1‑100 (MQ=255)
                                                                                                cGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAAttgc   <  2:1554021/99‑3 (MQ=255)
                                                                                                cGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCa  <  1:2206472/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
TGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCA  >  NC_000913/3332725‑3332920

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: