Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 721,004 T→A 17.9% 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.9% 368.2 / 66.0 152intergenic (+164/+52)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base T (56/67);  new base A (16/11);  total (72/80)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 158 total aligned reads displayed.

AAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGC  >  NC_000913/720905‑721102
                                                                                                   |                                                                                                  
aaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtat                                                                                                    >  1:3699712/1‑100 (MQ=255)
   tGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtattta                                                                                                 <  1:3046891/100‑1 (MQ=255)
     aTCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttatt                                                                                               <  1:2323442/100‑1 (MQ=255)
      tCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttatta                                                                                              <  2:2123502/100‑1 (MQ=255)
         aTCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattatta                                                                                           <  2:1114260/100‑1 (MQ=255)
           cTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACg                                                                                         <  2:866270/100‑1 (MQ=255)
             tAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                                                                                       >  1:1322392/1‑100 (MQ=255)
             tAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                                                                                       >  1:2082005/1‑100 (MQ=255)
              aaTTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTAAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCt                                                                                      <  2:2051546/100‑1 (MQ=255)
               atttatttAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                                                                                     <  1:874371/100‑1 (MQ=255)
                tttatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCttt                                                                                    >  1:1768455/1‑100 (MQ=255)
                tttatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCttt                                                                                    <  1:3159952/100‑1 (MQ=38)
                 ttatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttt                                                                                   <  2:685491/100‑2 (MQ=38)
                   atttaATGATGAGCTTTTTACTGAGTAATATAAAATATTGAATTGTTATTGTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCttttta                                                                                 <  2:1530453/100‑4 (MQ=37)
                    tttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                <  2:2072970/100‑1 (MQ=255)
                    tttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttan                                                                                <  2:1559050/100‑5 (MQ=37)
                     ttaATGATGAGCTTTTTACTTAGTAATATAAATTATTGAATTGTTATTTTTGTGTGTTGTTTAAGGTAAAAAACCGTATTTATTATTACGGCTttaatta                                                                               <  1:1436264/100‑1 (MQ=255)
                      taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                                                                              <  2:3681405/100‑1 (MQ=255)
                      taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                                                                              >  2:2606144/1‑100 (MQ=255)
                      taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                              <  1:2045797/100‑7 (MQ=25)
                       aaTGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCTTATTTATTATTACGGCTTtaattaat                                                                             <  1:2513465/100‑1 (MQ=255)
                       aaTGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                                                                             <  1:3155074/100‑1 (MQ=255)
                       aaTGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatctt                                                                             <  2:3321099/100‑8 (MQ=25)
                       aaTGATGAGCTGCTTACTTAGTACTATACAATATTGAATTGTTATTTTTGTGTGTTTTTTAAGATAAAAAGCCGTATTTATTTTTAAGGCTtta                                                                                   >  2:2353628/1‑94 (MQ=255)
                        atgattAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTTTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAGTACGGCTTTAATtagta                                                                            <  2:3910556/100‑1 (MQ=255)
                        atgatgAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                            >  1:1155716/1‑100 (MQ=255)
                          gatgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatcttaaa                                                                          <  2:1751511/100‑11 (MQ=21)
                            tgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAag                                                                        <  1:1991148/100‑1 (MQ=255)
                             gaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg                                                                       >  1:2505501/1‑100 (MQ=255)
                              aGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg                                                                       <  2:2952240/99‑1 (MQ=255)
                               gCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca                                                                     <  1:1217243/100‑1 (MQ=255)
                                cTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcag                                                                    >  1:1553393/1‑100 (MQ=255)
                                 tttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcagg                                                                   <  2:779298/100‑1 (MQ=255)
                                  ttttACTTAGTAATATAAAATATTGCATCGCTATTTTTGTGTTCTTTTAAAGATAAAAAGCCGTCTTTATTATTACTGCTTTAATTAATAAAaggcagtc                                                                  >  2:2573159/1‑98 (MQ=255)
                                  ttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcaggc                                                                  <  1:3302847/100‑1 (MQ=255)
                                     tACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGt                                                               <  1:1031481/100‑1 (MQ=255)
                                      aCTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTa                                                              >  2:3994749/1‑100 (MQ=255)
                                       cttaGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTATGGCTTTAATTAATAAAAGGAAGGCTGTAt                                                             >  2:3739287/4‑100 (MQ=255)
                                        ttaGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAtt                                                            <  2:2291399/98‑1 (MQ=255)
                                        ttaGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAt                                                             <  2:3965297/97‑1 (MQ=255)
                                         taGTAATATAAAATATTGAAGTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAGTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAtt                                                            <  2:1664735/98‑1 (MQ=255)
                                           gTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTaaa                                                         <  1:568499/100‑1 (MQ=255)
                                             aaTATAAAATTTTGAATTGTTATTTTTGTGTGTTGTTTAAGTTTATAGGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAt                                                       <  1:1309864/100‑1 (MQ=255)
                                             aaTATAAAATATTGAGTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCATCCTGTATTAAAAt                                                       >  2:1576925/1‑100 (MQ=255)
                                             aaTATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAt                                                       <  1:2419440/100‑1 (MQ=255)
                                             aaTATAAAATATTGAATTGTTATTTTTGTGTGTTCTTTAAGATAAAAAGCCGTATTTATAATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAt                                                       <  2:928879/100‑1 (MQ=255)
                                              atatAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAtt                                                      <  1:2917685/100‑1 (MQ=255)
                                                ataAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTaa                                                    <  2:3094434/100‑1 (MQ=255)
                                                 taAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTATATTAAAATTAat                                                   >  2:1504539/1‑100 (MQ=255)
                                                  aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                  <  1:980437/100‑1 (MQ=255)
                                                  aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                  <  1:2744427/100‑1 (MQ=255)
                                                  aaaaTATTGAATTGTTAGTTGTGTGTGTTGTTTAAGATAAAAATACGTATTTATTATTGCGGATTTAATTAATAAAAAGCAGGCTGTATTAAAATTAata                                                  >  2:1846862/1‑100 (MQ=255)
                                                   aaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAatat                                                 <  1:183699/100‑1 (MQ=255)
                                                   aaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAatat                                                 >  1:725721/1‑100 (MQ=255)
                                                       atTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCaa                                             <  1:2730547/100‑1 (MQ=255)
                                                         tGAATTGTTATTTTTGTGTGTTGTTTAAAATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                           <  2:2667603/100‑1 (MQ=255)
                                                          gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGc                                          >  2:1411619/1‑100 (MQ=255)
                                                          gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGc                                          <  1:1457840/100‑1 (MQ=255)
                                                          gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGc                                          <  1:1741658/100‑1 (MQ=255)
                                                          gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAACGCAGGCTGTATTAAAATTAATATTCAAAGc                                          >  2:1369333/1‑100 (MQ=255)
                                                          gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCATATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAACATTAATATTCAAAGc                                          >  2:3562876/1‑100 (MQ=255)
                                                           aaTTTTTATTTTTGTGTGTTGTTTAAGATATAAAGCTGTATTTATTATTAGGGCTTTAATTAATAAAGGGCAGGCGGTTTTAATATTAATATTCAAAg                                           <  2:848374/98‑1 (MQ=255)
                                                           aaTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCa                                         >  1:2445185/1‑100 (MQ=255)
                                                           aaTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCa                                         >  1:2975562/1‑100 (MQ=255)
                                                             ttGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATa                                       >  1:570909/1‑100 (MQ=255)
                                                               gTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATaaa                                     >  2:1990484/1‑100 (MQ=255)
                                                               gTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATaaa                                     >  2:177697/1‑100 (MQ=255)
                                                                ttATTTTTTTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAc                                    <  2:1983621/100‑1 (MQ=255)
                                                                ttATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAc                                    <  2:3379354/100‑1 (MQ=255)
                                                                 tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                   >  1:437589/1‑100 (MQ=255)
                                                                 tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAcc                                   <  2:1661200/100‑1 (MQ=255)
                                                                 tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGAAGGCTGTATTAAAATTAATATTCAAAGCATAAAc                                    >  2:463741/1‑99 (MQ=255)
                                                                   tttttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGa                                 >  2:515659/1‑100 (MQ=255)
                                                                            gttgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATa                        >  1:1536184/1‑100 (MQ=255)
                                                                              tgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAc                       >  2:1189246/1‑99 (MQ=255)
                                                                               gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCg                     >  1:1362585/1‑100 (MQ=255)
                                                                               gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCg                     <  2:3099271/100‑1 (MQ=255)
                                                                                tttattaattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCgg                    >  1:3583058/11‑100 (MQ=25)
                                                                                tttAAGATAAAAAGCCGTATTTATTATTACGACTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCgn                    <  2:305770/100‑2 (MQ=255)
                                                                                 ttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                   >  2:2341995/1‑100 (MQ=255)
                                                                                  tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                  >  1:1096559/1‑100 (MQ=255)
                                                                                  tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATAGTCAAAGCATAAACCGATCGCCAATACCGGtt                  <  1:2782372/100‑1 (MQ=255)
                                                                                   attaattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGttt                 >  2:2627164/8‑100 (MQ=25)
                                                                                     gATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                <  2:618896/99‑1 (MQ=255)
                                                                                      aTAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg              <  2:4124461/100‑1 (MQ=255)
                                                                                       tAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt             >  1:3015112/1‑100 (MQ=255)
                                                                                       tAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt             >  1:3485727/1‑100 (MQ=255)
                                                                                       attaaaGCCGTAATAATAAATACGGCTTTACTTAATAAAAGGCAGGCTGTAATAAAATTAATATTCAAAGCATAAACCGATAGCCAAAACCTGTTTCact             >  2:3418106/4‑98 (MQ=25)
                                                                                       attaaaGCCGTAATAATAAATACGGCTTTAATTAATACAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt             >  2:2213128/4‑100 (MQ=37)
                                                                                       attaaaGCCGTAATAATAAATACGGCTTTAATTAATACAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt             >  2:2216471/4‑100 (MQ=37)
                                                                                        ttaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTa            >  2:385352/3‑100 (MQ=37)
                                                                                        aaaaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAACGGTTTCAGTa            >  2:318112/1‑100 (MQ=255)
                                                                                         taaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa           >  2:1285925/2‑100 (MQ=37)
                                                                                         aaaaGCCGTATTTATTATTACGTCTTTAATTAACAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATCGCCAATACCGATTTCAGt             >  2:2108212/1‑98 (MQ=38)
                                                                                          aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt          >  1:2206954/1‑100 (MQ=37)
                                                                                               cGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGaaa      <  2:4086079/99‑1 (MQ=255)
                                                                                                gTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCATTACTGaaacc    <  1:3604395/100‑3 (MQ=255)
                                                                                                gTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATg    >  1:3884857/1‑100 (MQ=255)
                                                                                                 tatttattATTACGGCTTTTATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATgg   >  2:2122228/1‑100 (MQ=255)
                                                                                                  atttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGc  >  1:3002130/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
AAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGC  >  NC_000913/720905‑721102

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