Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,060,264 A→C 15.9% intergenic (+34/‑183) typA → / → yihL GTP‑binding protein/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,060,2640AC15.9% 387.2 / 61.8 154intergenic (+34/‑183)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (50/78);  new base C (17/8);  total (68/86)
Fisher's exact test for biased strand distribution p-value = 1.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 157 total aligned reads displayed.

TCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTG  >  NC_000913/4060167‑4060364
                                                                                                 |                                                                                                    
tCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATg                                                                                                    <  2:2647825/100‑1 (MQ=255)
 cGTAAACGTCTCCGGACTGAAAACGATCGTCGCCGCGTCAACCGCGAACAGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATgg                                                                                                   <  1:939713/100‑1 (MQ=255)
 cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATgg                                                                                                   <  1:296206/100‑1 (MQ=255)
  gTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                  >  1:2450290/1‑100 (MQ=255)
  gTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCAAGCGATTGc                                                                                                  >  1:1694109/1‑100 (MQ=255)
   tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCa                                                                                                 <  2:788715/100‑1 (MQ=255)
    aaaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAg                                                                                                <  1:3663582/100‑1 (MQ=255)
     aaCGTCACCTGACGGAAAACGATCGTCGCCGCGGCAACCGCGCAACGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAgg                                                                                               <  1:2874791/100‑1 (MQ=255)
     aaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAgg                                                                                               >  2:638166/1‑100 (MQ=255)
     aaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAgg                                                                                               <  2:2515017/100‑1 (MQ=255)
      aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGt                                                                                              <  2:2480159/100‑1 (MQ=255)
       cGTCACCTAACGAAAAACGATCGTCACCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGt                                                                                              <  1:3470902/99‑1 (MQ=255)
        gTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttt                                                                                            <  2:81342/100‑1 (MQ=255)
        gTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTCAGTCATAAAAAACCTGCCAGCGATGGCAGGttt                                                                                            <  1:2275437/100‑1 (MQ=255)
         tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttt                                                                                           >  2:1659146/1‑100 (MQ=255)
          cACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttt                                                                                          >  2:116685/1‑100 (MQ=255)
           aCCTGACGGAAAACGATCGTCTCCGCTCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAATACCTGCCATCGCTGGCAGGtttttt                                                                                         <  1:2549523/100‑1 (MQ=38)
           aCCTGACGGAAAACGATCGTCGCCGCGTCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttt                                                                                          <  1:724173/99‑1 (MQ=255)
           aCCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttt                                                                                         >  2:3050261/1‑100 (MQ=255)
           aCCTGACAGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttt                                                                                         >  2:1995250/1‑100 (MQ=255)
            ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGGTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                        >  1:2979854/1‑100 (MQ=255)
            ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttttta                                                                                        <  2:2633037/100‑2 (MQ=255)
             cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttat                                                                                       <  1:2521170/100‑3 (MQ=255)
             cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttttt                                                                                       <  2:3385522/100‑1 (MQ=255)
             cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttttt                                                                                       >  1:224608/1‑100 (MQ=255)
             cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttttt                                                                                       >  2:1592002/1‑100 (MQ=255)
              tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttg                                                                                      <  2:1406745/100‑2 (MQ=255)
               gACGGAAAACGATTGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgc                                                                                     <  1:1231200/100‑3 (MQ=255)
               gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgc                                                                                     <  1:210411/100‑3 (MQ=255)
                aCGGAAACCGATCGTTGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcc                                                                                    <  1:1623966/100‑4 (MQ=255)
                aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcc                                                                                    <  2:2026149/100‑4 (MQ=255)
                aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcc                                                                                    <  1:1240980/100‑4 (MQ=255)
                aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcc                                                                                    <  1:1429259/100‑4 (MQ=255)
                aCGGAAAACGATCGTCGCCGCGCCAACCGCCCACCTAAAGACGTTTAATTTCGTTTTTTAGTCATAACCAACCTGCCCGCGATGGCAGGttttttttgcc                                                                                    <  1:1428271/100‑4 (MQ=255)
                 cGGAAAACGATCGTCGTCGCGCAAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTTATAAAAAACCTGCCAGCGATGGCAGGtgttttttgact                                                                                   >  1:439016/1‑96 (MQ=255)
                 cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcct                                                                                   >  2:1806930/1‑96 (MQ=255)
                  ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTg                                                                                 >  2:3461212/1‑100 (MQ=255)
                  ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTg                                                                                 >  2:3142401/1‑100 (MQ=255)
                   gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGt                                                                                <  2:3296361/100‑1 (MQ=255)
                     aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCg                                                                              <  2:898084/100‑1 (MQ=255)
                      aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGt                                                                             <  1:1058065/100‑1 (MQ=255)
                        cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGa                                                                           <  1:369105/100‑1 (MQ=255)
                        cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGa                                                                           >  1:545495/1‑100 (MQ=255)
                            cgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTATTCATAAGAAACCTACCAGCGGTGGCACGttttttttgacggtc                                                                                >  1:3755111/1‑85 (MQ=255)
                            cgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAgcgc                                                                       <  1:108956/100‑1 (MQ=255)
                                gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGcaca                                                                    <  1:2357205/99‑1 (MQ=255)
                                gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGCCGTGAGCGcactc                                                                   <  2:2581579/100‑3 (MQ=255)
                                 ccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                                  <  1:638162/100‑1 (MQ=255)
                                     gcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATg                                                              >  2:3688041/1‑100 (MQ=255)
                                     gcCAACCGCGCACCGAAAGACGATTAATTGCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGTCAGGTTTGTTTTGCCTGTCGTGAGCGCACACCATTg                                                              <  1:3218264/100‑1 (MQ=255)
                                       cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCg                                                            >  2:3503853/1‑100 (MQ=255)
                                        aaCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                           <  2:1947075/100‑1 (MQ=255)
                                         aCCGCGCACCGAAAGACGATTAATTTCGGTTTTTAGTCATAAAAAACCTGCCACCGCTGGCAGGTTTTTTTTTCCTGTCGTGAGCGCACACCAATGCGa                                                           >  1:3153330/1‑99 (MQ=38)
                                           cgcgCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTa                                                        <  1:2664895/100‑1 (MQ=255)
                                              gcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGGAtct                                                     >  2:3618613/1‑100 (MQ=255)
                                               cACCGAAAGACGATTAATTTCGGTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTTGTGAGCGCACACCAATGCGAGTAtctc                                                    <  2:3036015/100‑1 (MQ=255)
                                                aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTcc                                                   <  2:1477698/100‑1 (MQ=255)
                                                 ccGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                                  <  2:2929053/100‑1 (MQ=255)
                                                 ccGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                                  >  2:534298/1‑100 (MQ=255)
                                                  cGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCt                                                 <  2:3649360/100‑1 (MQ=255)
                                                   gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCtt                                                <  2:1907225/100‑1 (MQ=255)
                                                   gAAAGACGATCAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTTTCCTGAGCGCACACCAATGCGAGTATCTCCCtt                                                <  1:2492779/100‑1 (MQ=255)
                                                     aaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtat                                              >  2:1793207/1‑100 (MQ=255)
                                                         cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGt                                          >  1:963189/1‑100 (MQ=255)
                                                          gATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTc                                         >  2:226040/1‑100 (MQ=255)
                                                           aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGGAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCt                                        <  1:2664998/100‑1 (MQ=255)
                                                           aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGCGTATCTCCCTTATAAGTCt                                        >  1:2414456/1‑100 (MQ=255)
                                                              aaTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                     <  2:311238/100‑1 (MQ=255)
                                                               aTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                    <  2:1919646/100‑1 (MQ=255)
                                                                tggcGTTTGTAAGTAATAAACGACCGGCGAGCGATGGCAGGGTTTTTTTGCCTGTCGTGAGCGCACCCCAATGCGCGTATCTCCCTTATAAGTCTGTgg                                    <  1:3276358/96‑1 (MQ=255)
                                                                 ttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGTGAGTATCTTCCTTATAAGTCTGTGttt                                  >  1:2359799/1‑100 (MQ=255)
                                                                 ttCATTTTTTAGTCATAAAAACCCTGCCAGCGATCCCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATATGTCTGTgg                                    <  1:839190/98‑1 (MQ=255)
                                                                          tAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                         <  1:3510493/100‑1 (MQ=255)
                                                                           aGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                        >  1:3676477/1‑100 (MQ=255)
                                                                            gTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                       >  2:3314854/1‑100 (MQ=255)
                                                                             taaTAAAAAACCTCCCAGCGATGGCAGGTTTATTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                      >  2:2565147/3‑100 (MQ=255)
                                                                             tCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                         <  2:2420544/97‑1 (MQ=255)
                                                                             gcaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                      >  2:2911174/5‑100 (MQ=38)
                                                                              caaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                     >  2:376436/4‑100 (MQ=38)
                                                                              caaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                     >  2:449020/4‑100 (MQ=38)
                                                                              caaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                     >  1:988077/4‑100 (MQ=38)
                                                                              cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCATTATAAGTCTGCGGTTACGTCAATagaga                     >  1:929086/1‑100 (MQ=255)
                                                                               aTAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                    <  1:1926350/100‑1 (MQ=255)
                                                                               aTAAAAAACCTGCCAGCGATGGCAGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                     <  1:893920/99‑1 (MQ=255)
                                                                                aaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                   >  1:2031204/2‑100 (MQ=38)
                                                                                aaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTATCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                   >  1:2023127/2‑100 (MQ=38)
                                                                                aaaaaaaCCTGCCATCGCGGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATGAGTCTGTGGTTACGTCAATAGAGAGc                   >  2:1384051/2‑100 (MQ=38)
                                                                                  aaaaaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                   >  1:600818/1‑98 (MQ=255)
                                                                                   aaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt                >  2:2479143/1‑100 (MQ=39)
                                                                                    aaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg               >  2:1701890/1‑100 (MQ=39)
                                                                                      aCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc             <  2:2961005/100‑1 (MQ=255)
                                                                                      aCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc             <  1:483109/100‑1 (MQ=255)
                                                                                       ccTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTca            >  2:1230430/1‑100 (MQ=255)
                                                                                         tGCCAGCGATGGCAGGTTTTTTTTGTTTTTTTTTAGCGCACACCAATGCGAGTATTTCCTTTATAAGTCTGTGGTTACGTCTATAGAGCGCTTTATcaca          >  1:3391576/1‑100 (MQ=255)
                                                                                          gCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa         >  2:3612806/1‑100 (MQ=255)
                                                                                           ccAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt        <  2:1730968/100‑1 (MQ=255)
                                                                                               cGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATc    <  1:2825446/100‑1 (MQ=255)
                                                                                                 aTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTg  <  2:600752/100‑1 (MQ=255)
                                                                                                 aTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTg  <  1:1869950/100‑1 (MQ=255)
                                                                                                 aTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTg  <  1:1736344/100‑1 (MQ=255)
                                                                                                 |                                                                                                    
TCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTG  >  NC_000913/4060167‑4060364

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: