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)
                       ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTg                                                                              >  1:3379183/1‑100 (MQ=255)
                       ggAAAACGAGCGGCGCCGGTCCAGCCGCGTAACGAAAGCCGATTAATTTCGTGTGTTAGTCATAAAATACCTGCCAGCGATGGCAGGttttttt                                                                                     <  2:3236732/94‑1 (MQ=255)
                        gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGt                                                                             >  1:1733493/1‑100 (MQ=255)
                              gATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAg                                                                       <  2:1699289/100‑1 (MQ=255)
                               cgatctgccgcGCCAACCGCGCACCGAAAGTCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAgc                                                                      <  1:2995110/94‑1 (MQ=255)
                                tcgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAgcg                                                                     <  2:1093363/100‑1 (MQ=255)
                                tcgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCATGTCGTGAgcg                                                                     <  2:1638543/100‑1 (MQ=255)
                                tcgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCGGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAg                                                                       <  1:2753996/99‑1 (MQ=255)
                                 cgtcgCCGCGCCAACCGCGCACTGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAgcgc                                                                    <  1:1485521/100‑1 (MQ=255)
                                  gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTAAAAGCAATGGCAGGAATTTTTTGCCTGTCGTGAGCGca                                                                   >  1:2196104/1‑100 (MQ=255)
                                   tcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGcac                                                                  <  1:2257241/100‑1 (MQ=255)
                                    cgccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGcaca                                                                 <  1:1103681/100‑1 (MQ=255)
                                    cgccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGcaca                                                                 >  2:47332/1‑100 (MQ=255)
                                     gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAgcagatcg                                                                >  2:2995083/1‑94 (MQ=255)
                                     gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCacac                                                                >  2:829087/1‑100 (MQ=255)
                                      ccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                               <  2:2918193/100‑1 (MQ=255)
                                        gcgcCAACCGCGCAACGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCaa                                                             >  2:927589/1‑100 (MQ=255)
                                           ccATCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATg                                                           <  1:595077/99‑1 (MQ=255)
                                           ccAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCATGTAGTGAGCGCACACCAATGc                                                          >  1:3114709/1‑100 (MQ=255)
                                            cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGc                                                          >  2:2024510/1‑99 (MQ=255)
                                            cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCg                                                         >  2:58965/1‑100 (MQ=255)
                                            cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGACTGTCGTGAGCGCACACCAATGCg                                                         >  2:2138487/1‑100 (MQ=255)
                                             aaCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                        >  2:1824498/1‑100 (MQ=255)
                                             aaCCGCGCACCGAAAGACGATTAAAATCATTTTTTAGTGATAAAAAACATGACAGCGATGACAGGTTATTTTAACTAGTAGTGAACGCACACCAATgaga                                                        >  2:1016195/1‑100 (MQ=255)
                                              aCCGCGTACCGAAAGACGATTAATTTCGTTGTTTAGTCATAAAAAACCTGCCAGCGATGGCAGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                       <  2:755609/100‑1 (MQ=255)
                                              aCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                       <  2:610004/100‑1 (MQ=255)
                                              aCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                       <  2:2827211/100‑1 (MQ=255)
                                              aCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                       <  1:3598677/100‑1 (MQ=255)
                                                  cgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtc                                                   <  2:73504/100‑1 (MQ=255)
                                                    cACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtctc                                                 >  2:2445727/1‑100 (MQ=255)
                                                     aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTcc                                                >  1:599722/1‑100 (MQ=255)
                                                        ggaaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCtt                                             >  1:101655/3‑100 (MQ=255)
                                                        gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCtt                                             >  2:3807252/1‑100 (MQ=255)
                                                        gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGAGCACACCAATGCGAGTATATCCCtt                                             >  1:767126/1‑100 (MQ=255)
                                                          aaGACGATTAATTTCGTTTTTTAGTCATAAAAAATCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtat                                           <  1:3447621/100‑1 (MQ=255)
                                                              cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAg                                        <  2:3719003/99‑1 (MQ=255)
                                                              cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGt                                       <  1:2922220/100‑1 (MQ=255)
                                                              cGATTAATCTCGTTTTTTAGTCATAAAAAACCTTCCAGCGATGGAAAGTTTTTGTTGCATGTCGGGAGAGCACCCCAATACGAGTATCTCCCGtatacg                                        >  2:393617/1‑97 (MQ=255)
                                                                 ttttttttGTTTTTTTGTTATTATAATTTTGTCATCGGTGGCAGGTTTTTTTTGCCGGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTc                                      <  1:533449/94‑1 (MQ=255)
                                                                 ttAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtg                                    <  1:1008336/100‑1 (MQ=255)
                                                                        cGTTTTTTAGTCATAAAAATCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                             <  1:2355386/100‑1 (MQ=255)
                                                                        cGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                             <  1:2138506/100‑1 (MQ=255)
                                                                         gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                            <  2:3100558/100‑1 (MQ=255)
                                                                         gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                            >  1:2800691/1‑100 (MQ=255)
                                                                         gTTTTTTAGTCATAAAAAACCTGCCAGATATGTCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                            <  1:2287015/100‑1 (MQ=255)
                                                                                 taactAAAAAACCTGCCAGCAATGGAAGGTTATTTTTGCCTGTCGTGAGTGCAAACCAATGCGAGTATCTCTCTTATAAGTATGTTGTTAAGTCAATag                     >  2:1660230/5‑99 (MQ=255)
                                                                                 gTCATAAAAAACCTGCCAGCGAGGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAGCTCCCTTATAAGTCTGTGGTTACGTCAATag                     >  2:3898826/1‑99 (MQ=255)
                                                                                  tCATAAAAAGCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                   <  2:3437656/100‑1 (MQ=255)
                                                                                  tCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGGGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                    >  2:1950276/1‑99 (MQ=255)
                                                                                   cATAAAAAACCTGTCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAATTCTGTGGTTAAGTCAATata                    >  2:549861/1‑96 (MQ=255)
                                                                                   cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                  <  2:2479935/100‑1 (MQ=255)
                                                                                   cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGACAATagaaa                  >  1:3680476/1‑98 (MQ=255)
                                                                                    aaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                 >  1:3344987/3‑100 (MQ=38)
                                                                                     tAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                  >  1:2837469/1‑98 (MQ=255)
                                                                                     aaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                >  1:3730464/2‑100 (MQ=38)
                                                                                     aaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                >  1:2725461/2‑100 (MQ=38)
                                                                                     aaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTAACGTCAATAGAGAGc                >  1:3875849/2‑100 (MQ=38)
                                                                                      aaaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTTTTTTCTTGAGCGCACAACAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCt               >  2:191447/1‑100 (MQ=255)
                                                                                        aaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt             >  1:2670077/1‑100 (MQ=39)
                                                                                        aaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt             >  1:2570555/1‑100 (MQ=38)
                                                                                        aaaaCCTGCCATCGCTGGCAGGTTTTTTATGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt             >  1:1625109/1‑100 (MQ=38)
                                                                                             cTGCCAGCGATGGCAGGTTTTTTTTGCTTGTCGTGAGCGCACACCAATGCGAGTATCTCTCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcac        >  2:2556234/1‑100 (MQ=255)
                                                                                              tGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcaca       >  1:818201/1‑100 (MQ=255)
                                                                                               gcaaGCGACGGCAGGTTTTTTTTGCATGTCGTGAGCGCACACCAATGCGAATACATCCCTTACAAATCTATGGTTACGTCAGTAGACAGCTTTGTgaca       >  2:2724732/4‑99 (MQ=255)
                                                                                               gCCAGCGCTTGCTGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGTCAGTATCTCGCTTATAAGTCTTTTTTTACTTCAATAGAGAGCTTTGt           <  2:3028774/95‑1 (MQ=255)
                                                                                               gCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcataa      >  2:2163903/1‑97 (MQ=255)
                                                                                                ccAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt     >  2:3286128/1‑100 (MQ=255)
                                                                                                 cAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAtt    >  2:1757853/1‑100 (MQ=255)
                                                                                                 cAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCACTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAtt    >  1:90162/1‑100 (MQ=255)
                                                                                                   gCGATGGCAGGTTTTTTTTGCCTGTCGTGAGTTTACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt     <  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: