Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,423,570 A→C intergenic (‑31/+166) yfcI ← / ← hisP transposase_31 family protein/histidine ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,423,5700AC99.3% 492.1 / ‑4.1 142intergenic (‑31/+166)yfcI/hisPtransposase_31 family protein/histidine ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (62/79);  total (62/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

Only 100 of 142 total aligned reads displayed.

TGGATGGCGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGACGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTATTTTTAGCGGACTGTAACTTTTAC  >  NC_000913/2423477‑2423665
                                                                                             |                                                                                               
tGGATGGCGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGAc                                                                                           <  2:3485529/100‑1 (MQ=255)
 ggATGGCGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACa                                                                                          <  1:2955694/100‑1 (MQ=255)
 ggATGGCGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACa                                                                                          >  1:416609/1‑100 (MQ=255)
     ggCGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATaa                                                                                      >  2:749830/1‑100 (MQ=255)
     ggCGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATaa                                                                                      <  1:3683170/100‑1 (MQ=255)
     ggCGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATaa                                                                                      <  1:1584237/100‑1 (MQ=255)
       cGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAAcc                                                                                    >  2:491788/1‑100 (MQ=255)
        gTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCa                                                                                   >  1:26380/1‑100 (MQ=255)
        gTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCa                                                                                   >  1:3169601/1‑100 (MQ=255)
          aaaaaaGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACg                                                                                 <  2:1927336/100‑1 (MQ=255)
            aaaaGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGAt                                                                               >  1:2468285/1‑100 (MQ=255)
            aaaaGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGAt                                                                               >  1:2264050/1‑100 (MQ=255)
            aaaaGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGAt                                                                               <  2:2689213/100‑1 (MQ=255)
            aaaaGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGAt                                                                               >  1:1331963/1‑100 (MQ=255)
             aaaGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGAtt                                                                              <  2:84799/100‑1 (MQ=255)
             aaaGATTCAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCTACAATAACCACGAtt                                                                              <  2:3035450/100‑1 (MQ=255)
              aaGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTa                                                                             <  2:460986/100‑1 (MQ=255)
              aaGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTa                                                                             <  1:1166557/100‑1 (MQ=255)
                gATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCTAATTTAAGATGCCGCGACAATAACCACGATTAac                                                                           <  1:3785158/100‑1 (MQ=255)
                 aTTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACACTAACCACGATTAaca                                                                          >  2:1201952/1‑100 (MQ=255)
                 aTTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAaca                                                                          <  1:3269018/100‑1 (MQ=255)
                  tttAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAacac                                                                         <  1:723964/100‑1 (MQ=255)
                   ttAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACAcc                                                                        >  1:2580547/1‑100 (MQ=255)
                    tAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACAccc                                                                       >  2:822636/1‑100 (MQ=255)
                     aaaTACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACAcccc                                                                      >  2:1117343/1‑99 (MQ=255)
                     aaaTACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCa                                                                      >  2:3122674/1‑100 (MQ=255)
                     aaaTACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCa                                                                      >  1:1749006/1‑100 (MQ=255)
                     aaaTACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCa                                                                      <  1:212536/100‑1 (MQ=255)
                      aaTACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCaa                                                                     <  2:2514753/100‑1 (MQ=255)
                      aaTACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCaa                                                                     >  2:228404/1‑100 (MQ=255)
                       aTACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAAt                                                                    >  1:3714292/1‑100 (MQ=255)
                         aCCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATgc                                                                  <  2:2729965/100‑1 (MQ=255)
                         aCCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATgc                                                                  <  2:2306045/100‑1 (MQ=255)
                         aCCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGACGCGACAATAACCACGATTAACACCCAATgc                                                                  <  1:1084731/100‑1 (MQ=255)
                          ccGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATgca                                                                 >  2:2510198/1‑100 (MQ=255)
                          ccGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATgca                                                                 <  2:1628431/100‑1 (MQ=255)
                           cGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACCATAACCACGCTTAACACCCAATgccg                                                                >  2:2100331/1‑98 (MQ=255)
                           cGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATgcag                                                                >  1:1236804/1‑100 (MQ=255)
                           cGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATgcag                                                                <  1:2313918/100‑1 (MQ=255)
                            gCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATgcagc                                                               <  2:2236912/100‑1 (MQ=255)
                            gCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATgcagc                                                               <  1:1913266/100‑1 (MQ=255)
                             catcatGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATgcagca                                                              <  1:2017261/100‑1 (MQ=255)
                               tcatGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACACAATGCCgcagg                                                            <  2:692205/100‑1 (MQ=255)
                                catGCGGCGTGGAAGTTGTCGATAGTGACATGGCGTTCCGTCACCCCCCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAgcaggc                                                           >  2:1608823/1‑100 (MQ=255)
                                   gcggcgCGGAAGTTTTTGATATTTGAAGGGCGTGCCGTCACCCCGTAAATATAATACCCGGCGTCAACAACCACGATTAACACCCAATGCAgccgccagg                                                        <  2:3757795/100‑2 (MQ=255)
                                     ggcgTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGttt                                                      <  1:1030741/100‑1 (MQ=255)
                                       cgTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGttttt                                                    <  2:2043162/100‑1 (MQ=255)
                                       cgTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGttttt                                                    >  2:1688779/1‑100 (MQ=255)
                                        gTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTa                                                   <  2:3891850/100‑1 (MQ=255)
                                        gTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTa                                                   <  1:595935/100‑1 (MQ=255)
                                          ggAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATc                                                 <  2:3056972/100‑1 (MQ=255)
                                          ggAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATc                                                 >  1:2531975/1‑100 (MQ=255)
                                            aaGTTGTCGATATTGTCATGGCGTTCCGTCACCACTCAAATTTAAGATGCCGCGACACTAACAACGATTAACACCCAATGCAGCAGGCAGTCTTTATCtt                                               >  2:1774307/1‑100 (MQ=255)
                                             aGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCttt                                              >  2:3358549/1‑100 (MQ=255)
                                              gTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTg                                             <  2:3054317/100‑1 (MQ=255)
                                               ttGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTgg                                            <  2:173641/100‑1 (MQ=255)
                                                 gTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGAc                                          >  1:4002183/1‑100 (MQ=255)
                                                 gTCGATATTGTCATGGCGTTCCGTCACACCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGAc                                          <  1:2163583/100‑1 (MQ=255)
                                                   cGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAg                                        >  2:3497786/1‑100 (MQ=255)
                                                    gATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGc                                       <  1:3195358/100‑1 (MQ=255)
                                                         tGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAAtt                                  >  1:2307494/1‑100 (MQ=255)
                                                          ttCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTc                                 <  1:556073/99‑1 (MQ=255)
                                                           tCATGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCa                                >  2:544885/1‑100 (MQ=255)
                                                           tCAGGGCGTACCGTCACCCCTCAACGTTAGGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCa                                <  2:1844390/100‑1 (MQ=255)
                                                             aTGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGa                              >  1:3908320/1‑100 (MQ=255)
                                                             aTGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGa                              <  2:758016/100‑1 (MQ=255)
                                                             aTGGCGTTCCGTCACCCCTCAAATTGAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGa                              <  2:2815853/100‑1 (MQ=255)
                                                              tGGCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGaa                             <  2:3586839/100‑1 (MQ=255)
                                                               ggCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAg                            <  1:2285372/100‑1 (MQ=255)
                                                               ggCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAg                            <  1:46096/100‑1 (MQ=255)
                                                               ggCGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAg                            <  1:1372973/100‑1 (MQ=255)
                                                               ggCGTTCCGTCACCCCTCAAATTTAAGATGCCGCAACAATAACCACGATTAACCCCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAg                            <  2:2036030/100‑1 (MQ=255)
                                                                 cGTTCCGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTa                          >  2:567529/1‑100 (MQ=255)
                                                                     ccGTCACCCCTCTAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGACTTCAGAAGTAtgtt                      >  1:1839817/1‑100 (MQ=255)
                                                                      cGTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTAttttt                     >  2:3611512/1‑100 (MQ=255)
                                                                       gTCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTATTTTTa                    >  1:1524979/1‑100 (MQ=255)
                                                                        tCACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTATTTTTAg                   <  1:334076/100‑1 (MQ=255)
                                                                         cACCCCTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTATTTTTAGc                  <  2:1763402/100‑1 (MQ=255)
                                                                          aCCGCTCAAATTTAAGATGCGGCGACAATAACCAGCATTAACACCCAATGCAGCAGGCAGTTTTTATCTCGAGACAGCGAATTCAGAAGTATTTTTcgcg                 <  2:2479472/100‑1 (MQ=255)
                                                                           ccccTGAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCGGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTATTTTTAGCgg                <  1:3869119/100‑1 (MQ=255)
                                                                            cccTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTATTTTTAGCGGa               >  1:3596236/1‑100 (MQ=255)
                                                                            cccTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTATTTTTAGCGGa               >  1:2230465/1‑100 (MQ=255)
                                                                            cacTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAAAGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTATTTGTAGCGGa               >  2:3711876/3‑100 (MQ=255)
                                                                             ccTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACTCAATTGAGCAGGCAGTTTATATCTTTGGACATCGAATTCAGAAGTATCTTTAGCTGAc              >  2:3406289/1‑100 (MQ=255)
                                                                             ccTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTATTTTTAGCGGAc              <  2:3970720/100‑1 (MQ=255)
                                                                             ccTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTATTTTTAGCGGAc              <  1:242873/100‑1 (MQ=255)
                                                                             ccTCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTATTTTTAGCGGAc              <  2:2292980/100‑1 (MQ=255)
                                                                               tCAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTATTTTTAGCGGACTg            <  2:2975187/100‑1 (MQ=255)
                                                                                cAAATTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTGTCTATCTTTGGACAGCGAATTCAGAAGTATTTTTAGCGGACTGt           >  1:1590585/1‑100 (MQ=255)
                                                                                 aaaTTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTATTTTTAGCGGACTGTa          <  1:2567163/100‑1 (MQ=255)
                                                                                  aaTTTAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTATTTTTAGCGGACTGTaa         >  2:3982182/1‑100 (MQ=255)
                                                                                    tttAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTATTTTTAGCGGACTGTAACt       <  1:965250/100‑1 (MQ=255)
                                                                                     ttAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCCGTCTTTATCTTTGGACAGCGAATTCAGAAGTATTTTTAGCGGACTGTAACtt      >  1:3452759/1‑100 (MQ=255)
                                                                                     ttAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTATTTTTAGCGGACTGTAACtt      >  1:1657146/1‑100 (MQ=255)
                                                                                      tAAGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTATTTTTAGCGGACTGTAACttt     <  1:16101/100‑1 (MQ=255)
                                                                                       aaGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTATTTTTAGCGGACTGTAACtttt    <  1:2632117/100‑1 (MQ=255)
                                                                                        aGATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTATTTTTAGCGGACTGTAACTTTTa   <  1:1491882/100‑1 (MQ=255)
                                                                                         gATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTATTTTTAGCGGACTGTAACTTTTAc  <  2:3374472/100‑1 (MQ=255)
                                                                                         gATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTATTTTTAGCGGACTGTAACTTTTAc  <  2:3254959/100‑1 (MQ=255)
                                                                                         gATGCCGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTATTTTTAGCGGACTGTAACTTTTAc  >  1:17837/1‑100 (MQ=255)
                                                                                             |                                                                                               
TGGATGGCGTAAAAAAGATTTAAATACCGCATCATGCGGCGTGGAAGTTGTCGATATTGTCATGGCGTTCCGTCACCCCTCAAATTTAAGATGACGCGACAATAACCACGATTAACACCCAATGCAGCAGGCAGTTTTTATCTTTGGACAGCGAATTCAGAAGTATTTTTAGCGGACTGTAACTTTTAC  >  NC_000913/2423477‑2423665

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: