Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,332,823 T→C 16.3% intergenic (‑88/+39) 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,8230TC16.3% 249.4 / 44.9 117intergenic (‑88/+39)yhbE/rpmAEamA family inner membrane putative transporter/50S ribosomal subunit protein L27
Reads supporting (aligned to +/- strand):  ref base T (57/40);  new base C (3/16);  total (60/57)
Fisher's exact test for biased strand distribution p-value = 7.80e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.29e-02

Only 100 of 120 total aligned reads displayed.

GCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCAC  >  NC_000913/3332726‑3332921
                                                                                                 |                                                                                                  
gctgctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACgtgt                                                                                                  <  2:1506723/100‑1 (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:313945/100‑1 (MQ=255)
       cATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCgggg                                                                                           <  2:2376307/100‑1 (MQ=255)
       cATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCgggg                                                                                           <  2:2500928/100‑1 (MQ=255)
       cATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCgggg                                                                                           <  1:198948/100‑1 (MQ=255)
       cATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCgggg                                                                                           <  2:811718/100‑1 (MQ=255)
        aTACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGc                                                                                          <  1:1657065/100‑1 (MQ=255)
        aTACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGc                                                                                          <  2:2025937/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)
         tACCGTATTCCCTGCAAAGCGCATTTGCGAATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGCt                                                                                         <  2:2653877/100‑1 (MQ=39)
          aCCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCtt                                                                                        >  1:2760486/1‑100 (MQ=255)
          aCCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGCtt                                                                                        <  2:2457995/100‑1 (MQ=255)
          aCCGTATGCCCTGCAAAGCGGATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCtt                                                                                        >  2:1766239/1‑100 (MQ=255)
           ccGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCttt                                                                                       >  1:3196/1‑100 (MQ=255)
           ccGTATTCCCTGCAAAGCGCATTTGAGGATTAACGTCCCACCATGCGCATATTTACCGGTTTTAGAATGTAAAAATCCCCTCAACGTGTTGCGGGGCttt                                                                                       >  2:1593965/1‑100 (MQ=255)
            cGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTc                                                                                      >  2:2108926/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)
               aTTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCCCCATGCGTATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATc                                                                                   <  1:891411/100‑1 (MQ=255)
                  gccTGCAACGCGCATTTGCAGATTAACGTTGCGCCGTGCGCATTTTTACCGGTTATCGAATGTAAAAAGCGCCGGAACGTGTTGCGGGTCTTTGATGCg                                                                                 >  2:2238685/2‑99 (MQ=255)
                   ccTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGtt                                                                               <  2:2644693/100‑1 (MQ=255)
                      gCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTAcc                                                                            >  1:1059265/1‑100 (MQ=255)
                      gCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTAcc                                                                            >  1:932414/1‑100 (MQ=255)
                        aaaGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCgg                                                                          <  2:1238458/100‑1 (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                                                                <  2:2320357/100‑1 (MQ=255)
                                  tGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGaaaa                                                                >  1:2705261/1‑100 (MQ=255)
                                   gCGGATTAACGTCGCACCATGCACATTTTTACCGGTTATCAAATGTAAAAAGCCCCGCAACATTTTGCGGGGATTTCATCCGTTACCGAGACGCGaaaaa                                                               >  2:771590/1‑100 (MQ=255)
                                    cGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAAc                                                              <  1:62645/100‑1 (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)
                                       aTTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTa                                                           >  2:363206/1‑100 (MQ=255)
                                        ttAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTAt                                                          <  1:2815166/100‑1 (MQ=255)
                                          aaCGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTc                                                        >  2:283997/1‑100 (MQ=255)
                                          aaCGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTc                                                        <  2:1020568/100‑1 (MQ=255)
                                            cGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAg                                                      >  1:635323/1‑100 (MQ=255)
                                             gTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGc                                                     <  2:1182625/100‑1 (MQ=255)
                                                gCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTc                                                  >  2:1379203/1‑100 (MQ=255)
                                                 cACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCg                                                 <  1:1087015/100‑1 (MQ=255)
                                                  aCCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCACAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGa                                                <  2:2834442/100‑1 (MQ=255)
                                                       gcgcATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGAGACGCGAAAAAATTATTCAGCTTCGATGCTg                                           >  1:308532/1‑100 (MQ=255)
                                                        cgcATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGa                                          >  1:2002319/1‑100 (MQ=255)
                                                        cgcATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGa                                          >  1:752510/1‑100 (MQ=255)
                                                          cATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGAtn                                        <  1:2369835/100‑2 (MQ=255)
                                                          cATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATa                                        >  2:679225/1‑100 (MQ=255)
                                                           aTTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATaa                                       >  2:2543139/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:1389853/1‑100 (MQ=255)
                                                                    ggTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGt                              >  2:1221232/1‑100 (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)
                                                                            acggatgaaaGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCggg                      <  2:2850114/93‑1 (MQ=255)
                                                                             aTGTAAAAAGCCCCGCAACGTGTTGCGGTGCTTTCATCCGTTACCGTGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCggg                      <  1:1484321/99‑1 (MQ=255)
                                                                                  aaaattcccGCAACGGGTTGCTTGGCGTTCAGCTGTTACCGGGACTCGAATGACTTATTCTGCTTTGATGCTCATAAATTTACGGTTTTTCGGGGCttt                 <  2:555257/93‑1 (MQ=255)
                                                                                  aaaaGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTa                >  1:2797074/1‑100 (MQ=255)
                                                                                    aaGCCCCGCAACGTGTTGCGAGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAAc              >  2:126075/1‑100 (MQ=255)
                                                                                     aGCCCCGCAACGTGTTGCGCGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGAAGCTGATAAATTTAAGAAATTTAGAGCCTTAAACt             >  1:905654/1‑100 (MQ=255)
                                                                                     aGCCCCGCAACACGTTGCGGGGCTTTCATCCGTTACCGGGACGCTAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACt             <  1:801473/100‑1 (MQ=39)
                                                                                      gCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACtt            >  1:170092/1‑100 (MQ=255)
                                                                                      gCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACtt            >  1:2460053/1‑100 (MQ=255)
                                                                                      gCCCCGCAACACGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACt             >  1:56020/1‑99 (MQ=255)
                                                                                       ccccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTc           >  1:2243262/1‑100 (MQ=255)
                                                                                        cccGCTACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGTCCTTTCACTTCg          >  2:1493906/1‑100 (MQ=255)
                                                                                        cccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGTGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCg          >  2:217607/1‑100 (MQ=255)
                                                                                        cccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCg          >  2:767918/1‑100 (MQ=255)
                                                                                         ccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGa         <  2:2386299/100‑1 (MQ=255)
                                                                                         ccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGa         >  2:2380707/1‑100 (MQ=255)
                                                                                          cGCAACGTGTTGCGGGTCTTTCATCCGTTACCGTGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGaa        <  1:2816761/100‑1 (MQ=255)
                                                                                              acttgtTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTc    <  2:298768/97‑1 (MQ=255)
                                                                                              aCGTGTTGCGGGGATTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTc    >  2:793430/1‑100 (MQ=255)
                                                                                               cGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAAttgc    <  2:1554021/99‑3 (MQ=255)
                                                                                               cGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCa   <  1:2206472/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
GCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCAC  >  NC_000913/3332726‑3332921

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: