Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,863,777 C→G 15.5% S62C (TCC→TGC)  ygbL → putative class II aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,863,7770CG15.5% 413.9 / 53.0 158S62C (TCC→TGC) ygbLputative class II aldolase
Reads supporting (aligned to +/- strand):  ref base C (75/57);  new base G (17/8);  total (93/65)
Fisher's exact test for biased strand distribution p-value = 3.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 171 total aligned reads displayed.

CCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGC  >  NC_000913/2863678‑2863876
                                                                                                   |                                                                                                   
ccGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGCCACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTc                                                                                                     <  2:2171532/100‑1 (MQ=255)
ccGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGCCAGGTTCATGCCTGGGCAAACACGATCCGCAGCGGCTTTc                                                                                                     >  1:998295/1‑100 (MQ=255)
ccGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTc                                                                                                     <  2:1895417/100‑1 (MQ=255)
ccGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTc                                                                                                     <  1:674466/100‑1 (MQ=255)
ccGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTc                                                                                                     <  1:2281886/100‑1 (MQ=255)
 cGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTcc                                                                                                    <  1:310767/100‑1 (MQ=255)
   gTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCaa                                                                                                  <  1:27996/100‑1 (MQ=255)
   gTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCAGGGCCATCTCGATCCGCAGCGGCTTTccca                                                                                                  >  2:2481637/1‑98 (MQ=255)
    ttCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCa                                                                                                   <  1:1816642/98‑1 (MQ=255)
     tCGGCTTGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAg                                                                                                <  1:945537/100‑1 (MQ=255)
       ggCTGGCAATCTGTCGCTGCTTTTACCTTACGGGAATTTACTGGCGACACCGACAGGAACATGCCTGAGCAATCTCGATCCGCAGCGGCTTCCCAAAg                                                                                                >  2:196472/1‑98 (MQ=255)
        gCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTcg                                                                                             <  1:1148011/100‑1 (MQ=255)
         cTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAg                                                                                                >  2:3017176/1‑96 (MQ=255)
          tGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTcgcc                                                                                           >  1:266381/1‑100 (MQ=255)
             cAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCcgcg                                                                                        >  1:2544070/1‑100 (MQ=255)
              aaTCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCGGCGACTTTGGAAAGCCGCTGCgg                                                                                       <  1:1012957/100‑1 (MQ=37)
                tCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGAt                                                                                     <  2:445639/100‑1 (MQ=255)
                tCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCACAGCGGCTTTCCAAAGTCGCCGCGGAt                                                                                     <  2:2933970/100‑1 (MQ=255)
                  tGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATgg                                                                                   <  1:3780774/100‑1 (MQ=255)
                  ngTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATgg                                                                                   >  2:3307489/2‑100 (MQ=255)
                   gTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACATCGATAGGTTCATGCCTGGGCAATCTCGATCCGCGGAGACTTTGGAAAGCCGCTGCGGAtcga                                                                                  <  1:2385718/100‑4 (MQ=21)
                   gTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGc                                                                                  <  2:1071663/100‑1 (MQ=255)
                    tCGCTGCTTTTACCTGACTGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCg                                                                                 <  1:1948089/100‑1 (MQ=255)
                    tCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCg                                                                                 <  2:2157572/100‑1 (MQ=255)
                      gctgctTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCGGCGACTTTGGAAAGCCGCTGCGGAtcgagat                                                                               <  1:3005227/100‑7 (MQ=21)
                      gctgctTTTACCTGAAGGGAATTTACTGGAGAAACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGATTTCAAAAGTCGCCGCTGATg                                                                                    >  2:1764950/1‑95 (MQ=255)
                       ctgctTTTACCTCAAGGGAATTTAATGTCGACAACGACAGGTTCATGCCTGGTGGATATCGGTTCACAGCGGCTGTCAAAAGgcgcgctgtatt                                                                                    >  2:1868845/1‑86 (MQ=255)
                              tACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAg                                                                       >  2:3131695/1‑100 (MQ=255)
                              tACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAg                                                                       >  1:2984282/1‑100 (MQ=255)
                                ccTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGTGCAATCTCGATCCGCAGAAGCTTTCCAAAGTCGAAGCGGATGGCGAATGGTTAAgt                                                                      >  1:94038/1‑99 (MQ=255)
                                ccTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAgtg                                                                     <  1:1215511/100‑1 (MQ=255)
                                ccTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAgtg                                                                     >  1:1079340/1‑100 (MQ=255)
                                 cTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAgtgn                                                                    <  1:3207091/100‑2 (MQ=255)
                                   gccGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGAACCGCAGCGGCGGTGCAAAACCCCCGCGGATGGCGAATTGTTAActggtg                                                                  >  2:452214/3‑100 (MQ=255)
                                   gACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAgtggtg                                                                  <  2:1451663/100‑1 (MQ=255)
                                   gACGGGAATTTACGGGCGACACCGACAGGATCAAGCCTGGGCAATCACGATCCGCAGCAGCTTTGAAAAAACGCCGCGGATGGCGAATGGTTAAgtggcg                                                                  >  2:2413634/1‑98 (MQ=255)
                                    aCGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATAGCGACTGTTTAAttggt                                                                   >  2:2599631/1‑98 (MQ=255)
                                     cGGGAATTTACTGGCGACACCGACAGGTTCACACCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAAACCCCCCGGATGGCGAACGCTTAAGTGGTGAc                                                                >  2:1238671/1‑100 (MQ=255)
                                       ggAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCTCCGCGGATGGCGAATGGTTAAGTGGTGACaa                                                              >  1:1880469/1‑100 (MQ=255)
                                       ggAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACaa                                                              >  1:3889581/1‑100 (MQ=255)
                                        gAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCACTCTGGATCCGCAGCGGCTTTGCAAAGTCGCAGCCGCTGGAGAGTggctgagggg                                                                    >  2:1132030/1‑85 (MQ=255)
                                        gAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAACGGTTTACCAATGACGCAGCGGCTGACGAATGGTTCATTGATGAcc                                                               >  2:3828860/1‑97 (MQ=255)
                                          aTTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAAcc                                                           <  2:41539/100‑1 (MQ=255)
                                             tACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTc                                                        >  1:2905924/1‑100 (MQ=255)
                                             tACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTc                                                        <  2:3627316/100‑1 (MQ=255)
                                               cTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGa                                                      >  1:3656623/1‑100 (MQ=255)
                                                 ggCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGaaa                                                    >  1:1808432/1‑100 (MQ=255)
                                                   cGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTGCAAAGTCGCCGCTGATGGCGAATTGTTAAGTGGTGACAAACCCTCGAAaga                                                  <  2:160090/100‑1 (MQ=255)
                                                   cGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAaga                                                  <  2:340808/100‑1 (MQ=255)
                                                     acacCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAgg                                                >  1:2172747/1‑100 (MQ=255)
                                                      cacCGACAGGTTCATGCCTGTGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGt                                               <  2:3861025/100‑1 (MQ=255)
                                                      cacCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGt                                               <  2:1899898/100‑1 (MQ=255)
                                                         cGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGc                                             >  2:2104702/1‑99 (MQ=255)
                                                           aCAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGctct                                          >  2:24404/1‑100 (MQ=255)
                                                            cAGGTTCATGCCTGGGCAATCTCGATTCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCtt                                         <  1:1239809/100‑1 (MQ=255)
                                                             aGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCttt                                        <  2:2536006/100‑1 (MQ=255)
                                                               gTTCATGCCTGGGCAATCTCGATCCTCAGCGGTTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCa                                      <  1:290804/100‑1 (MQ=255)
                                                               gTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCa                                      <  2:839559/100‑1 (MQ=255)
                                                               gTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCa                                      >  1:2101444/1‑100 (MQ=255)
                                                               gTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCa                                      >  1:794319/1‑100 (MQ=255)
                                                               gTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCa                                      <  1:3507985/100‑1 (MQ=255)
                                                                  cATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCt                                   <  2:738027/100‑1 (MQ=255)
                                                                  cATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATCGTTAAGTGGTCACAAACCCTCGAAAGAGTTGCTCTTTCAttt                                   <  1:2640134/100‑3 (MQ=255)
                                                                    tGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTgg                                 <  1:3131725/100‑1 (MQ=255)
                                                                    tGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAACTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTgg                                 <  2:2406364/100‑1 (MQ=255)
                                                                        tGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCt                             >  1:531761/1‑100 (MQ=255)
                                                                        tGGGCAATCTCGATCCGCAGCAGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGgc                                >  2:1345703/1‑97 (MQ=255)
                                                                            cAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACGAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTAt                         >  2:3465949/1‑100 (MQ=255)
                                                                               tctcGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGATGCTCTTTCATCTGGCGCTGTATCGc                      >  2:541601/1‑100 (MQ=255)
                                                                                ctcGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGca                     >  2:459729/1‑100 (MQ=255)
                                                                                ctcGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGca                     >  1:2283067/1‑100 (MQ=255)
                                                                                ctcGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGca                     <  2:2973281/100‑1 (MQ=255)
                                                                                cgccaTCCGCGGCGACTTTGGAAAGCCGCTGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGca                     >  2:1985750/5‑100 (MQ=25)
                                                                                cgccaTCCGCGGCGACTTTGGAAAGCCGCTGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTACTCTTTCATCTGGCGCTGTATCGca                     >  1:1664777/5‑100 (MQ=21)
                                                                                 tcGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGcaa                    >  2:1435437/1‑100 (MQ=255)
                                                                                 gccaTCCGCGGCGACTTTGGAAAGCCGCTGCGGATGGCGAATGGTTAAGTGGTGACAAACACTCGAAAGAGGAGCTCTTTCATCTGGCGCTGTATCGcaa                    >  1:1907717/4‑100 (MQ=21)
                                                                                  cGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGcaac                   <  1:2747250/100‑1 (MQ=255)
                                                                                  cGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGcaac                   <  2:2861572/100‑1 (MQ=255)
                                                                                   gATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGcaaca                  >  1:100871/1‑100 (MQ=255)
                                                                                   gATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGcaaca                  >  1:3246454/1‑100 (MQ=255)
                                                                                   gATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGcaaca                  >  1:3734418/1‑100 (MQ=255)
                                                                                   caTCCGCGGCGACTTTGGAAAGCCGCTGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGcaaca                  >  2:2155675/2‑100 (MQ=25)
                                                                                   caTCCGCGGCGACTTTGGAAAGCCGCTGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGcaaca                  >  1:3504506/2‑100 (MQ=37)
                                                                                    aTCCTCAGCGGCTTTCCAAAGTCGTCTCGTATTTCGAATGGTTAAGTTGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGcaac                   <  1:162181/98‑1 (MQ=255)
                                                                                    aTCCGCGGCGACTTTGGAAAGCCGCTGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGcaacaa                 >  2:264239/1‑100 (MQ=37)
                                                                                    aTCCGCGGCGACTTTGGAAAGCCGCTGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGcaacaa                 >  1:3628578/1‑100 (MQ=25)
                                                                                    aTCCGCGGCGACTTTGGAAAGCCGCTGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGcaacaa                 >  1:3065866/1‑100 (MQ=37)
                                                                                    aTCCGCGGCGACTTTGGAAAGCCGCTGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTAGCGCTGTATCGcaacaa                 >  1:3452076/1‑100 (MQ=25)
                                                                                      ccGCGGCGACTTTGGAAAGCCGCTGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAAAAAt                >  1:2752354/1‑99 (MQ=37)
                                                                                       cGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATcc              >  1:1169723/1‑100 (MQ=255)
                                                                                       cGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAAAAATcc              >  1:548021/1‑100 (MQ=255)
                                                                                            cGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCt         >  2:565279/1‑100 (MQ=255)
                                                                                             ggCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACGATCCGCGCt         <  2:379800/99‑1 (MQ=255)
                                                                                              gCTTTCCAAAGTCGTCGCGGATGGCGAATGGTAAAGTGGTGACAATCCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCt         <  2:2053238/98‑1 (MQ=255)
                                                                                              gCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGt       >  2:1960573/1‑100 (MQ=255)
                                                                                               cTTTCCACAGTAGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTACATCTGGCGCTGTATCGCAACAATCCGCGCTGTa      >  2:2047239/1‑100 (MQ=255)
                                                                                                  tCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGTCAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAg   <  1:2481660/100‑1 (MQ=255)
                                                                                                  tCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGAGCTATTTCATCTAGCGCTGTATCGCAACAATCCGCGCTGTAAAg   >  1:1284070/1‑100 (MQ=255)
                                                                                                   ccAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGc  >  2:2208455/1‑100 (MQ=255)
                                                                                                   ccAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGc  >  2:2805498/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CCGGTTCGGCTGGCAATCTGTCGCTGCTTTTACCTGACGGGAATTTACTGGCGACACCGACAGGTTCATGCCTGGGCAATCTCGATCCGCAGCGGCTTTCCAAAGTCGCCGCGGATGGCGAATGGTTAAGTGGTGACAAACCCTCGAAAGAGGTGCTCTTTCATCTGGCGCTGTATCGCAACAATCCGCGCTGTAAAGC  >  NC_000913/2863678‑2863876

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