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)
                  ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTg                                                                                 >  2:3461212/1‑100 (MQ=255)
                  ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTg                                                                                 >  2:3142401/1‑100 (MQ=255)
                   gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGt                                                                                <  2:3296361/100‑1 (MQ=255)
                     aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCg                                                                              <  2:898084/100‑1 (MQ=255)
                      aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGt                                                                             <  1:1058065/100‑1 (MQ=255)
                        cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGa                                                                           <  1:369105/100‑1 (MQ=255)
                        cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGa                                                                           >  1:545495/1‑100 (MQ=255)
                            cgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTATTCATAAGAAACCTACCAGCGGTGGCACGttttttttgacggtc                                                                                >  1:3755111/1‑85 (MQ=255)
                            cgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAgcgc                                                                       <  1:108956/100‑1 (MQ=255)
                                gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGcaca                                                                    <  1:2357205/99‑1 (MQ=255)
                                gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGCCGTGAGCGcactc                                                                   <  2:2581579/100‑3 (MQ=255)
                                 ccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACAcc                                                                  <  1:638162/100‑1 (MQ=255)
                                     gcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATg                                                              >  2:3688041/1‑100 (MQ=255)
                                     gcCAACCGCGCACCGAAAGACGATTAATTGCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGTCAGG‑TTTGTTTTGCCTGTCGTGAGCGCACACCATTg                                                              <  1:3218264/100‑1 (MQ=255)
                                       cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCg                                                            >  2:3503853/1‑100 (MQ=255)
                                        aaCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                           <  2:1947075/100‑1 (MQ=255)
                                         aCCGCGCACCGAAAGACGATTAATTTCGGTTTTTAGTCATAAAAAACCTGCCACCGCTGGCAGGTTTTTTTTT‑CCTGTCGTGAGCGCACACCAATGCGa                                                           >  1:3153330/1‑99 (MQ=38)
                                           cgcgCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTa                                                        <  1:2664895/100‑1 (MQ=255)
                                              gcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGGAtct                                                     >  2:3618613/1‑100 (MQ=255)
                                               cACCGAAAGACGATTAATTTCGGTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTTGTGAGCGCACACCAATGCGAGTAtctc                                                    <  2:3036015/100‑1 (MQ=255)
                                                aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTcc                                                   <  2:1477698/100‑1 (MQ=255)
                                                 ccGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                                  <  2:2929053/100‑1 (MQ=255)
                                                 ccGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                                  >  2:534298/1‑100 (MQ=255)
                                                  cGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCt                                                 <  2:3649360/100‑1 (MQ=255)
                                                   gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCtt                                                <  2:1907225/100‑1 (MQ=255)
                                                   gAAAGACGATCAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTTTCCTGAGCGCACACCAATGCGAGTATCTCCCtt                                                <  1:2492779/100‑1 (MQ=255)
                                                     aaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtat                                              >  2:1793207/1‑100 (MQ=255)
                                                         cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGt                                          >  1:963189/1‑100 (MQ=255)
                                                          gATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTc                                         >  2:226040/1‑100 (MQ=255)
                                                           aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGGAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCt                                        <  1:2664998/100‑1 (MQ=255)
                                                           aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGCGTATCTCCCTTATAAGTCt                                        >  1:2414456/1‑100 (MQ=255)
                                                              aaTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                     <  2:311238/100‑1 (MQ=255)
                                                               aTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                    <  2:1919646/100‑1 (MQ=255)
                                                                tggcGTTTGTAAGTAATAAACGACCGGCGAGCGATGGCAGG‑GTTTTTTTGCCTGTCGTGAGCGCACCCCAATGCGCGTATCTCCCTTATAAGTCTGTgg                                    <  1:3276358/96‑1 (MQ=255)
                                                                 ttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGTGAGTATCTTCCTTATAAGTCTGTGttt                                  >  1:2359799/1‑100 (MQ=255)
                                                                 ttCATTTTTTAGTCATAAAAACCCTGCCAGCGATCCCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATATGTCTGTgg                                    <  1:839190/98‑1 (MQ=255)
                                                                          tAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                         <  1:3510493/100‑1 (MQ=255)
                                                                           aGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                        >  1:3676477/1‑100 (MQ=255)
                                                                            gTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                       >  2:3314854/1‑100 (MQ=255)
                                                                             taaTAAAAAACCTCCCAGCGATGGCAGG‑TTTATTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                      >  2:2565147/3‑100 (MQ=255)
                                                                             tCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                         <  2:2420544/97‑1 (MQ=255)
                                                                             gcaaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                      >  2:2911174/5‑100 (MQ=38)
                                                                              caaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                     >  2:376436/4‑100 (MQ=38)
                                                                              caaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                     >  2:449020/4‑100 (MQ=38)
                                                                              caaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                     >  1:988077/4‑100 (MQ=38)
                                                                              cATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCATTATAAGTCTGCGGTTACGTCAATagaga                     >  1:929086/1‑100 (MQ=255)
                                                                               aTAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                    <  1:1926350/100‑1 (MQ=255)
                                                                               aTAAAAAACCTGCCAGCGATGGCA‑GTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                     <  1:893920/99‑1 (MQ=255)
                                                                                aaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                   >  1:2031204/2‑100 (MQ=38)
                                                                                aaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTATCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                   >  1:2023127/2‑100 (MQ=38)
                                                                                aaaaaaaCCTGCCATCGCGGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATGAGTCTGTGGTTACGTCAATAGAGAGc                   >  2:1384051/2‑100 (MQ=38)
                                                                                  aaaaaCCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                   >  1:600818/1‑98 (MQ=255)
                                                                                   aaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt                >  2:2479143/1‑100 (MQ=39)
                                                                                    aaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg               >  2:1701890/1‑100 (MQ=39)
                                                                                      aCCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc             <  2:2961005/100‑1 (MQ=255)
                                                                                      aCCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc             <  1:483109/100‑1 (MQ=255)
                                                                                       ccTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTca            >  2:1230430/1‑100 (MQ=255)
                                                                                         tGCCAGCGATGGCAGG‑TTTTTTTTGTTTTTTTTTAGCGCACACCAATGCGAGTATTTCCTTTATAAGTCTGTGGTTACGTCTATAGAGCGCTTTATcaca          >  1:3391576/1‑100 (MQ=255)
                                                                                          gCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa         >  2:3612806/1‑100 (MQ=255)
                                                                                           ccAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt        <  2:1730968/100‑1 (MQ=255)
                                                                                               cGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATc    <  1:2825446/100‑1 (MQ=255)
                                                                                                 aTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTg  <  2:600752/100‑1 (MQ=255)
                                                                                                 aTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTg  <  1:1869950/100‑1 (MQ=255)
                                                                                                 aTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTg  <  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: ‑