Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,662,589 T→C D678D (GAT→GAC ynfF → putative selenate reductase YnfF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,662,5890TC98.8% 310.0 / ‑4.5 80D678D (GAT→GACynfFputative selenate reductase YnfF
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (34/45);  total (34/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

GTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGATCAACCCCATTGATGCT  >  NC_000913/1662454‑1662721
                                                                                                                                       |                                                                                                                                    
gTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACcccc                                                                                                                                  >  1:1387535/1‑140 (MQ=255)
 ttGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCt                                                                                                                                 >  1:1260567/1‑140 (MQ=255)
  tGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTg                                                                                                                                >  2:1022212/1‑140 (MQ=255)
     aaaCACCTTCGGGGAAAATCGAAATTTATTCTGAGCGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCgg                                                                                                                             <  1:1234000/140‑1 (MQ=255)
     aaaCACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCgg                                                                                                                             <  1:850134/140‑1 (MQ=255)
     aaaCACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCgg                                                                                                                             <  1:1428259/140‑1 (MQ=255)
       acacCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGaa                                                                                                                           <  2:756582/140‑1 (MQ=255)
       acacCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGaa                                                                                                                           <  1:467042/140‑1 (MQ=255)
          ccTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAAc                                                                                                                        <  2:310585/140‑1 (MQ=255)
          ccTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAAc                                                                                                                        <  2:57723/140‑1 (MQ=255)
              cGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTAt                                                                                                                    >  1:7559/1‑140 (MQ=255)
               ggggAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATc                                                                                                                   <  1:930688/140‑1 (MQ=255)
                    aaaTCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTg                                                                                                              <  2:889840/140‑1 (MQ=255)
                    aaaTCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTg                                                                                                              <  1:68640/140‑1 (MQ=255)
                          aaaTTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAACCCTATCCACTGCAGTTa                                                                                                        <  1:423486/140‑1 (MQ=255)
                           aaTTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTaa                                                                                                       <  1:5512/140‑1 (MQ=255)
                              ttATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACgg                                                                                                    <  2:1413190/140‑1 (MQ=255)
                                aTTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGc                                                                                                  <  1:498219/140‑1 (MQ=255)
                                aTTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGc                                                                                                  <  1:464654/140‑1 (MQ=255)
                                  tCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCtt                                                                                                >  2:518054/1‑140 (MQ=255)
                                           cTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACaa                                                                                       <  2:1387535/140‑1 (MQ=255)
                                                 aaGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGcgcg                                                                                 <  2:39491/140‑1 (MQ=255)
                                                   gATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTa                                                                               >  2:781658/1‑140 (MQ=255)
                                                     ttGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTAcc                                                                             >  1:1130475/1‑140 (MQ=255)
                                                     ttGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTAcc                                                                             >  1:1352585/1‑140 (MQ=255)
                                                      tGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTAccc                                                                            <  2:296108/140‑1 (MQ=255)
                                                      tGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTAccc                                                                            >  2:1269550/1‑140 (MQ=255)
                                                             aCCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTAc                                                                              >  2:750019/1‑131 (MQ=255)
                                                             aCCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTAc                                                                              <  1:750019/131‑1 (MQ=255)
                                                              ccTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGc                                                                    >  2:607696/1‑140 (MQ=255)
                                                              ccTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGc                                                                    >  1:978479/1‑140 (MQ=255)
                                                              ccTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGc                                                                    <  1:781658/140‑1 (MQ=255)
                                                              ccTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGc                                                                    <  1:607696/140‑1 (MQ=255)
                                                                  ggAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACg                                                                >  1:850223/1‑140 (MQ=255)
                                                                   gAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACgg                                                               >  2:1261549/1‑140 (MQ=255)
                                                                   gAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACgg                                                               <  1:917923/140‑1 (MQ=255)
                                                                   gAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACgg                                                               <  2:880911/140‑1 (MQ=255)
                                                                     aTTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCa                                                             <  1:667076/140‑1 (MQ=255)
                                                                      ttGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCaa                                                            >  2:1232154/1‑140 (MQ=255)
                                                                       tGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACgg                                                               >  2:1441154/1‑136 (MQ=255)
                                                                       tGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACgg                                                               <  1:1441154/136‑1 (MQ=255)
                                                                       tGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAat                                                           >  1:712424/1‑140 (MQ=255)
                                                                       tGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAat                                                           >  1:86585/1‑140 (MQ=255)
                                                                       tGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAat                                                           >  2:1158209/1‑140 (MQ=255)
                                                                        gAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAata                                                          <  2:978479/140‑1 (MQ=255)
                                                                        gAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAata                                                          >  1:288875/1‑140 (MQ=255)
                                                                         aaaaaaGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAatat                                                         <  2:1260567/140‑1 (MQ=255)
                                                                         aaaaaaGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAatat                                                         <  1:518054/140‑1 (MQ=255)
                                                                           aaaaGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTg                                                       <  2:7559/140‑1 (MQ=255)
                                                                            aaaGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGa                                                      <  2:600446/140‑1 (MQ=255)
                                                                            aaaGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGa                                                      <  2:288875/140‑1 (MQ=255)
                                                                                 tGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATgtgt                                                 <  1:505205/140‑1 (MQ=255)
                                                                                 tGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATgtgt                                                 <  2:458640/140‑1 (MQ=255)
                                                                                 tGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATgtgt                                                 <  2:86585/140‑1 (MQ=255)
                                                                                      ttattCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTAcag                                            <  2:1352585/140‑1 (MQ=255)
                                                                                          tCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTg                                                       >  1:1236814/1‑125 (MQ=255)
                                                                                          tCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTg                                                       <  2:1236814/125‑1 (MQ=255)
                                                                                           cATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGc                                       <  2:23552/140‑1 (MQ=255)
                                                                                           cATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGc                                       <  1:1232154/140‑1 (MQ=255)
                                                                                           cATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGc                                       <  2:166187/140‑1 (MQ=255)
                                                                                                   tCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCcac                               >  1:1556587/1‑140 (MQ=255)
                                                                                                    ccTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCcaca                              >  2:320258/1‑140 (MQ=255)
                                                                                                    ccTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGCAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCcaca                              >  2:1511000/1‑140 (MQ=255)
                                                                                                     cTGCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACaa                             <  1:1022212/140‑1 (MQ=255)
                                                                                                       gCGTATACCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAaga                           >  1:1339682/1‑140 (MQ=255)
                                                                                                          tataCCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGgt                        >  2:317840/1‑140 (MQ=255)
                                                                                                           ataCCCCAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCGTTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGgtg                       >  2:1088119/1‑140 (MQ=255)
                                                                                                                ccAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGAt                  <  1:1158209/140‑1 (MQ=255)
                                                                                                                 cAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGATc                 <  1:1224569/140‑1 (MQ=255)
                                                                                                                 cAGGTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGATc                 <  2:1339682/140‑1 (MQ=255)
                                                                                                                   ggTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGATCaa               >  1:524958/1‑140 (MQ=255)
                                                                                                                    gTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTATAGCAGGCCTGCCCACAAGAGGTGTGGATCAAc              >  1:406081/1‑140 (MQ=255)
                                                                                                                    gTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACaa                             <  1:670986/125‑1 (MQ=255)
                                                                                                                    gTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACaa                             >  2:670986/1‑125 (MQ=255)
                                                                                                                    gTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGATc                 <  1:227469/137‑1 (MQ=255)
                                                                                                                    gTTTTGATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGATc                 >  2:227469/1‑137 (MQ=255)
                                                                                                                         gATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGATCAACCCCAt         >  1:1330090/1‑140 (MQ=255)
                                                                                                                         gATGGCTGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGATCAACCCCAt         >  1:1477249/1‑140 (MQ=255)
                                                                                                                               tGGGACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGATCAACCCCATTGATGc   >  1:303811/1‑140 (MQ=255)
                                                                                                                                gggACGACCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCGTTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGATCAACCCCATTGATGCt  <  1:1088119/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                    
GTTGAAAACACCTTCGGGGAAAATCGAAATTTATTCTGAACGACTGGCGAAGATTGCAGATACCTGGGAATTGAAAAAAGATGAAATTATTCATCCCCTTCCTGCGTATACCCCAGGTTTTGATGGCTGGGACGATCCCCTGCGGAAAACCTATCCACTGCAGTTAACGGGCTTCCATTACAAAGCGCGTACCCACTCCAGCTACGGCAATATTGATGTGTTACAGCAGGCCTGCCCACAAGAGGTGTGGATCAACCCCATTGATGCT  >  NC_000913/1662454‑1662721

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: