Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 721,006 T→A 15.9% intergenic (+166/+50) 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,0060TA15.9% 361.8 / 55.7 138intergenic (+166/+50)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base T (44/72);  new base A (10/12);  total (54/84)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

Only 100 of 153 total aligned reads displayed.

ATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTG  >  NC_000913/720907‑721105
                                                                                                   |                                                                                                   
aTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtattt                                                                                                     <  2:2255296/100‑1 (MQ=255)
 tGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtattta                                                                                                    <  2:1245179/100‑1 (MQ=255)
 tGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtattta                                                                                                    <  2:1252555/100‑1 (MQ=255)
  gATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttat                                                                                                   <  2:2232970/100‑1 (MQ=255)
  gATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttat                                                                                                   >  1:246029/1‑100 (MQ=255)
   aTCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttatt                                                                                                  <  2:3310055/100‑1 (MQ=255)
   aTCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttatt                                                                                                  <  2:1602881/100‑1 (MQ=255)
     cAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                                                                                                <  2:150535/100‑1 (MQ=255)
     cAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                                                                                                <  2:2801879/100‑1 (MQ=255)
     cAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                                                                                                >  1:3688690/1‑100 (MQ=255)
         cTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACg                                                                                            <  2:1151129/100‑1 (MQ=255)
         cTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACg                                                                                            <  1:3815168/100‑1 (MQ=255)
           tAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                                                                                          <  2:1736581/100‑1 (MQ=255)
              tttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCttt                                                                                       <  1:1050899/100‑1 (MQ=38)
               ttatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttt                                                                                      <  2:1972432/100‑2 (MQ=38)
                  tttaATGATGATCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttat                                                                                   <  2:2187086/100‑5 (MQ=37)
                  tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                                                                                   <  2:3321899/100‑1 (MQ=255)
                  tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttat                                                                                   <  1:1756991/100‑5 (MQ=38)
                    taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                                                                                 >  1:3155778/1‑100 (MQ=255)
                    taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                                 <  2:2635973/100‑7 (MQ=37)
                    taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                                                                                 <  2:2756019/100‑7 (MQ=37)
                     aaTGAGGAGCTTTTTACTCAGTAATATAAAAAATGGAATTGTTATTTTTGTGTGTGGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatctt                                                                                <  1:3245187/100‑8 (MQ=21)
                      atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                               <  2:3737110/100‑1 (MQ=255)
                      atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                               <  2:3197637/100‑1 (MQ=255)
                      atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                               >  1:2546662/1‑100 (MQ=255)
                      atgatgAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                                                                               <  1:1728791/100‑1 (MQ=255)
                       tgatgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTTTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaataa                                                                              <  1:1196757/100‑1 (MQ=255)
                         atgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaaa                                                                            >  1:2581732/1‑100 (MQ=255)
                          tgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTCTTATTACGGCTTTAATTAATaaaa                                                                            >  1:3754029/1‑99 (MQ=255)
                             gCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca                                                                        <  2:2460327/100‑1 (MQ=255)
                              cTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcag                                                                       >  2:290545/1‑100 (MQ=255)
                               tttttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTCTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca                                                                        >  1:2112695/1‑98 (MQ=255)
                                ttttACTAAGTAATATAAAATATTAAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcaggc                                                                     <  1:2255389/100‑1 (MQ=255)
                                     cTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTa                                                                 <  2:1119096/99‑1 (MQ=255)
                                     cTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTAt                                                                <  1:738800/100‑1 (MQ=255)
                                       cAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                              >  2:2446576/1‑100 (MQ=255)
                                       cAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTa                                                              >  2:1309675/1‑100 (MQ=255)
                                          tAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTaaaa                                                           <  2:1080694/100‑1 (MQ=255)
                                          tAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTaaaa                                                           <  1:979386/100‑1 (MQ=255)
                                            atatAAAATATTGAATTGTTATTTTTTTGTTTTGTTTAAAATAAAAAGCTGTATTTATTATTACGCCTTTAATTAATAAAAGGCAGGCTGTATTAAAAt                                                          <  1:1563410/99‑1 (MQ=255)
                                            atatAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAAtt                                                         >  2:1551012/1‑100 (MQ=255)
                                            atatAAAATATTGAATTGTTATTTTTGTGTGTTGCTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAATGAAATAATGTATTaaaact                                                         >  1:3171001/1‑98 (MQ=255)
                                               taAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAat                                                      <  2:2626091/100‑1 (MQ=255)
                                               taAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAat                                                      <  2:1669943/100‑1 (MQ=255)
                                               taAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAat                                                      <  1:2942348/100‑1 (MQ=255)
                                                aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGTCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                     <  1:296554/100‑1 (MQ=255)
                                                aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                     >  2:1339915/1‑100 (MQ=255)
                                                aaaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                     >  2:2054493/1‑100 (MQ=255)
                                                 aaaTATTAAATTTTTATTTTTGTGTTGTGTTTAAAATAAAAAGCCGTATTTAATATAACGGATTTAATTAATAAAAGGCAGTCTGTATTAAAATTAatat                                                    <  1:1713829/100‑1 (MQ=255)
                                                  aaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATAtt                                                   >  1:1070043/1‑100 (MQ=255)
                                                    tatTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCa                                                 <  1:2100905/100‑1 (MQ=255)
                                                       tGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                              <  2:2062985/100‑1 (MQ=255)
                                                        gAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTaaagatcggaag                                             >  1:2054081/1‑90 (MQ=255)
                                                         aaTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                              <  2:3429021/98‑1 (MQ=255)
                                                         aaTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGc                                             <  1:2794293/99‑1 (MQ=255)
                                                         aaTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCa                                            <  2:765839/100‑1 (MQ=255)
                                                             gTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATaaa                                        >  1:899938/1‑100 (MQ=255)
                                                                 tttttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGa                                    >  2:2478376/1‑100 (MQ=255)
                                                                 tttttGAGTGTTGGTTAAGATAAAAAGGCGTATTTATTATTACGGCTTGAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAACGa                                    <  1:1850343/100‑1 (MQ=255)
                                                                  ttttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGAt                                   <  2:1238362/100‑1 (MQ=255)
                                                                  ttttGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGAt                                   <  1:2416018/100‑1 (MQ=255)
                                                                      gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAg                                 <  1:2397361/98‑1 (MQ=255)
                                                                      gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                               <  1:1309667/100‑1 (MQ=255)
                                                                      gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                               <  2:3814963/100‑1 (MQ=255)
                                                                      gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                               <  1:85930/100‑1 (MQ=255)
                                                                      gtgtgtTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                               <  2:398689/100‑1 (MQ=255)
                                                                          gttgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATa                           >  2:2828960/1‑100 (MQ=255)
                                                                           ttgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAc                          <  2:2221638/100‑1 (MQ=255)
                                                                           ttgttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAc                          >  1:1547291/1‑100 (MQ=255)
                                                                             gttTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCg                        >  1:3778573/1‑100 (MQ=255)
                                                                              tttattaattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCgg                       >  2:703858/11‑100 (MQ=25)
                                                                              tttAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCgg                       >  1:2237979/1‑100 (MQ=255)
                                                                                tAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                     >  1:3298580/1‑100 (MQ=255)
                                                                                 aaGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGttt                    >  2:512413/1‑100 (MQ=255)
                                                                                 aaGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGttt                    >  1:3015393/1‑100 (MQ=255)
                                                                                 aaGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGttt                    >  2:2754069/1‑100 (MQ=255)
                                                                                  aGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTc                   <  2:1823829/100‑1 (MQ=255)
                                                                                   gATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCTAATATCGGTTTCa                  >  1:3471942/1‑100 (MQ=255)
                                                                                   gATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCa                  <  1:1245989/100‑1 (MQ=255)
                                                                                    aattaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCa                  >  1:1472886/5‑99 (MQ=37)
                                                                                    aTAAAAAGCCGTATTTATTATTCCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTTAAAGCATAAACCGATAGCCAATACCGGTTTCAg                 <  2:921042/100‑1 (MQ=255)
                                                                                     tAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                <  2:1426109/100‑1 (MQ=255)
                                                                                     tAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                >  2:3385213/1‑100 (MQ=255)
                                                                                     attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                >  2:3101077/4‑100 (MQ=37)
                                                                                     attaaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                >  1:7885/4‑100 (MQ=37)
                                                                                       aaaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa              >  2:1458491/1‑100 (MQ=255)
                                                                                       aaaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa              <  2:1289358/100‑1 (MQ=255)
                                                                                       aaaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa              <  1:2861600/100‑1 (MQ=255)
                                                                                        aaaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt             <  2:197224/100‑1 (MQ=255)
                                                                                        aaaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAt             >  2:2317528/1‑100 (MQ=37)
                                                                                         aaGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATg            <  2:3733476/100‑1 (MQ=255)
                                                                                         aaGCCGTAATAATAAATACGGCTTTAATTAATAAAAGGCAGACTGTATAAAAATTAATATTCAAAACATAAAACCATAACCAATAACGGTTTCAGCAATg            >  1:1315174/1‑100 (MQ=255)
                                                                                           gCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGaa          <  2:2258819/100‑1 (MQ=255)
                                                                                            ccGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGaaa         >  1:1384691/1‑100 (MQ=255)
                                                                                             cGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAAt        <  1:3584641/100‑1 (MQ=255)
                                                                                              gTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATg       <  1:259925/100‑1 (MQ=255)
                                                                                                 tttattATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTAAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCg    >  1:437336/1‑100 (MQ=255)
                                                                                                  ttattattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATGCCGGTTTCAGTAACGAAATAGc     >  1:3114054/1‑98 (MQ=255)
                                                                                                   tattattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTg  >  2:1011806/1‑100 (MQ=255)
                                                                                                   tattattaCGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTg  <  2:529024/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
ATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTG  >  NC_000913/720907‑721105

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: