Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,662,532 (A)6→5 coding (1977/2424 nt) ynfF → putative selenate reductase YnfF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,662,5270A.99.6% 1346.3 / ‑3.0 257coding (1972/2424 nt)ynfFputative selenate reductase YnfF
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (126/130);  minor base T (0/1);  total (126/131)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 295 total aligned reads displayed.

ATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGG  >  NC_000913/1662427‑1662623
                                                                                                    |                                                                                                
aTTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGa                                                                                                  <  1:2531881/101‑1 (MQ=255)
 ttCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGaa                                                                                                 >  1:2532227/1‑101 (MQ=255)
  tCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGaaa                                                                                                <  1:2531695/101‑1 (MQ=255)
  tCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGaaa                                                                                                <  2:2531863/101‑1 (MQ=255)
   ccGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGaaaa                                                                                               <  2:2531759/101‑1 (MQ=255)
   ccGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGaaaa                                                                                               >  1:2532231/1‑101 (MQ=255)
   ccGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGaaaa                                                                                               <  2:2531656/101‑1 (MQ=255)
   ccGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGaaaa                                                                                               >  2:2532229/1‑101 (MQ=255)
   ccGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGaaaa                                                                                               <  1:2531561/101‑1 (MQ=255)
     gTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGaaaaag                                                                                             <  2:2531630/101‑2 (MQ=255)
      tGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGaaaaaga                                                                                            >  2:2532234/1‑99 (MQ=255)
       gAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGaaaaagat                                                                                           <  2:2531744/101‑4 (MQ=255)
        aaGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGaaaaagatg                                                                                          <  2:2531502/101‑5 (MQ=255)
        aaGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGaaaaagatg                                                                                          >  1:2532238/1‑97 (MQ=255)
           aTCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGaaa                                                                                      >  2:2532241/1‑101 (MQ=255)
             ccACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAatt                                                                                    <  2:2531836/101‑1 (MQ=255)
              cacaGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAANGATGAAatta                                                                                   <  1:2531753/101‑1 (MQ=255)
              cacaGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAatta                                                                                   >  2:2532245/1‑101 (MQ=255)
                caGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAattatt                                                                                 >  1:2532247/1‑101 (MQ=255)
                  ggCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCa                                                                               <  1:2531705/101‑1 (MQ=255)
                    cAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATc                                                                             >  1:2532251/1‑101 (MQ=255)
                     aaaCCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATcc                                                                            >  2:2532253/1‑101 (MQ=255)
                      aaCCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATccc                                                                           <  2:2531813/101‑1 (MQ=255)
                       aCCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATcccc                                                                          >  1:2532256/1‑101 (MQ=255)
                        cccGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCt                                                                         <  2:2531845/101‑1 (MQ=255)
                         ccGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCtt                                                                        <  1:2531700/101‑1 (MQ=255)
                         ccGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCtt                                                                        >  1:2532261/1‑101 (MQ=255)
                         ccGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCtt                                                                        <  1:2531842/101‑1 (MQ=255)
                         ccGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCtt                                                                        >  2:2532260/1‑101 (MQ=255)
                          cGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTc                                                                       <  1:2532121/101‑1 (MQ=255)
                          cGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTc                                                                       <  2:2531790/101‑1 (MQ=255)
                           gTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCa                                                                               <  2:2531853/92‑1 (MQ=255)
                           gTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTcc                                                                      <  2:2531688/101‑1 (MQ=255)
                           gTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTcc                                                                      <  1:2531723/101‑1 (MQ=255)
                           gTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTcc                                                                      <  2:2531887/101‑1 (MQ=255)
                            ttGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCt                                                                     <  1:2531986/101‑1 (MQ=255)
                            ttGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCt                                                                     <  1:2531968/101‑1 (MQ=255)
                            ttGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCt                                                                     >  1:2532266/1‑101 (MQ=255)
                             tGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTg                                                                    <  2:2531879/101‑1 (MQ=255)
                              gAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGc                                                                   <  1:2532075/101‑1 (MQ=255)
                              gAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGc                                                                   <  2:2531907/101‑1 (MQ=255)
                              gAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGc                                                                   >  2:2532270/1‑101 (MQ=255)
                               aaaaCACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCg                                                                  <  1:2531771/101‑1 (MQ=255)
                               aaaaCACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCg                                                                  >  1:2532272/1‑101 (MQ=255)
                                aaaCACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGt                                                                 >  1:2532276/1‑101 (MQ=255)
                                aaaCACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGt                                                                 >  2:2532278/1‑101 (MQ=255)
                                  acacCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGtat                                                               <  1:2531924/101‑1 (MQ=255)
                                  acacCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGtat                                                               >  1:2532281/1‑101 (MQ=255)
                                       ttCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCg                                                                  >  1:2532285/1‑93 (MQ=255)
                                        tCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCa                                                         >  2:2532286/1‑101 (MQ=255)
                                        tCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCa                                                         <  1:2531810/101‑1 (MQ=255)
                                          ggggAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAgg                                                       >  2:2532289/1‑101 (MQ=255)
                                          ggggAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAgg                                                       <  1:2531829/101‑1 (MQ=255)
                                          ggggAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAgg                                                       <  1:2531822/101‑1 (MQ=255)
                                          ggggAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAgg                                                       <  2:2531959/101‑1 (MQ=255)
                                           gggAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGt                                                      <  2:2532167/101‑1 (MQ=255)
                                           gggAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGt                                                      <  1:2531936/101‑1 (MQ=255)
                                            ggAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGtt                                                     >  1:2532292/1‑101 (MQ=255)
                                             gAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGttt                                                    >  2:2532295/1‑101 (MQ=255)
                                             gAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTACACCCCAGGttt                                                    <  1:2531944/101‑1 (MQ=255)
                                               aaaTCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTg                                                  >  2:2532298/1‑101 (MQ=255)
                                                aaTCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGa                                                 >  1:2532301/1‑101 (MQ=255)
                                                 aTCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGAt                                                <  1:2531911/101‑1 (MQ=255)
                                                   cGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATgg                                              <  1:2531680/101‑1 (MQ=255)
                                                   cGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATgg                                              >  2:2532305/1‑101 (MQ=255)
                                                    gAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGc                                             >  2:2532308/1‑101 (MQ=255)
                                                      aaTTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTg                                           >  1:2532311/1‑101 (MQ=255)
                                                       aTTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTgg                                          <  2:2531895/101‑1 (MQ=255)
                                                        tttATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTggg                                         >  2:2532314/1‑101 (MQ=255)
                                                         ttATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGa                                        <  1:2531798/101‑1 (MQ=255)
                                                          tATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGAc                                       >  1:2532316/1‑101 (MQ=255)
                                                          tATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGAc                                       <  1:2532032/101‑1 (MQ=255)
                                                            ttCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGa                                     >  2:2532318/1‑101 (MQ=255)
                                                                 aaCGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATcccc                                <  2:2531932/101‑1 (MQ=255)
                                                                 aaCGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATcccc                                >  1:2532323/1‑101 (MQ=255)
                                                                  aCGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCt                               >  2:2532325/1‑101 (MQ=255)
                                                                  aCGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCt                               <  2:2531557/101‑1 (MQ=255)
                                                                  aCGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCt                               <  1:2531775/101‑1 (MQ=255)
                                                                   cGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTg                              <  2:2532015/101‑1 (MQ=255)
                                                                     aCTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCg                            >  1:2532329/1‑101 (MQ=255)
                                                                       tGGCGAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGa                          >  2:2532331/1‑101 (MQ=255)
                                                                           gAAGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAAc                      >  1:2532334/1‑101 (MQ=255)
                                                                            aaGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAAcc                     >  2:2532338/1‑101 (MQ=255)
                                                                             aGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCt                    <  2:2531917/101‑1 (MQ=255)
                                                                             aGATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCt                    >  1:2532340/1‑101 (MQ=255)
                                                                              gATTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTa                   <  1:2531869/101‑1 (MQ=255)
                                                                               aTTGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTAt                  >  2:2532343/1‑101 (MQ=255)
                                                                                 tGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATcc                >  1:2532345/1‑101 (MQ=255)
                                                                                 tGCAGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATcc                >  2:2532349/1‑101 (MQ=255)
                                                                                    aGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACt             <  2:2532087/101‑1 (MQ=255)
                                                                                    aGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACt             >  1:2532351/1‑101 (MQ=255)
                                                                                    aGATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACt             <  1:2531905/101‑1 (MQ=255)
                                                                                     gATACCTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTg            <  2:2532005/101‑1 (MQ=255)
                                                                                         ccTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGt        >  1:2532356/1‑101 (MQ=255)
                                                                                          cTGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGtt       <  1:2531663/101‑1 (MQ=255)
                                                                                           tGGGAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTa      >  2:2532358/1‑101 (MQ=255)
                                                                                            gggAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTaa     >  1:2532363/1‑101 (MQ=255)
                                                                                            gggAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTaa     >  2:2532362/1‑101 (MQ=255)
                                                                                              gAATTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACg   <  1:2531963/101‑1 (MQ=255)
                                                                                               aaTTGAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACgg  <  1:2531888/101‑1 (MQ=255)
                                                                                                    |                                                                                                
ATTCCGTGAAGATCCACAGGCAAACCCGTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGG  >  NC_000913/1662427‑1662623

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