Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,060,280 (T)9→8 100% intergenic (+50/‑167) 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,2720T.94.2% 417.7 / 18.5 103intergenic (+42/‑175)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (2/4);  new base . (45/52);  total (47/56)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

Only 100 of 136 total aligned reads displayed.

CGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAA  >  NC_000913/4060174‑4060367
                                                                                                  |                                                                                               
cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtt                                                                                                >  1:292442/1‑100 (MQ=255)
  tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGGCAAAAAAAACCTGCCATCGCTGGCAGGtttt                                                                                              <  2:217112/100‑1 (MQ=38)
   cACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttt                                                                                             <  2:955304/100‑1 (MQ=255)
   cACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttt                                                                                             <  2:2942476/100‑1 (MQ=255)
   cACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttt                                                                                             <  2:2706428/100‑1 (MQ=255)
   cACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAAAAAAAACCTGCCATCGCTGGCAGGttttt                                                                                             <  1:1918620/100‑1 (MQ=38)
   cACCTGACCGAAAGCGATCGTTGTCACGTGACCCTCTAACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttt                                                                                              <  1:2672840/99‑1 (MQ=19)
    aCCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttttt                                                                                            <  1:1090299/100‑1 (MQ=255)
    aCCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttt                                                                                            <  1:1645533/100‑1 (MQ=255)
     ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                           <  1:1537523/100‑1 (MQ=255)
     ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                           <  1:2105707/100‑1 (MQ=255)
      cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAGTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttttt                                                                                          <  2:1111477/100‑1 (MQ=255)
      cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttat                                                                                          <  2:1391937/100‑3 (MQ=255)
      cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttttt                                                                                          >  2:1575007/1‑100 (MQ=255)
      cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttttt                                                                                          <  2:324421/100‑1 (MQ=255)
       tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatg                                                                                         <  2:2424329/100‑4 (MQ=38)
       tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttg                                                                                         <  1:92309/100‑2 (MQ=255)
       tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttttgc                                                                                         <  1:971633/100‑3 (MQ=255)
        gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatga                                                                                        <  2:2848015/100‑5 (MQ=38)
        gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACAGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgc                                                                                        <  1:2438650/100‑3 (MQ=255)
         aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcc                                                                                       <  1:1400582/100‑4 (MQ=255)
         aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcc                                                                                       <  2:2281181/100‑4 (MQ=255)
         aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcc                                                                                       <  1:71874/100‑4 (MQ=255)
            gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATCAAAAAGCTGCCAGCTATGGCAGG‑TTTTTTTTGCCTGt                                                                                   >  1:2058264/1‑100 (MQ=255)
            gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGt                                                                                   >  2:2035068/1‑100 (MQ=255)
                aCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTg                                                                               <  2:1267562/100‑1 (MQ=255)
                aCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTg                                                                               <  1:6170/100‑1 (MQ=255)
                 cGATCGTCGCCGCGCCAACCGCTCACCGGAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTc                                                                                  <  2:2769755/96‑1 (MQ=255)
                 cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGt                                                                                >  1:707224/1‑99 (MQ=255)
                 cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGa                                                                              >  1:1142186/1‑100 (MQ=255)
                   aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAgc                                                                            >  1:2032587/1‑100 (MQ=255)
                   aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTGTTGCCTGTCGTGAgc                                                                            >  2:3014821/1‑100 (MQ=255)
                       tcgcgccgcCAACCGCGCAGCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGca                                                                         <  1:2427732/93‑1 (MQ=255)
                       tcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGcac                                                                        >  1:2505570/1‑100 (MQ=255)
                       tcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTGTTTTTGCCTGTCGTGAGCGcac                                                                        >  1:882212/1‑100 (MQ=255)
                        cgccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAATACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGcaca                                                                       <  1:545365/100‑1 (MQ=255)
                        cgccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGcaca                                                                       <  2:2320548/100‑1 (MQ=255)
                         gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTTTTTTTTTTTTACGTGCGAGAGGTGGCAGG‑TTTTTTTTGCCTGTCGTGAGGGCacac                                                                      >  1:2677300/1‑100 (MQ=255)
                         gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGTCAGG‑TTTTTTTTGCCTGTCGTGAGCGcactc                                                                      <  2:2058276/100‑3 (MQ=255)
                         gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCacac                                                                      >  1:2302288/1‑100 (MQ=255)
                         gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCacac                                                                      <  2:1206929/100‑1 (MQ=255)
                           cgcaCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCa                                                                    >  2:649119/1‑100 (MQ=255)
                                cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCg                                                               >  2:1026276/1‑100 (MQ=255)
                                 aaCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                              <  2:1711434/100‑1 (MQ=255)
                                  aCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                             <  1:414534/100‑1 (MQ=255)
                                  aCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                             <  2:2751595/100‑1 (MQ=255)
                                   ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                             <  1:1795268/100‑1 (MQ=255)
                                        cACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtctc                                                       <  1:1675091/100‑1 (MQ=255)
                                        cACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtctc                                                       >  2:1910847/1‑100 (MQ=255)
                                          ccGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTcc                                                      >  2:2641354/1‑100 (MQ=255)
                                          ccGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                                     >  2:1637880/1‑100 (MQ=255)
                                           cGAAAGACTATTAATTTCGTTTTTTCTTCATAAAAAGCGTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTcc                                                      <  1:2582275/98‑1 (MQ=255)
                                           cGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCTAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCt                                                    <  1:2140745/100‑1 (MQ=255)
                                            gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCtt                                                   <  1:977484/100‑1 (MQ=255)
                                               aGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtata                                                <  2:1657512/100‑1 (MQ=255)
                                               aGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtata                                                <  2:428514/100‑1 (MQ=255)
                                                gACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATaa                                               <  1:155949/100‑1 (MQ=255)
                                                gACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATaa                                               <  1:3025491/100‑1 (MQ=255)
                                                gACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATaa                                               <  1:314057/100‑1 (MQ=255)
                                                 aCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAg                                              <  1:1818509/100‑1 (MQ=255)
                                                 aCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAg                                              <  1:1342724/100‑1 (MQ=255)
                                                 aCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGACGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAg                                              <  1:2152035/100‑1 (MQ=255)
                                                  cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAg                                              <  1:1896135/99‑1 (MQ=255)
                                                  cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGt                                             >  2:3057172/1‑100 (MQ=255)
                                                      tAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgt                                         <  1:2507904/100‑1 (MQ=255)
                                                       aaTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                        >  2:1741979/1‑100 (MQ=255)
                                                       aaTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                        <  2:1134391/100‑1 (MQ=255)
                                                        aTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                       <  1:2281050/100‑1 (MQ=255)
                                                            cGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                                   <  2:460293/100‑1 (MQ=255)
                                                             gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                                  >  2:386508/1‑100 (MQ=255)
                                                               tttttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGATCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                                <  2:1766695/100‑1 (MQ=255)
                                                                ttttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCa                               >  2:2018014/1‑100 (MQ=255)
                                                                 nttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCaa                              >  1:2001085/2‑100 (MQ=255)
                                                                  ttAGTCATAAAAAACCTGTCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAAt                             <  2:494450/100‑1 (MQ=255)
                                                                   tAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                            <  1:2490721/100‑1 (MQ=255)
                                                                   tAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACAACAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                            <  2:2908485/100‑1 (MQ=255)
                                                                    aGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                            <  2:707229/100‑1 (MQ=255)
                                                                    aGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                           >  2:2922346/1‑100 (MQ=255)
                                                                    aGTCATAAAAAACCTGCCAGCGATGCCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                           >  1:2252445/1‑100 (MQ=255)
                                                                     ggcaaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                          >  2:1569019/6‑100 (MQ=37)
                                                                       caaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                        >  2:1805958/4‑100 (MQ=38)
                                                                       caaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                        >  2:2162572/4‑100 (MQ=38)
                                                                       cATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                        <  2:2142379/100‑1 (MQ=255)
                                                                        aTAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                       >  1:1449003/1‑100 (MQ=255)
                                                                         tAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                      >  2:2542998/1‑100 (MQ=255)
                                                                            aaaaCCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt                   <  2:895165/100‑1 (MQ=255)
                                                                             aaaCCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGGGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagatcggaaga                  >  2:2406464/1‑92 (MQ=255)
                                                                               aCCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc                <  2:2617122/100‑1 (MQ=255)
                                                                               aCCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc                <  1:1297540/100‑1 (MQ=255)
                                                                               aCCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc                <  1:735077/100‑1 (MQ=255)
                                                                                ccTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc                <  2:1561367/99‑1 (MQ=255)
                                                                                  tGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcac              <  1:379463/100‑1 (MQ=255)
                                                                                  tGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcaca             <  2:1978923/100‑1 (MQ=255)
                                                                                  tGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcaca             <  2:136576/100‑1 (MQ=255)
                                                                                   gCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa            <  2:1189789/100‑1 (MQ=255)
                                                                                    ccAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt           >  1:403391/1‑100 (MQ=255)
                                                                                     cAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt           <  2:2613962/99‑1 (MQ=255)
                                                                                          aTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTg     <  1:885797/100‑1 (MQ=255)
                                                                                            ggCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCa   >  1:853415/1‑100 (MQ=255)
                                                                                             gCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCaa  <  2:703614/100‑1 (MQ=255)
                                                                                                  |                                                                                               
CGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCAA  >  NC_000913/4060174‑4060367

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