Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,630,921 A→T 6.1% intergenic (+56/+37) yfgI → / ← guaA nalidixic acid resistance protein, putative periplasmic protein/GMP synthetase (glutamine aminotransferase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,630,9210AT 6.1% 436.5 / 12.3 128intergenic (+56/+37)yfgI/guaAnalidixic acid resistance protein, putative periplasmic protein/GMP synthetase (glutamine aminotransferase)
Reads supporting (aligned to +/- strand):  ref base A (61/59);  new base T (2/6);  total (63/65)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 138 total aligned reads displayed.

GCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATCGTAC  >  NC_000913/2630776‑2631069
                                                                                                                                                 |                                                                                                                                                    
gCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACa                                                                                                                                                   <  4:744098/149‑1 (MQ=255)
   aaCAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACagag                                                                                                                                                >  2:577430/1‑149 (MQ=255)
   aaCAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACagag                                                                                                                                                >  3:467866/1‑149 (MQ=255)
   aaCAGAAATATGACTGAAGGCGGTAGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACagag                                                                                                                                                >  3:265864/1‑149 (MQ=255)
    aCAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAgg                                                                                                                                               >  2:181137/1‑149 (MQ=255)
    aCAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCAGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAgg                                                                                                                                               <  1:497382/149‑1 (MQ=255)
     cAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTTACTACAGAggg                                                                                                                                              <  1:743944/149‑1 (MQ=255)
       gAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGtt                                                                                                                                            >  4:871900/1‑149 (MQ=255)
          atatGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTg                                                                                                                                         <  3:650660/149‑1 (MQ=255)
          atatGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTg                                                                                                                                         <  2:558949/149‑1 (MQ=255)
             tGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGttc                                                                                                                                      <  4:115399/149‑1 (MQ=255)
               aCTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGttcat                                                                                                                                    <  1:591171/149‑1 (MQ=255)
                cTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGttcatt                                                                                                                                   <  4:610416/149‑1 (MQ=255)
                 tGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTTACTACAGAggg                                                                                                                                              <  1:343809/137‑1 (MQ=255)
                 tGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGtttt                                                                                                                                          <  3:557729/141‑1 (MQ=255)
                   aaGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTtcattcat                                                                                                                                >  1:311033/1‑149 (MQ=255)
                    aGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATa                                                                                                                               >  1:877929/1‑149 (MQ=255)
                            gcgATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAggg                                                                                                                       >  4:761534/1‑149 (MQ=255)
                              gATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTc                                                                                                                     >  1:413652/1‑149 (MQ=255)
                              gATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTc                                                                                                                     >  4:976985/1‑149 (MQ=255)
                                 cTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCaaa                                                                                                                  <  3:367164/149‑1 (MQ=255)
                                  tgatgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTTACTACAGAggg                                                                                                                                              <  4:574680/120‑1 (MQ=255)
                                  tgatgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAAt                                                                                                                 >  3:96484/1‑149 (MQ=255)
                                  tgatgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAAACCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAgg                                                                                                                                               <  1:526692/119‑1 (MQ=255)
                                    atgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCa                                                                                                               <  1:6778/149‑1 (MQ=255)
                                    atgaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCa                                                                                                               >  1:146084/1‑149 (MQ=255)
                                      gaAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCAtt                                                                                                             <  3:412908/149‑1 (MQ=255)
                                        aaaGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTcc                                                                                                           >  2:971119/1‑149 (MQ=255)
                                         aaGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTccc                                                                                                          <  3:871900/149‑1 (MQ=255)
                                           gCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCAc                                                                                                        <  4:724752/149‑1 (MQ=255)
                                             gTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGtttt                                                                                                                                          >  1:413738/1‑113 (MQ=255)
                                             gTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGtttt                                                                                                                                          <  2:413738/113‑1 (MQ=255)
                                                ggCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAAt                                                                                                   <  2:273474/149‑1 (MQ=255)
                                                   aaaGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGt                                                                                                >  1:105664/1‑149 (MQ=255)
                                                     aGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGttc                                                                                                                                      <  1:285901/109‑1 (MQ=255)
                                                     aGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGttc                                                                                                                                      >  2:285901/1‑109 (MQ=255)
                                                              gaACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCt                                                                                     >  2:49827/1‑149 (MQ=255)
                                                              gaACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCt                                                                                     <  3:251614/149‑1 (MQ=255)
                                                               aaCAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGttcat                                                                                                                                    >  2:1014617/1‑101 (MQ=255)
                                                               aaCAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGttcat                                                                                                                                    <  1:1014617/101‑1 (MQ=255)
                                                                 cAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGc                                                                                  <  3:731008/149‑1 (MQ=255)
                                                                       gtaCAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGTTAAAAAACCCTCTGTAGTAACAGAGAGCGGTGTTCAGTCATACTGCAGCGTCTAATCATGCCCACTCACTGGTAGCTGGCGGCTTGCCggtta                                                                            <  2:700530/147‑5 (MQ=255)
                                                                       ggACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGa                                                                            <  2:413652/149‑1 (MQ=255)
                                                                        gACATCTAGTTTAATTAATCGATGTTAGTAACTTCAATGCTGTAATTCTTGAAGATAAAAAACCCTCTGAAGTAACAGCGGGTTATCTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTCGCTTGTGGCTTGCCGCTGAt                                                                           <  2:763804/149‑1 (MQ=255)
                                                                        gACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGAt                                                                           >  3:622659/1‑149 (MQ=255)
                                                                            gCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTTAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAACTCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCa                                                                       <  4:41843/149‑1 (MQ=255)
                                                                            gCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCa                                                                       <  4:670725/149‑1 (MQ=255)
                                                                            gCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCa                                                                       <  1:77320/149‑1 (MQ=255)
                                                                                     attaatCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATa                                                                     <  3:791817/142‑1 (MQ=255)
                                                                                     attaatCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCAGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCAcg                                                              >  4:671934/1‑149 (MQ=255)
                                                                                        aatCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCgg                                                           <  4:881750/149‑1 (MQ=255)
                                                                                         atCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCggg                                                          >  1:1010224/1‑149 (MQ=255)
                                                                                          tCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGTGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGa                                                         >  3:18394/1‑149 (MQ=255)
                                                                                          tCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGa                                                         <  1:1036776/149‑1 (MQ=255)
                                                                                          tCGATGTTAGTAACTTCAATCCTATAAGTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAAGGGTAGCTGGCGGCTTGCCGCTGAGGTCATACACCACGCGGGa                                                         >  2:1048359/1‑149 (MQ=255)
                                                                                            gATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATa                                                                     <  3:1008678/135‑1 (MQ=255)
                                                                                            gATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATa                                                                     >  4:1008678/1‑135 (MQ=255)
                                                                                               gTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGttc                                                                                                                                      <  4:818001/67‑1 (MQ=255)
                                                                                               gTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGttc                                                                                                                                      >  3:818001/1‑67 (MQ=255)
                                                                                               gTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATAc                                                    <  2:215325/149‑1 (MQ=255)
                                                                                                  aGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCa                                                                                                         <  2:920504/93‑1 (MQ=255)
                                                                                                  aGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGt                                                 <  1:1005527/149‑1 (MQ=255)
                                                                                                   gTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATCGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGtt                                                >  1:580653/1‑149 (MQ=255)
                                                                                                   gTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGtt                                                <  2:492342/149‑1 (MQ=255)
                                                                                                        ttCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACtt                                           >  4:281630/1‑149 (MQ=255)
                                                                                                        ttCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACtt                                           >  1:596169/1‑149 (MQ=255)
                                                                                                        ttCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACtt                                           >  1:695207/1‑149 (MQ=255)
                                                                                                        ttCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACtt                                           <  4:96484/149‑1 (MQ=255)
                                                                                                           aaTCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCAt                                        >  1:491487/1‑149 (MQ=255)
                                                                                                            aTCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCAtt                                       <  1:223860/149‑1 (MQ=255)
                                                                                                                tataATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATa                                   >  4:103558/1‑149 (MQ=255)
                                                                                                                 ataATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCAtt                                       >  2:301147/1‑144 (MQ=255)
                                                                                                                 ataATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCAtt                                       <  1:301147/144‑1 (MQ=255)
                                                                                                                 ataATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATaa                                  <  3:230833/149‑1 (MQ=255)
                                                                                                                   aaTTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATg                                >  2:672131/1‑149 (MQ=255)
                                                                                                                       cTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGt                            <  3:103558/149‑1 (MQ=255)
                                                                                                                        ttGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTg                                                                                           <  3:562734/85‑1 (MQ=255)
                                                                                                                        ttGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTg                                                                                           >  4:562734/1‑85 (MQ=255)
                                                                                                                        ttGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGtt                           <  2:580653/149‑1 (MQ=255)
                                                                                                                        ttGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGtt                           >  1:655875/1‑149 (MQ=255)
                                                                                                                             gATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGaaa                      <  2:372448/149‑1 (MQ=255)
                                                                                                                              aTAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAAc                     <  2:368907/149‑1 (MQ=255)
                                                                                                                                 aaaaaCCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAAcgcg                  >  2:906368/1‑149 (MQ=255)
                                                                                                                                   aaaCCCTCTGTTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCAtt                                       >  4:177660/1‑126 (MQ=255)
                                                                                                                                   aaaCCCTCTGTTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCAtt                                       <  3:177660/126‑1 (MQ=255)
                                                                                                                                    aaCCCTCTGTTACTACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATAcc                                                   <  2:841542/113‑1 (MQ=255)
                                                                                                                                      cccTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGa             <  4:467866/149‑1 (MQ=255)
                                                                                                                                       ccTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAg            >  1:1040711/1‑149 (MQ=255)
                                                                                                                                        ctctGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAgg           >  4:502634/1‑149 (MQ=255)
                                                                                                                                         tctGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGa          >  3:937199/1‑149 (MQ=255)
                                                                                                                                         tctGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGa          >  2:181999/1‑149 (MQ=255)
                                                                                                                                             tagtaACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATc      >  3:333311/1‑149 (MQ=255)
                                                                                                                                             tagtaACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATc      <  1:181999/149‑1 (MQ=255)
                                                                                                                                             tagtaACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATc      <  2:410808/149‑1 (MQ=255)
                                                                                                                                               gtaACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATc      <  1:1033430/147‑1 (MQ=255)
                                                                                                                                               gtaACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATc      >  2:1033430/1‑147 (MQ=255)
                                                                                                                                                taACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATCGTa   <  3:137736/149‑1 (MQ=255)
                                                                                                                                                 aaCAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATCGTAc  <  4:265864/149‑1 (MQ=255)
                                                                                                                                                 aaCAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATCGTAc  <  4:118676/149‑1 (MQ=255)
                                                                                                                                                 aaCAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATCGTAc  <  1:794022/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GCCAACAGAAATATGACTGAAGGCGGTCGCGATCTGATGAAAAGCGTGGGCAAAGCAGAAGAGAACAAATCGGACAGCTGGTTTAATTAATCGATGTTAGTAACTTCAATCCTATAATTCTTGAAGATAAAAAACCCTCTGTAGTAACAGAGGGTTTTGTTCATTCATAGTGCAGGGTCAAATCATTCCCACTCAATGGTAGCTGGCGGCTTGCCGCTGATGTCATACACCACGCGGGAAATACCGTTCACTTCATTGATAATGCGGTTGGAAACGCGACCGAGGAAATCGTAC  >  NC_000913/2630776‑2631069

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