Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,060,264 A→C 15.8% 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.8% 477.1 / 86.3 188intergenic (+34/‑183)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (66/92);  new base C (17/13);  total (83/105)
Fisher's exact test for biased strand distribution p-value = 1.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.57e-01

Only 100 of 189 total aligned reads displayed.

TTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCT  >  NC_000913/4060166‑4060363
                                                                                                  |                                                                                                    
ttCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGAt                                                                                                     >  1:2321276/1‑100 (MQ=255)
  cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATgg                                                                                                   >  1:1207451/1‑100 (MQ=255)
  cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATgg                                                                                                   >  1:2202073/1‑100 (MQ=255)
    tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGCTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCa                                                                                                 >  2:2570616/1‑100 (MQ=255)
    tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCa                                                                                                 <  2:3602423/100‑1 (MQ=255)
    tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCa                                                                                                 >  2:2940554/1‑100 (MQ=255)
    tAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCCGCCAGCGATGGCa                                                                                                 <  2:2538140/100‑1 (MQ=255)
     aaaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAg                                                                                                <  2:189566/100‑1 (MQ=255)
       aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAggn                                                                                              <  2:1140073/100‑2 (MQ=255)
       aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGt                                                                                             <  1:1469290/100‑1 (MQ=255)
       aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGt                                                                                             <  2:553040/100‑1 (MQ=255)
       aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACAGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGt                                                                                             >  1:1994503/1‑100 (MQ=255)
        cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtt                                                                                            <  2:827315/100‑1 (MQ=255)
          tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttt                                                                                          >  2:2740296/1‑100 (MQ=255)
           cACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttt                                                                                         <  2:3090349/100‑1 (MQ=255)
             tcTGACGGAATACGATCGTCTCGTCTCCTACCTTAGAACTAAATAGGATTAATTTCGTTTTTTTGTCATAACATACCTGCCATCGCTGGCATGtgtttga                                                                                       <  1:2628024/99‑3 (MQ=255)
             ccggACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCAGCGATGGCAGGtttttttt                                                                                      <  1:1108300/97‑1 (MQ=255)
             ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                       >  2:3183332/1‑100 (MQ=255)
             ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                       <  2:1789672/100‑1 (MQ=255)
             ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                       <  1:1912054/100‑1 (MQ=255)
             ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGAACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttttt                                                                                        <  1:758998/99‑1 (MQ=255)
              cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttttt                                                                                      <  1:1624306/100‑1 (MQ=255)
               tGACGGAAAACGATCGTCGCCGCTTCAACTGTTTATCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatg                                                                                     <  2:59483/100‑4 (MQ=21)
               tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttn                                                                                     <  1:276393/100‑2 (MQ=255)
                gACGGAAAACGATCGTCGCCGCGCTAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTg                                                                                    <  1:3024638/100‑1 (MQ=255)
                gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatga                                                                                    <  1:3390481/100‑5 (MQ=38)
                gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTg                                                                                    <  1:3492672/100‑1 (MQ=255)
                gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTg                                                                                    <  1:3462869/100‑1 (MQ=255)
                  cgatctaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGcc                                                                                  <  2:2430683/94‑1 (MQ=255)
                  ccgAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatgact                                                                                  <  1:2261470/98‑7 (MQ=37)
                  cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGcc                                                                                  >  2:1643818/1‑100 (MQ=255)
                  cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGcc                                                                                  <  1:61324/100‑1 (MQ=255)
                   ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCt                                                                                 >  2:2082131/1‑100 (MQ=255)
                   ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAAAAAAAACCTGCCATCGCTGGCAGGttttttatgacta                                                                                 <  1:3586074/100‑8 (MQ=25)
                    gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTg                                                                                >  1:2670041/1‑100 (MQ=255)
                    gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTg                                                                                >  1:3028158/1‑100 (MQ=255)
                        aCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGccagatcg                                                                            >  1:2430669/1‑94 (MQ=255)
                        aCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGt                                                                            >  2:3324197/1‑100 (MQ=255)
                         cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTg                                                                           <  2:2952664/100‑1 (MQ=255)
                         cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTg                                                                           >  2:2154122/1‑100 (MQ=255)
                         cGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTg                                                                           <  1:2740277/100‑1 (MQ=255)
                           aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAg                                                                         <  2:1679115/100‑1 (MQ=255)
                           aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAg                                                                         <  2:1604112/100‑1 (MQ=255)
                            tcgacgCCGCGCCCACCGCGCACCGAAACACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAg                                                                         <  1:361884/100‑1 (MQ=255)
                             cgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAACCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAgcg                                                                       <  1:1058609/100‑1 (MQ=255)
                              gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATGTCGTTGTTTAGTCATAAACAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCGCTCGTAAacgc                                                                      >  2:1407513/1‑100 (MQ=255)
                               tcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGca                                                                     <  2:3202174/100‑1 (MQ=255)
                                 gccgcGCCAAGCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAAAAAAAACCTGCCAGCGATGGCAGTTTTTTTTTTGCCTGTCGTGAGCGcaca                                                                   <  2:2968804/100‑1 (MQ=255)
                                 gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGcaca                                                                   >  2:3507306/1‑100 (MQ=255)
                                   cgcgCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGcaca                                                                   >  1:2494183/1‑98 (MQ=255)
                                   cgcgCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                                 >  2:505118/1‑100 (MQ=255)
                                   cgcgCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                                 <  2:2320294/100‑1 (MQ=255)
                                      gcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAAt                                                              >  1:3587497/1‑100 (MQ=255)
                                        cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGc                                                            <  2:1216700/100‑1 (MQ=255)
                                        cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGc                                                            >  2:3441498/1‑100 (MQ=255)
                                         aaCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCg                                                           <  2:1774655/100‑1 (MQ=255)
                                           ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAACCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                         <  2:436566/100‑1 (MQ=255)
                                            cgcgCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGt                                                        <  1:2066069/100‑1 (MQ=255)
                                             gcgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAACCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTa                                                       <  1:2970914/100‑1 (MQ=255)
                                             gcgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTa                                                       >  1:3109929/1‑100 (MQ=255)
                                              cgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAt                                                      >  1:1748973/1‑100 (MQ=255)
                                               gcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtc                                                     <  1:658838/100‑1 (MQ=255)
                                                cACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtct                                                    >  1:1257535/1‑100 (MQ=255)
                                                  ccGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTcc                                                  >  2:1347312/1‑100 (MQ=255)
                                                        gACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTAGTGAGCGCAGACCAATGCGAGTAtctacat                                                >  1:132157/1‑92 (MQ=255)
                                                         aCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATaa                                           <  2:3728979/100‑1 (MQ=255)
                                                          ctaTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGATGTCATGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAg                                          <  2:2631641/98‑1 (MQ=255)
                                                            aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTc                                        <  1:413669/100‑1 (MQ=255)
                                                                aTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                    <  2:2706381/100‑1 (MQ=255)
                                                                 tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                    >  1:975905/1‑99 (MQ=255)
                                                                 tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                   <  2:3103607/100‑1 (MQ=255)
                                                                 tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                   >  2:1960201/1‑100 (MQ=255)
                                                                 tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                   <  2:942279/100‑1 (MQ=255)
                                                                  ttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGt                                  <  2:1748980/100‑1 (MQ=255)
                                                                     gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                               <  2:1488005/100‑1 (MQ=255)
                                                                      ttttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                              <  2:3149691/100‑1 (MQ=255)
                                                                      ttttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                              <  2:1876176/100‑1 (MQ=255)
                                                                        ttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                            <  2:3250113/100‑1 (MQ=255)
                                                                           tAGTCATAAAACACCGGCCCGCGAGGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAAt                         <  1:879242/100‑1 (MQ=255)
                                                                            aGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                        <  2:201495/100‑1 (MQ=255)
                                                                             gcaaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                       >  2:723761/6‑100 (MQ=38)
                                                                             gcaaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                        >  2:1119899/6‑99 (MQ=37)
                                                                              caaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                      >  2:280689/5‑100 (MQ=38)
                                                                               cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTTTTTTTGTAAGCGAAAAACAAATCTATTATATCACTTATAGGTTTTTGGGTACGTCAATaga                      >  1:1101294/1‑99 (MQ=255)
                                                                                  aaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                  >  2:3416846/1‑100 (MQ=255)
                                                                                  aaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                  >  1:1433888/1‑100 (MQ=255)
                                                                                  aaaaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCGCCCTTATAAGTCTGTGGTTACGTCAATagagag                   <  1:1140744/99‑1 (MQ=255)
                                                                                    aaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt                >  2:3188343/1‑100 (MQ=255)
                                                                                    aaaaCCTGCCAGCGCTGGCAGGTTTTTTTTTGACTGTCGTTAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt                >  2:2692103/1‑100 (MQ=255)
                                                                                     aaaCCTGCCATCGCTGGCAGGTTTTTTTTTTTTTTTTTTCAGCTCACACCAATGCGAGTATCTCCCTTATAAATGTGTCGTTACGTTAAtttaatcctat               >  1:371467/1‑90 (MQ=12)
                                                                                     aaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt               >  2:3610737/1‑100 (MQ=255)
                                                                                          tGCCAGCGATGGCAGGTTTTTTTTTTTTTTTTTTTTTTGCACACAAATGCGAGTATCTCTATTTTAAGTTTGTGGTTACGTCAATATAGAGCTTTGtta           >  1:1858272/1‑97 (MQ=255)
                                                                                           gCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcaca         <  1:3690376/100‑1 (MQ=255)
                                                                                             cAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTCTAAGTCAGTGGTTACGTCAATAGAGACCATTGTCACaa        >  1:2363473/1‑99 (MQ=255)
                                                                                             cAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt       <  1:3089894/100‑1 (MQ=255)
                                                                                              aGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt       <  2:2202084/99‑1 (MQ=255)
                                                                                                cGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTa     <  1:2667790/99‑1 (MQ=255)
                                                                                                 gATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATc   <  2:3427363/100‑1 (MQ=255)
                                                                                                 gATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATc   <  2:2651356/100‑1 (MQ=255)
                                                                                                  aTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCt  <  2:1828496/100‑1 (MQ=255)
                                                                                                  |                                                                                                    
TTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCT  >  NC_000913/4060166‑4060363

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: