Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,060,261 G→T 15.7% intergenic (+31/‑186) 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,2610GT15.7% 466.1 / 66.1 186intergenic (+31/‑186)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (81/74);  new base T (15/15);  total (97/89)
Fisher's exact test for biased strand distribution p-value = 8.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-01

Only 100 of 191 total aligned reads displayed.

CGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTAT  >  NC_000913/4060162‑4060361
                                                                                                   |                                                                                                    
cGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAg                                                                                                      >  1:3760513/1‑100 (MQ=255)
cGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAg                                                                                                      >  2:2598617/1‑100 (MQ=255)
 gTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGc                                                                                                     >  1:661319/1‑100 (MQ=255)
     tCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATg                                                                                                 >  1:1551505/1‑100 (MQ=255)
      cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGTGCTTTCTGTGCACTGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATgg                                                                                                <  1:436286/100‑1 (MQ=255)
      cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATgg                                                                                                <  1:2570116/100‑1 (MQ=255)
      cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATgg                                                                                                >  1:1893350/1‑100 (MQ=255)
        tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCa                                                                                              >  1:3639739/1‑100 (MQ=255)
            cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtt                                                                                          >  1:1162137/1‑100 (MQ=255)
             gTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtt                                                                                          <  1:1852899/99‑1 (MQ=255)
             gTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttt                                                                                         >  1:531786/1‑100 (MQ=255)
               cACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttt                                                                                       <  1:1994523/100‑1 (MQ=255)
                aCCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttt                                                                                        >  2:56559/1‑98 (MQ=255)
                 ccTGACGGAAAACGTTCGTTGTCGCGCTTACCGTGCAACGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttttt                                                                                      <  1:171345/99‑1 (MQ=25)
                 ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                     <  1:3284928/100‑1 (MQ=255)
                 ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                     <  1:3452143/100‑1 (MQ=255)
                 ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                     <  2:2207932/100‑1 (MQ=255)
                  cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGTACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttttta                                                                                     <  1:344061/99‑2 (MQ=38)
                  cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttttt                                                                                    >  1:2827237/1‑100 (MQ=255)
                   ttgcGGAAATCGATCGTCGTCGGGCTGACCTCGCACTGGGAGGCGGTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGGGATGGCAGGtttttttt                                                                                    <  2:3640851/96‑1 (MQ=255)
                   tGACGGAAAACGATCGTCGCCGTGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatg                                                                                   <  2:140708/100‑4 (MQ=38)
                   tGACGGAAAACGATCGTCGCCGCGCCAACCGCTCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatg                                                                                   <  1:180461/100‑4 (MQ=38)
                   tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatg                                                                                   <  1:208462/100‑4 (MQ=255)
                   tGACGGAAAACGATCGTCGCCGCGCCAACCACGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttg                                                                                   <  2:2646438/100‑2 (MQ=255)
                    gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgc                                                                                  >  1:2212335/1‑98 (MQ=255)
                    gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttg                                                                                   >  1:3883889/1‑98 (MQ=255)
                     aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcc                                                                                 <  2:674419/100‑4 (MQ=255)
                     aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcc                                                                                 <  2:2776413/100‑4 (MQ=255)
                      cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatgact                                                                                <  2:3768296/100‑7 (MQ=37)
                      cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcat                                                                                >  1:1638558/1‑96 (MQ=255)
                       ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTg                                                                              >  1:3379183/1‑100 (MQ=255)
                       ggAAAACGAGCGGCGCCGGTCCAGCCGCGTAACGAAAGCCGATTAATTTCGTGTGTTAGTCATAAAATACCTGCCAGCGATGGCAGGttttttt                                                                                     <  2:3236732/94‑1 (MQ=255)
                        gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGt                                                                             >  1:1733493/1‑100 (MQ=255)
                              gATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAg                                                                       <  2:1699289/100‑1 (MQ=255)
                               cgatctgccgcGCCAACCGCGCACCGAAAGTCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAgc                                                                      <  1:2995110/94‑1 (MQ=255)
                                tcgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAgcg                                                                     <  2:1093363/100‑1 (MQ=255)
                                tcgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCATGTCGTGAgcg                                                                     <  2:1638543/100‑1 (MQ=255)
                                tcgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCGGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAg                                                                       <  1:2753996/99‑1 (MQ=255)
                                 cgtcgCCGCGCCAACCGCGCACTGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAgcgc                                                                    <  1:1485521/100‑1 (MQ=255)
                                  gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTAAAAGCAATGGCAGG‑AATTTTTTGCCTGTCGTGAGCGca                                                                   >  1:2196104/1‑100 (MQ=255)
                                   tcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGcac                                                                  <  1:2257241/100‑1 (MQ=255)
                                    cgccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGcaca                                                                 <  1:1103681/100‑1 (MQ=255)
                                    cgccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGcaca                                                                 >  2:47332/1‑100 (MQ=255)
                                     gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAgcagatcg                                                                >  2:2995083/1‑94 (MQ=255)
                                     gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCacac                                                                >  2:829087/1‑100 (MQ=255)
                                      ccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACAcc                                                               <  2:2918193/100‑1 (MQ=255)
                                        gcgcCAACCGCGCAACGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCaa                                                             >  2:927589/1‑100 (MQ=255)
                                           ccATCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATg                                                           <  1:595077/99‑1 (MQ=255)
                                           ccAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCATGTAGTGAGCGCACACCAATGc                                                          >  1:3114709/1‑100 (MQ=255)
                                            cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGc                                                          >  2:2024510/1‑99 (MQ=255)
                                            cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCg                                                         >  2:58965/1‑100 (MQ=255)
                                            cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGACTGTCGTGAGCGCACACCAATGCg                                                         >  2:2138487/1‑100 (MQ=255)
                                             aaCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                        >  2:1824498/1‑100 (MQ=255)
                                             aaCCGCGCACCGAAAGACGATTAAAATCATTTTTTAGTGATAAAAAACATGACAGCGATGACAGG‑TTATTTTAACTAGTAGTGAACGCACACCAATgaga                                                        >  2:1016195/1‑100 (MQ=255)
                                              aCCGCGTACCGAAAGACGATTAATTTCGTTGTTTAGTCATAAAAAACCTGCCAGCGATGGCA‑GTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                       <  2:755609/100‑1 (MQ=255)
                                              aCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                       <  2:610004/100‑1 (MQ=255)
                                              aCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                       <  2:2827211/100‑1 (MQ=255)
                                              aCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                       <  1:3598677/100‑1 (MQ=255)
                                                  cgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtc                                                   <  2:73504/100‑1 (MQ=255)
                                                    cACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtctc                                                 >  2:2445727/1‑100 (MQ=255)
                                                     aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTcc                                                >  1:599722/1‑100 (MQ=255)
                                                        ggaaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTT‑CCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCtt                                             >  1:101655/3‑100 (MQ=255)
                                                        gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCtt                                             >  2:3807252/1‑100 (MQ=255)
                                                        gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGAGCACACCAATGCGAGTATATCCCtt                                             >  1:767126/1‑100 (MQ=255)
                                                          aaGACGATTAATTTCGTTTTTTAGTCATAAAAAATCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtat                                           <  1:3447621/100‑1 (MQ=255)
                                                              cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAg                                        <  2:3719003/99‑1 (MQ=255)
                                                              cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGt                                       <  1:2922220/100‑1 (MQ=255)
                                                              cGATTAATCTCGTTTTTTAGTCATAAAAAACCTTCCAGCGATGGAAAG‑TTTTTGTTGCATGTCGGGAGAGCACCCCAATACGAGTATCTCCCGtatacg                                        >  2:393617/1‑97 (MQ=255)
                                                                 ttttttttGTTTTTTTGTTATTATAATTTTGTCATCGGTGGCAGG‑TTTTTTTTGCCGGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTc                                      <  1:533449/94‑1 (MQ=255)
                                                                 ttAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtg                                    <  1:1008336/100‑1 (MQ=255)
                                                                        cGTTTTTTAGTCATAAAAATCCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                             <  1:2355386/100‑1 (MQ=255)
                                                                        cGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                             <  1:2138506/100‑1 (MQ=255)
                                                                         gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                            <  2:3100558/100‑1 (MQ=255)
                                                                         gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                            >  1:2800691/1‑100 (MQ=255)
                                                                         gTTTTTTAGTCATAAAAAACCTGCCAGATATGTCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                            <  1:2287015/100‑1 (MQ=255)
                                                                                 taactAAAAAACCTGCCAGCAATGGAAGG‑TTATTTTTGCCTGTCGTGAGTGCAAACCAATGCGAGTATCTCTCTTATAAGTATGTTGTTAAGTCAATag                     >  2:1660230/5‑99 (MQ=255)
                                                                                 gTCATAAAAAACCTGCCAGCGAGGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAGCTCCCTTATAAGTCTGTGGTTACGTCAATag                     >  2:3898826/1‑99 (MQ=255)
                                                                                  tCATAAAAAGCCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                   <  2:3437656/100‑1 (MQ=255)
                                                                                  tCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGGGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                    >  2:1950276/1‑99 (MQ=255)
                                                                                   cATAAAAAACCTGTCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAATTCTGTGGTTAAGTCAATata                    >  2:549861/1‑96 (MQ=255)
                                                                                   cATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                  <  2:2479935/100‑1 (MQ=255)
                                                                                   cATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGACAATagaaa                  >  1:3680476/1‑98 (MQ=255)
                                                                                    aaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                 >  1:3344987/3‑100 (MQ=38)
                                                                                     tAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                  >  1:2837469/1‑98 (MQ=255)
                                                                                     aaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                >  1:3730464/2‑100 (MQ=38)
                                                                                     aaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                >  1:2725461/2‑100 (MQ=38)
                                                                                     aaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTAACGTCAATAGAGAGc                >  1:3875849/2‑100 (MQ=38)
                                                                                      aaaaaaCCTGCCAGCGATGGCAGGTTTTTTTTT‑TTTTTCTTGAGCGCACAACAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCt               >  2:191447/1‑100 (MQ=255)
                                                                                        aaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt             >  1:2670077/1‑100 (MQ=39)
                                                                                        aaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt             >  1:2570555/1‑100 (MQ=38)
                                                                                        aaaaCCTGCCATCGCTGGCAGG‑TTTTTTATGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt             >  1:1625109/1‑100 (MQ=38)
                                                                                             cTGCCAGCGATGGCAGG‑TTTTTTTTGCTTGTCGTGAGCGCACACCAATGCGAGTATCTCTCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcac        >  2:2556234/1‑100 (MQ=255)
                                                                                              tGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcaca       >  1:818201/1‑100 (MQ=255)
                                                                                               gcaaGCGACGGCAGG‑TTTTTTTTGCATGTCGTGAGCGCACACCAATGCGAATACATCCCTTACAAATCTATGGTTACGTCAGTAGACAGCTTTGTgaca       >  2:2724732/4‑99 (MQ=255)
                                                                                               gCCAGCGCTTGCTGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGTCAGTATCTCGCTTATAAGTCTTTTTTTACTTCAATAGAGAGCTTTGt           <  2:3028774/95‑1 (MQ=255)
                                                                                               gCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcataa      >  2:2163903/1‑97 (MQ=255)
                                                                                                ccAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt     >  2:3286128/1‑100 (MQ=255)
                                                                                                 cAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAtt    >  2:1757853/1‑100 (MQ=255)
                                                                                                 cAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCACTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAtt    >  1:90162/1‑100 (MQ=255)
                                                                                                   gCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGTTTACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt     <  1:186089/97‑1 (MQ=255)
                                                                                                   |                                                                                                    
CGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTAT  >  NC_000913/4060162‑4060361

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