Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,285,193 T→A 6.1% W19R (TGG→AGG)  narJ → molybdenum‑cofactor‑assembly chaperone delta subunit of nitrate reductase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,285,1930TA 6.1% 357.9 / 10.5 99W19R (TGG→AGG) narJmolybdenum‑cofactor‑assembly chaperone delta subunit of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base T (45/48);  new base A (3/3);  total (48/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-01

Only 100 of 102 total aligned reads displayed.

GTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCG  >  NC_000913/1285045‑1285337
                                                                                                                                                    |                                                                                                                                                
gTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTAt                                                                                                                                                  >  1:13858/1‑149 (MQ=255)
 tGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATg                                                                                                                                                 <  1:574082/149‑1 (MQ=255)
 tGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATg                                                                                                                                                 <  2:382296/149‑1 (MQ=255)
    tGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGgca                                                                                                                                              <  1:169375/149‑1 (MQ=255)
    tGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGgca                                                                                                                                              <  2:662693/149‑1 (MQ=255)
        tGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATgg                                                                                                                                                <  2:1823634/143‑1 (MQ=255)
        tGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATgg                                                                                                                                                >  1:1823634/1‑143 (MQ=255)
         gCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATc                                                                                                                                         <  4:1595212/149‑1 (MQ=255)
           cACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATcaa                                                                                                                                       <  3:963006/149‑1 (MQ=255)
               gTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAag                                                                                                                                   >  2:714236/1‑149 (MQ=255)
                ttCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAaga                                                                                                                                  <  2:982689/149‑1 (MQ=255)
                ttCAGATACCAAATTCAATCTCTTCAACAGCCGTCATATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATAAACAaga                                                                                                                                  >  2:1110478/1‑149 (MQ=255)
                   aGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGAt                                                                                                                               <  3:1763237/149‑1 (MQ=255)
                                aaTCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATgg                                                                                                                                                <  1:1008441/119‑1 (MQ=255)
                                aaTCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATgg                                                                                                                                                >  2:1008441/1‑119 (MQ=255)
                                 aTCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATcaacaa                                                                                                                                    <  3:1577266/130‑1 (MQ=255)
                                 aTCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATcaacaa                                                                                                                                    >  4:1577266/1‑130 (MQ=255)
                                  tCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGc                                                                                                                >  1:1045554/1‑149 (MQ=255)
                                  tCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGc                                                                                                                >  3:1357041/1‑149 (MQ=255)
                                  tCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGc                                                                                                                >  1:1388911/1‑149 (MQ=255)
                                   cTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGa                                                                                                                    >  1:696943/1‑144 (MQ=255)
                                   cTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGa                                                                                                                    <  2:696943/144‑1 (MQ=255)
                                     gTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCATAAGGCAGCAt                                                                                                                                          >  4:1329133/1‑120 (MQ=255)
                                     gTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCATAAGGCAGCAt                                                                                                                                          <  3:1329133/120‑1 (MQ=255)
                                       tCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTAt                                                                                                                                                  >  3:502487/1‑110 (MQ=255)
                                       tCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTAt                                                                                                                                                  <  4:502487/110‑1 (MQ=255)
                                        cAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTc                                                                                                          <  1:714236/149‑1 (MQ=255)
                                         aaCAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCg                                                                                                         >  4:711900/1‑149 (MQ=255)
                                                   gtATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGc                                                                                               <  2:1388911/149‑1 (MQ=255)
                                                        gATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCaaaa                                                                                          >  4:1172746/1‑149 (MQ=255)
                                                           gCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAagag                                                                                       <  1:1110478/149‑1 (MQ=255)
                                                             cATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGa                                                                                     <  3:711900/149‑1 (MQ=255)
                                                              aTCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGAt                                                                                    >  2:162381/1‑149 (MQ=255)
                                                              aTCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATATGCCAAAAGAGGAt                                                                                    >  2:42999/1‑149 (MQ=255)
                                                                       aCCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATgcg                                                                           >  3:534446/1‑149 (MQ=255)
                                                                         cAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAaga                                                                                                                                  <  3:1475290/92‑1 (MQ=255)
                                                                         cAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAaga                                                                                                                                  >  4:1475290/1‑92 (MQ=255)
                                                                         cAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCt                                                                         >  4:1701439/1‑149 (MQ=255)
                                                                          aGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTg                                                                        <  4:534446/149‑1 (MQ=255)
                                                                           gCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATATGCCAAAAGAGGATGCCCATGCGCTgg                                                                       <  1:42999/149‑1 (MQ=255)
                                                                             aaaaCGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGt                                                                                                           <  2:559329/111‑1 (MQ=255)
                                                                             aaaaCGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGt                                                                                                           >  1:559329/1‑111 (MQ=255)
                                                                             aaaaCGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGc                                                                     <  2:13858/149‑1 (MQ=255)
                                                                              aaaCGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAaga                                                                                        >  3:128722/1‑129 (MQ=255)
                                                                              aaaCGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAaga                                                                                        <  4:128722/129‑1 (MQ=255)
                                                                                   gAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTc                                                               >  3:700057/1‑149 (MQ=255)
                                                                                     gCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCt                                                             <  4:700057/149‑1 (MQ=255)
                                                                                      ccGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATcaacaa                                                                                                                                    <  2:981263/77‑1 (MQ=255)
                                                                                      ccGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATcaacaa                                                                                                                                    >  1:981263/1‑77 (MQ=255)
                                                                                      ccGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTcc                                                              <  4:880092/147‑1 (MQ=255)
                                                                                      ccGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTcc                                                              >  3:880092/1‑147 (MQ=255)
                                                                                              aTGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGc                                                                                               >  3:1484234/1‑106 (MQ=255)
                                                                                              aTGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGc                                                                                               <  4:1484234/106‑1 (MQ=255)
                                                                                              aTGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGcc                                                                                 >  4:1248940/1‑120 (MQ=255)
                                                                                              aTGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGcc                                                                                 <  3:1248940/120‑1 (MQ=255)
                                                                                               tGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAaga                                                                                        >  1:113992/1‑112 (MQ=255)
                                                                                               tGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAaga                                                                                        <  2:113992/112‑1 (MQ=255)
                                                                                                  tCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAacg                                                >  2:483115/1‑149 (MQ=255)
                                                                                                   cGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAacga                                               <  1:162381/149‑1 (MQ=255)
                                                                                                    gAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTAGGCATTTTCCTGCGCGATTTAacgac                                              >  3:1869517/1‑149 (MQ=255)
                                                                                                     aaCTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAacgacg                                             >  3:1029087/1‑149 (MQ=255)
                                                                                                     aaCTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAacgacg                                             >  1:964164/1‑149 (MQ=255)
                                                                                                         cGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATgg                                         <  2:964164/149‑1 (MQ=255)
                                                                                                          gTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGa                                        <  3:1172746/149‑1 (MQ=255)
                                                                                                          gTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGa                                        >  3:846483/1‑149 (MQ=255)
                                                                                                           tGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGAt                                       <  3:1701439/149‑1 (MQ=255)
                                                                                                            gATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATc                                      <  4:846483/149‑1 (MQ=255)
                                                                                                               tGTATCGCGTCTCCTTGAATATCCGGGTGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGc                                   <  4:123616/149‑1 (MQ=255)
                                                                                                                 tATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTc                                 <  4:1029087/149‑1 (MQ=255)
                                                                                                                        tctcCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGccc                          <  1:483115/149‑1 (MQ=255)
                                                                                                                          tcCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAg                        >  4:639498/1‑149 (MQ=255)
                                                                                                                             ttGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGc                                   <  1:1473172/135‑1 (MQ=255)
                                                                                                                             ttGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGc                                   >  2:1473172/1‑135 (MQ=255)
                                                                                                                                  tatCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTaa                                                  >  3:544470/1‑115 (MQ=255)
                                                                                                                                  tatCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTaa                                                  <  4:544470/115‑1 (MQ=255)
                                                                                                                                   atCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTgg                                                                       >  3:1334223/1‑93 (MQ=255)
                                                                                                                                   atCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTgg                                                                       <  4:1334223/93‑1 (MQ=255)
                                                                                                                                    tCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACa              <  2:1045554/149‑1 (MQ=255)
                                                                                                                                    tCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACa              <  4:1357041/149‑1 (MQ=255)
                                                                                                                                      cGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTgg                                                                       >  2:656976/1‑90 (MQ=255)
                                                                                                                                      cGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTgg                                                                       <  1:656976/90‑1 (MQ=255)
                                                                                                                                       ggATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCg           >  1:526321/1‑149 (MQ=255)
                                                                                                                                         aTGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGaa         <  2:526321/149‑1 (MQ=255)
                                                                                                                                          tgctgcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAAc        <  3:639498/149‑1 (MQ=255)
                                                                                                                                          tgctgcCATAATGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTc                                                               >  4:171522/1‑94 (MQ=255)
                                                                                                                                          tgctgcCATAAGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTc                                                               <  3:171522/94‑1 (MQ=255)
                                                                                                                                          tgctgcCATAAGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGaa         >  4:981591/1‑148 (MQ=255)
                                                                                                                                          tgctgcCATAAGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGaa         <  3:981591/148‑1 (MQ=255)
                                                                                                                                             tgcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTgg                                                                       >  3:1688897/1‑83 (MQ=255)
                                                                                                                                             tgcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTgg                                                                       <  4:1688897/83‑1 (MQ=255)
                                                                                                                                             tgcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTaa                                                  <  1:484980/104‑1 (MQ=255)
                                                                                                                                             tgcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTaa                                                  >  2:484980/1‑104 (MQ=255)
                                                                                                                                             tgcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAAcgacga                                            <  3:258359/110‑1 (MQ=255)
                                                                                                                                             tgcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAAcgacga                                            >  4:258359/1‑110 (MQ=255)
                                                                                                                                             tgcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGAt                                       >  3:996697/1‑115 (MQ=255)
                                                                                                                                             tgcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGt     >  2:947383/1‑149 (MQ=255)
                                                                                                                                             tgcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGt     >  1:1371676/1‑149 (MQ=255)
                                                                                                                                             tgcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGt     >  1:171406/1‑149 (MQ=255)
                                                                                                                                             tgcCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGt     >  2:801759/1‑149 (MQ=255)
                                                                                                                                                cTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCg  >  1:1581482/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GTGATGGCTGCCACGGTTCAGATACCAAATTCAATCTGTTCAACAGCCGTCGTATCGATGCCATCGATGTGACCAGCAAAACGGAGCCGCATCCATGATCGAACTCGTGATTGTATCGCGTCTCCTTGAATATCCGGATGCTGCCTTATGGCAGCATCAACAAGAGATGTTTGAGGCGATTGCCGCGTCGAAAAATCTGCCAAAAGAGGATGCCCATGCGCTGGGCATTTTCCTGCGCGATTTAACGACGATGGATCCGCTCGATGCCCAGGCGCAGTACAGCGAACTGTTCG  >  NC_000913/1285045‑1285337

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