Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,289,094 T→A 6.9% I444N (ATT→AAT)  nrfA → nitrite reductase, formate‑dependent, cytochrome

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,289,0940TA 6.9% 449.3 / 20.1 131I444N (ATT→AAT) nrfAnitrite reductase, formate‑dependent, cytochrome
Reads supporting (aligned to +/- strand):  ref base T (57/65);  new base A (5/4);  total (62/69)
Fisher's exact test for biased strand distribution p-value = 7.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

Only 100 of 132 total aligned reads displayed.

GAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATGGAGTGA  >  NC_000913/4288949‑4289242
                                                                                                                                                 |                                                                                                                                                    
gaagaaGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTgg                                                                                                                                                   >  4:682796/1‑149 (MQ=255)
 aagaagGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGc                                                                                                                                                  <  4:831829/149‑1 (MQ=255)
  agaagGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGcc                                                                                                                                                 <  2:542510/149‑1 (MQ=255)
   gaagGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGACCATTGGCCt                                                                                                                                                >  4:649972/1‑149 (MQ=255)
     agGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGa                                                                                                                                              >  2:763464/1‑149 (MQ=255)
      ggTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGaa                                                                                                                                             >  2:662458/1‑149 (MQ=255)
       gTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGGGCCCGGATGTCTCAACCAAAGAGAAAGCCCATCAGGCCATTGGCCAGAAc                                                                                                                                            >  2:466014/1‑149 (MQ=255)
       gTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAAc                                                                                                                                            >  1:840415/1‑149 (MQ=255)
        tttACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCAATGGCCTgctgg                                                                                                                                           <  3:230947/149‑5 (MQ=255)
             ggATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGaa                                                                                                                                      <  4:422551/149‑1 (MQ=255)
              gATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAAc                                                                                                                                     >  4:153239/1‑149 (MQ=255)
              gATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAAc                                                                                                                                     >  1:681878/1‑149 (MQ=255)
                 gCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACaaa                                                                                                                                  <  3:229767/149‑1 (MQ=255)
                   tCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATc                                                                                                                                >  4:714074/1‑149 (MQ=255)
                      gTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAg                                                                                                                             <  4:125388/149‑1 (MQ=255)
                      gTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAg                                                                                                                             >  1:634311/1‑149 (MQ=255)
                       tACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAgg                                                                                                                            <  2:3002/149‑1 (MQ=255)
                       tACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAgg                                                                                                                            <  1:846365/149‑1 (MQ=255)
                        aCGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGc                                                                                                                           <  4:293982/149‑1 (MQ=255)
                             gATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATTACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCgaga                                                                                                                      >  3:334547/1‑149 (MQ=255)
                              aTGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGaa                                                                                                                     <  1:750106/149‑1 (MQ=255)
                                 gATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCa                                                                                                                  <  1:410967/149‑1 (MQ=255)
                                  aTAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAg                                                                                                                 <  1:763464/149‑1 (MQ=255)
                                  aTAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAg                                                                                                                 <  1:560558/149‑1 (MQ=255)
                                    aaaGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGa                                                                                                               <  2:789435/149‑1 (MQ=255)
                                    aaaGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGa                                                                                                               <  3:780590/149‑1 (MQ=255)
                                    aaaGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGa                                                                                                               <  2:181426/149‑1 (MQ=255)
                                      aGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACt                                                                                                             >  4:556564/1‑149 (MQ=255)
                                       gcggcgGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAg                                                                                                                 >  3:238276/1‑144 (MQ=255)
                                        cggcggATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTtc                                                                                                           <  4:588636/149‑1 (MQ=255)
                                         ggcggATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTtca                                                                                                          <  1:588618/149‑1 (MQ=255)
                                         ggcggATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTtca                                                                                                          >  2:26670/1‑149 (MQ=255)
                                         ggcggATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTtca                                                                                                          <  4:818627/149‑1 (MQ=255)
                                          gcggATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTtcat                                                                                                         <  4:98581/149‑1 (MQ=255)
                                                  acgcacCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGTCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGAc                                                                                                              >  2:678/1‑136 (MQ=255)
                                                  acgcacCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACgg                                                                                                 >  1:89173/1‑149 (MQ=255)
                                                    gcacCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCAATGGCCTgctgg                                                                                                                                           <  3:675718/105‑5 (MQ=255)
                                                      acCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCAATGGc                                                                                                                                                  >  2:222974/1‑96 (MQ=255)
                                                      acCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCAATGGc                                                                                                                                                  <  1:222974/96‑1 (MQ=255)
                                                         aaaCTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACAGTGATccc                                                                                          >  4:533077/1‑149 (MQ=255)
                                                               gcgcgcCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTg                                                                                    <  3:123524/149‑1 (MQ=255)
                                                                 gcgcCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTggg                                                                                  <  2:161889/149‑1 (MQ=255)
                                                                   gcCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGgaa                                                                                <  3:556564/149‑1 (MQ=255)
                                                                      tgctgGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGgaagaa                                                                             <  1:529561/149‑1 (MQ=255)
                                                                         tgGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAg                                                                          <  3:711142/149‑1 (MQ=255)
                                                                          ggCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAgg                                                                         <  3:649972/149‑1 (MQ=255)
                                                                          ggCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAgg                                                                         <  1:722581/149‑1 (MQ=255)
                                                                           tcGACCAAAGACATCCCCCATGTAATCCAGGTCTGGTATATCTCAACCAAAGAGAAAGCCCAGCAGGCAATTGGCCTAAACATGGAACAAATCAAGGACGAGAAGCATGACTTCAGCAAAGCGGTGATCCCGCACTGGGAAGAACAGGc                                                                        <  2:242154/148‑1 (MQ=255)
                                                                             gACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCAc                                                                      >  1:877957/1‑149 (MQ=255)
                                                                              aCCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGa                                                                                              <  3:787443/124‑1 (MQ=255)
                                                                                cAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACTTa                                                                   <  4:107452/149‑1 (MQ=255)
                                                                                cAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTa                                                                   >  3:624164/1‑149 (MQ=255)
                                                                                 aaaGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAAGGGTGATCCCGCAGTGGGAAGAACAGGCACGTaa                                                                  >  4:761408/1‑149 (MQ=255)
                                                                                  aaGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGc                                                                                                                   >  3:353162/1‑99 (MQ=255)
                                                                                  aaGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGc                                                                                                                   <  4:353162/99‑1 (MQ=255)
                                                                                    ggCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTaaaaa                                                               <  2:441229/149‑1 (MQ=255)
                                                                                     gCATCACCAATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAAc                                                              >  3:797263/1‑149 (MQ=255)
                                                                                        tCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGt                                                           <  1:117186/149‑1 (MQ=255)
                                                                                            ccATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGt                                                           <  4:643647/145‑1 (MQ=255)
                                                                                            ccATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGt                                                       >  4:66962/1‑149 (MQ=255)
                                                                                             cATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGtt                                                      >  2:414826/1‑149 (MQ=255)
                                                                                             cATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGtt                                                      >  1:557564/1‑149 (MQ=255)
                                                                                              aTGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTa                                                     >  1:203241/1‑149 (MQ=255)
                                                                                              aTGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTa                                                     >  1:857964/1‑149 (MQ=255)
                                                                                              aTGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTa                                                     >  2:371059/1‑149 (MQ=255)
                                                                                              aTGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTa                                                     >  1:346535/1‑149 (MQ=255)
                                                                                               tGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTaa                                                    >  1:804477/1‑149 (MQ=255)
                                                                                               tGAAATCCAGATCCCGAATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTaa                                                    >  3:488976/1‑149 (MQ=255)
                                                                                                   aTCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCa                                                >  2:165647/1‑149 (MQ=255)
                                                                                                    tCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGaa                                                                                                                     <  4:253650/79‑1 (MQ=255)
                                                                                                    tCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCaa                                               >  1:313810/1‑149 (MQ=255)
                                                                                                         aTCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAAcc                                          <  3:433599/149‑1 (MQ=255)
                                                                                                            ccGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGt                                                                                                <  2:850548/92‑1 (MQ=255)
                                                                                                            ccGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGt                                                                                                >  1:850548/1‑92 (MQ=255)
                                                                                                               gATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTc                                    <  2:313810/149‑1 (MQ=255)
                                                                                                                   tctcAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGcc                                >  2:738991/1‑149 (MQ=255)
                                                                                                                    ctcAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCt                               >  2:867479/1‑149 (MQ=255)
                                                                                                                    ctcAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCt                               >  4:400282/1‑149 (MQ=255)
                                                                                                                     tcAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTc                              >  2:46117/1‑149 (MQ=255)
                                                                                                                       aaCCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGt                                                                                                <  1:275373/81‑1 (MQ=255)
                                                                                                                       aaCCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGt                                                                                                >  2:275373/1‑81 (MQ=255)
                                                                                                                       aaCCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCt                                                         <  1:212300/120‑1 (MQ=255)
                                                                                                                       aaCCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGc                            >  1:70131/1‑149 (MQ=255)
                                                                                                                            aaGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTCCCGCCTCGCAAggg                       >  3:383298/1‑149 (MQ=255)
                                                                                                                             agagAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAgggg                      <  2:197897/149‑1 (MQ=255)
                                                                                                                              gagaAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGc                     <  4:337192/149‑1 (MQ=255)
                                                                                                                              gagaAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGc                     >  4:33210/1‑149 (MQ=255)
                                                                                                                                gaAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCgg                   >  3:253683/1‑149 (MQ=255)
                                                                                                                                    catgttcagGCCAATGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGaaaa               >  2:786746/7‑149 (MQ=255)
                                                                                                                                     cccAGCAGGCCATTGGCCTGAACATGGAACAAAGCAAGGCCGAGAAGCAGGACTTCATCAAACCGGTGATCTCGCAGTGGGAAGAACAGTCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCGGCCTCGCAAGGGGCGGAAAac              <  2:17011/149‑1 (MQ=255)
                                                                                                                                     atgttcagGCCAATGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAac              >  3:861691/6‑149 (MQ=255)
                                                                                                                                        ttcagGCCAATGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAAcaca           >  1:584280/3‑149 (MQ=255)
                                                                                                                                         gcagGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACaa          <  1:264079/149‑1 (MQ=255)
                                                                                                                                         gcagGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACaa          <  1:251915/149‑1 (MQ=255)
                                                                                                                                              ccATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATGGAg     <  2:681878/149‑1 (MQ=255)
                                                                                                                                              ccATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTCCCGCCTCGCAAGGGGCGGAAAACACAATGGAg     <  4:383298/149‑1 (MQ=255)
                                                                                                                                               cATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATGGAGt    <  3:600311/149‑1 (MQ=255)
                                                                                                                                                aTTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATGGAGTg   >  4:242413/1‑149 (MQ=255)
                                                                                                                                                aTTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATGGAGTg   <  1:738991/149‑1 (MQ=255)
                                                                                                                                                 ttGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATGGAGTGa  >  1:854900/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATGGAGTGA  >  NC_000913/4288949‑4289242

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