Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 721,004 T→A 17.5% intergenic (+164/+52) ybfK → / ← kdpE uncharacterized protein/response regulator in two‑component regulatory system with KdpD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913721,0040TA17.5% 479.4 / 88.0 194intergenic (+164/+52)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base T (70/89);  new base A (15/19);  total (85/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-02

Only 100 of 207 total aligned reads displayed.

AAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCG  >  NC_000913/720905‑721103
                                                                                                   |                                                                                                   
 aaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatt                                                                                                    >  2:3366860/1‑100 (MQ=255)
 aaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatt                                                                                                    <  2:1242886/100‑1 (MQ=255)
   tGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtattta                                                                                                  <  2:2354164/100‑1 (MQ=255)
    gATCAATCTTAATTTATTTAATGATGAGCTTTTTACTAAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGACGtatttat                                                                                                 >  1:4165275/1‑100 (MQ=255)
      tCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAATCCGtatttttt                                                                                                >  1:741850/1‑99 (MQ=255)
      tCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttatta                                                                                               <  2:1110427/100‑1 (MQ=255)
         aTCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAATATATTGTATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattatta                                                                                            >  1:2395748/1‑100 (MQ=255)
           cTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACg                                                                                          <  1:3801842/100‑1 (MQ=255)
           cTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACg                                                                                          >  2:768453/1‑100 (MQ=255)
             tcttttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACgg                                                                                         <  2:2998456/96‑1 (MQ=255)
             tAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                                                                                        <  1:311952/100‑1 (MQ=255)
              aaTTTATTTAATGATGAGCTTTTTTTTTTTTTTTTTAAAATATTGAATTGTTATTTTTTTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCt                                                                                       >  1:2740809/1‑100 (MQ=255)
              aaTTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCt                                                                                       >  2:1332767/1‑100 (MQ=255)
               atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtt                                                                                      <  2:697379/100‑1 (MQ=255)
               atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtn                                                                                      <  1:1861818/100‑2 (MQ=255)
                tttatttaATTATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                                                                                      <  1:1480602/99‑1 (MQ=255)
                 ttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTtta                                                                                    <  1:617866/100‑1 (MQ=255)
                 ttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTtta                                                                                    <  2:3504300/100‑1 (MQ=255)
                 ttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttt                                                                                    <  2:3808052/100‑2 (MQ=38)
                  tatttaATGATGAGCTTTTTACTAAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaa                                                                                   <  1:1042226/100‑1 (MQ=255)
                  tatttaAAGATGATCTTTTTAATCAGTAATATAAAATATTGAATTATTTTTTTTGTGTGTTATTAAAGAAAAAAAGCCGTATTTATTATTACGGCTtta                                                                                    <  1:4309277/99‑1 (MQ=255)
                   atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCTTATTTATTATTACGTCTctaat                                                                                  >  1:1143330/1‑100 (MQ=255)
                   atttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCttttta                                                                                  <  1:3195732/100‑4 (MQ=38)
                    tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                 <  1:1152030/100‑1 (MQ=255)
                    tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                 <  1:1705608/100‑1 (MQ=255)
                    tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                 <  2:1805101/100‑1 (MQ=255)
                     ttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatc                                                                                <  1:676403/100‑6 (MQ=37)
                      taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                               <  1:4075733/100‑7 (MQ=37)
                      taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                               <  1:2961621/100‑7 (MQ=37)
                      taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                               <  2:2031368/100‑7 (MQ=37)
                       aaTGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                                                                              <  2:3207513/100‑1 (MQ=255)
                       aaTGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatctt                                                                              <  2:4361417/100‑8 (MQ=25)
                         tgatgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatcttaa                                                                            <  2:657424/100‑10 (MQ=25)
                         tgatgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCAGTATTTATTATTACGGCTTTAATtaataa                                                                            >  2:357652/1‑100 (MQ=255)
                             gaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg                                                                        >  1:1733528/1‑100 (MQ=255)
                              aGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggc                                                                       <  2:3694452/100‑1 (MQ=255)
                              aGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggc                                                                       <  1:556763/100‑1 (MQ=255)
                              aGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggc                                                                       <  1:1230443/100‑1 (MQ=255)
                               gCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca                                                                      >  2:2176090/1‑100 (MQ=255)
                                 tttttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTCATAAAaggcagg                                                                    <  2:2886794/100‑1 (MQ=255)
                                     tACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGt                                                                <  2:4334277/100‑1 (MQ=255)
                                      aCTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTa                                                               <  2:304183/100‑1 (MQ=255)
                                       cTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTan                                                              <  1:2071140/100‑2 (MQ=255)
                                        tCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAtt                                                             <  2:2273999/100‑1 (MQ=255)
                                          aGTAATATAAAATATTGAATTGTTATGTATGTGTGTTGTTTAAGATAAAAATCGGTATTTATTATTATATCTTGAATTAATAAAAGGCAGGGAGTATTaa                                                           >  1:4329203/1‑100 (MQ=255)
                                           gTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTaaa                                                          <  1:2500329/100‑1 (MQ=255)
                                            tAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTaaaa                                                         <  1:2233829/100‑1 (MQ=255)
                                             aaTATAAAAGATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAt                                                        <  1:1464814/100‑1 (MQ=255)
                                              atatAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAtt                                                       <  1:95462/100‑1 (MQ=255)
                                               tataAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTa                                                      <  1:2254286/100‑1 (MQ=255)
                                               tataAAACATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTa                                                      >  2:2305248/1‑100 (MQ=255)
                                                ataAAATATTGAATTGTTATTTTTTTTTTTTTTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTTAAATTaa                                                     >  1:214463/1‑100 (MQ=255)
                                                 taAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAat                                                    <  1:1999722/100‑1 (MQ=255)
                                                  aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                   >  2:255110/1‑100 (MQ=255)
                                                  aaaaTATTGAATTGTTATTTTTGCGTGTTGTTTAAGATAAAAAGCCGTATTTATTTTTACGGCTTTAATTAATAAAAGGCAGTCTGTATTAAAATTAata                                                   >  1:892009/1‑100 (MQ=255)
                                                     atatTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTc                                                >  1:2788603/1‑100 (MQ=255)
                                                      tatTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCa                                               >  1:4225702/1‑100 (MQ=255)
                                                        ttGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCaaa                                             >  1:3328187/1‑100 (MQ=255)
                                                         tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                            <  2:4072064/100‑1 (MQ=255)
                                                         tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                            <  1:2655246/100‑1 (MQ=255)
                                                          gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGc                                           <  2:4000635/100‑1 (MQ=255)
                                                            aTTGTTATTTTTGTGTGTAGTTAAAAAAAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCa                                          <  1:1773109/99‑1 (MQ=255)
                                                               gTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATaaa                                      >  1:724629/1‑100 (MQ=255)
                                                               gTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATaaa                                      >  1:1719192/1‑100 (MQ=255)
                                                               gTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATaaa                                      <  2:2407390/100‑1 (MQ=255)
                                                                 tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                    >  2:4098269/1‑100 (MQ=255)
                                                                 tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                    >  1:3976280/1‑100 (MQ=255)
                                                                 tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATGACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                    <  1:3080347/100‑1 (MQ=255)
                                                                      ttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAg                               <  1:377518/100‑1 (MQ=255)
                                                                       tgtgttTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGc                              <  1:2302235/100‑1 (MQ=255)
                                                                         tgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCa                            >  2:1031335/1‑100 (MQ=255)
                                                                          gtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCaa                           >  2:716065/1‑100 (MQ=255)
                                                                           tgttgtTTAAGATAAAAAGCCGTATTTATTATTACGACTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAAt                          >  2:3093260/1‑100 (MQ=255)
                                                                           tgttgtTTAAGATAAAAAGCCGTATTTATCATTACGGCCTTAATTAATAAAAGGCATGCTGTATTAAAATTACTATTCAAAGCATAAACCGATAGc                              >  1:3468197/1‑96 (MQ=255)
                                                                            tttgctTAAGATAAAACGTAGTAGGTATTATTACGGCTTTAATTAATTCAGGGAAGGATGGACTGAAATTTATATTTAAATCATAAACCGAGAGcc                             <  1:1332233/95‑1 (MQ=255)
                                                                            gttgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATa                         <  2:1733534/100‑1 (MQ=255)
                                                                               gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAcc                       >  1:503970/1‑99 (MQ=255)
                                                                               gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCg                      >  2:1573809/1‑100 (MQ=255)
                                                                               gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCg                      >  2:1293870/1‑100 (MQ=255)
                                                                               gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCg                      >  2:1146131/1‑100 (MQ=255)
                                                                                tttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCgg                     >  2:4208544/1‑100 (MQ=255)
                                                                                 ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCgg                     <  1:1162324/99‑1 (MQ=255)
                                                                                 ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                    >  2:3919423/1‑100 (MQ=255)
                                                                                 ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAAAATTCAAAGCATAAACCGATAGCCAATACCgg                     >  2:4306767/1‑99 (MQ=255)
                                                                                  tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                    <  2:2700728/99‑1 (MQ=255)
                                                                                  aaaaaaaaaaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAcn                       <  1:284122/90‑2 (MQ=255)
                                                                                   aaGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGttt                  >  2:3313910/1‑100 (MQ=255)
                                                                                     gATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCa                >  2:110031/1‑100 (MQ=255)
                                                                                     gATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCa                >  2:19304/1‑100 (MQ=255)
                                                                                      aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg               <  2:887015/100‑1 (MQ=255)
                                                                                      aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg               <  2:3015258/100‑1 (MQ=255)
                                                                                      aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg               <  1:3662100/100‑1 (MQ=255)
                                                                                       attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt              >  2:2921590/4‑100 (MQ=37)
                                                                                       attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAACCAATACCGGTTTCAg               >  2:1427061/4‑99 (MQ=37)
                                                                                         aaaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCACAGCATAAACCGATAGCCAATACCGGTTTAAGTaa            >  1:127058/1‑100 (MQ=255)
                                                                                         aaaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa            >  2:1833370/1‑100 (MQ=255)
                                                                                          aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt           >  1:2539266/1‑100 (MQ=37)
                                                                                          aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt           >  2:4401222/1‑100 (MQ=37)
                                                                                              ccGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGaaa       <  2:2483742/100‑1 (MQ=255)
                                                                                                  atttattATTACGCCTTTAATTAATTAACGGCAGGCTGTATTAAAATTGATATCAAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAAggg    <  1:1882238/99‑1 (MQ=25)
                                                                                                   |                                                                                                   
AAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCG  >  NC_000913/720905‑721103

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: