Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,994,896 G→C pseudogene (197/633 nt) ygeK ← pseudogene, response regulator family, part of T3SS PAI ETT2 remnant; putative 2‑component transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,994,8960GC99.3% 468.0 / ‑3.3 136pseudogene (197/633 nt)ygeKpseudogene, response regulator family, part of T3SS PAI ETT2 remnant; putative 2‑component transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base C (66/69);  total (66/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

Only 100 of 138 total aligned reads displayed.

ATTTAGTTTCGCCGGGTGACTGGTTTTCCATACGCATCCTTTCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTA  >  NC_000913/2994748‑2995037
                                                                                                                                                    |                                                                                                                                             
aTTTAGTTTCGCCGGGTGACTGGTTTTCCATACGCATCCTTTCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGccc                                                                                                                                             >  4:1517216/1‑151 (MQ=255)
    aGTTTCGCCGGGTGACTGGTTTTCCATACGCATCCTTTCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGcccccc                                                                                                                                          <  1:865896/150‑1 (MQ=255)
      tttCGCCGGGTGACTGGTTTTCCATACGCATCCTTTCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCcatc                                                                                                                                       >  4:1436988/1‑151 (MQ=255)
       ttCGCCGGGTGACTGGTTTTCCATACGCATCCTTTCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCcatc                                                                                                                                       >  1:246535/1‑150 (MQ=255)
       ttCGCCGGGTGACTGGTTTTCCATACGCATCCTTTCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCcatc                                                                                                                                       >  4:1418798/1‑150 (MQ=255)
            cGGGTGACTGGTTTTCCATACGCATCCTTTCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCa                                                                                                                                 <  3:490759/151‑1 (MQ=255)
             gggTGACTGGTTTTCCATACGCATCCTTTCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCAc                                                                                                                                >  2:348538/1‑151 (MQ=255)
               gTGACTGGTTTTCCATACGCATCCTTTCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCt                                                                                                                              <  1:575768/151‑1 (MQ=255)
                 gACTGGTTTTCCATACGCATCCTTTCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTg                                                                                                                             <  4:1154420/150‑1 (MQ=255)
                  aCTGGTTTTCCATACGCATCCTTTCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCa                                                                                                                           <  1:195990/151‑1 (MQ=255)
                     ggTTTTCCATACGCATCCTTTCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCATCATTCCACCTGCAAAt                                                                                                                        <  4:262183/150‑1 (MQ=255)
                      gTTTTCCATACGCATCCTTTCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCAGTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAAt                                                                                                                        <  4:122063/150‑1 (MQ=255)
                          tCCATACGCATCCTTTCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATtctc                                                                                                                    <  1:111759/150‑1 (MQ=255)
                           ccATACGCATCCTTTCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTcc                                                                                                                   <  4:1211522/150‑1 (MQ=255)
                            cATACGCATCCTTTCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTa                                                                                                                 <  1:1334634/151‑1 (MQ=255)
                             aTACGCATCCTTTCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTa                                                                                                                 >  4:505931/1‑150 (MQ=255)
                             aTACGCATCCTTTCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTa                                                                                                                 >  1:1713398/1‑150 (MQ=255)
                               aCGCATCCTTTCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCCTTTTACGAGTTCGCCCCCCAGCATTCCACCTGCAAATTCTCCTaaaa                                                                                                              <  4:148008/151‑1 (MQ=255)
                               aCGCATCCTTTCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTaaa                                                                                                               >  4:1607288/1‑150 (MQ=255)
                                    tCCTTTCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAaaataaat                                                                                                         >  1:1064004/1‑151 (MQ=255)
                                     ccTTTCGCACCTGCTTCAAGCAATTTTGTTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAaaataaat                                                                                                         >  4:1042891/1‑150 (MQ=255)
                                        ttCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATAttt                                                                                                     >  1:1386394/1‑151 (MQ=255)
                                        ttCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATAtt                                                                                                      <  3:504888/150‑1 (MQ=255)
                                         tCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATAttt                                                                                                     <  1:419473/150‑1 (MQ=255)
                                          cGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAGATTCTCCTAAAATAAATATTTa                                                                                                    <  4:353299/150‑1 (MQ=255)
                                           gCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAgg                                                                                                  >  1:1342279/1‑151 (MQ=255)
                                            cACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGc                                                                                                 <  4:1071260/151‑1 (MQ=255)
                                             aCCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCt                                                                                                <  1:725131/151‑1 (MQ=255)
                                              ccTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTg                                                                                               <  4:553459/151‑1 (MQ=255)
                                               cTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGt                                                                                              >  4:1867599/1‑151 (MQ=255)
                                                    tCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGa                                                                                           <  4:1082907/149‑1 (MQ=255)
                                                    tCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGa                                                                                           <  2:411626/149‑1 (MQ=255)
                                                     cAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGAttt                                                                                        <  4:1125649/151‑1 (MQ=255)
                                                     cAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAACTCTCCTAAAATAAATATTTAGGCTGTAGAttt                                                                                        <  4:665497/151‑1 (MQ=255)
                                                      aaGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGAttt                                                                                        >  4:617981/1‑150 (MQ=255)
                                                      aaGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTg                                                                                       >  2:468206/1‑151 (MQ=255)
                                                       aGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGt                                                                                      <  1:149269/151‑1 (MQ=255)
                                                         cAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGtt                                                                                     >  4:824592/1‑150 (MQ=255)
                                                         cAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGtt                                                                                     >  1:1830131/1‑150 (MQ=255)
                                                           aTTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTac                                                                                   >  1:939052/1‑150 (MQ=255)
                                                            ttttATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTacac                                                                                 >  1:28284/1‑151 (MQ=255)
                                                                aTTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTcc                                                                              <  1:285601/150‑1 (MQ=255)
                                                                    aaTCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCtat                                                                         <  4:547378/151‑1 (MQ=255)
                                                                     aTCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCtata                                                                        <  4:631090/151‑1 (MQ=255)
                                                                       cAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAg                                                                       <  1:1316662/150‑1 (MQ=255)
                                                                         atatACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCtt                                                                    >  4:1900878/1‑151 (MQ=255)
                                                                          tataCGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCttt                                                                   >  1:493148/1‑151 (MQ=255)
                                                                          tataCGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCttt                                                                   >  2:111686/1‑151 (MQ=255)
                                                                          tataCGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCttt                                                                   >  4:1373217/1‑151 (MQ=255)
                                                                               ggCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTttcttc                                                               >  4:750836/1‑150 (MQ=255)
                                                                                gCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAg                                                             <  1:729325/151‑1 (MQ=255)
                                                                                 cATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGa                                                            <  1:1537303/151‑1 (MQ=255)
                                                                                  aTCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTttcttc                                                               >  2:452778/1‑147 (MQ=255)
                                                                                      ttGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGAt                                                           <  3:456248/147‑1 (MQ=255)
                                                                                      ttGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATc                                                          <  1:1263386/148‑1 (MQ=255)
                                                                                       tGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCttt                                                       <  1:960383/150‑1 (MQ=255)
                                                                                        gCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTtat                                                     >  1:1242866/1‑151 (MQ=255)
                                                                                        gCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTta                                                      <  3:192135/150‑1 (MQ=255)
                                                                                        gCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTta                                                      <  1:1372056/150‑1 (MQ=255)
                                                                                         cTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTtat                                                     >  4:1478710/1‑150 (MQ=255)
                                                                                         cTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTtat                                                     >  1:1473112/1‑150 (MQ=255)
                                                                                              aaTGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGcccc                                               <  1:824237/151‑1 (MQ=255)
                                                                                              aaTGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGcccc                                               >  1:557965/1‑151 (MQ=255)
                                                                                               aTGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGccccc                                              >  4:41075/1‑151 (MQ=255)
                                                                                                 gTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCcaa                                            <  4:1671377/151‑1 (MQ=255)
                                                                                                    ataataTGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCcaacaa                                         <  4:678145/151‑1 (MQ=255)
                                                                                                      aataTGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAAcc                                       >  4:1573174/1‑151 (MQ=255)
                                                                                                       atatGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCt                                      <  1:1372658/151‑1 (MQ=255)
                                                                                                         atGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCa                                    <  1:1520256/151‑1 (MQ=255)
                                                                                                           gAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCata                                  <  1:187985/151‑1 (MQ=255)
                                                                                                            aGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCatat                                 >  4:43582/1‑151 (MQ=255)
                                                                                                                tctatctTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAAt                             >  4:1758995/1‑151 (MQ=255)
                                                                                                                 ctatctTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATc                            >  4:851482/1‑151 (MQ=255)
                                                                                                                  tatctTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATc                            >  1:1441702/1‑150 (MQ=255)
                                                                                                                  tatctTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATc                            >  4:68878/1‑150 (MQ=255)
                                                                                                                  tatctTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATc                            >  4:1739571/1‑150 (MQ=255)
                                                                                                                  tatctTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATc                            >  2:329415/1‑150 (MQ=255)
                                                                                                                  tatctTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATc                            >  4:1042377/1‑150 (MQ=255)
                                                                                                                   atctTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATc                            <  2:215299/149‑1 (MQ=255)
                                                                                                                   atctTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATc                            <  1:690018/149‑1 (MQ=255)
                                                                                                                     ctTATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAgg                         >  4:233616/1‑150 (MQ=255)
                                                                                                                       tATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAgg                         <  4:1802056/148‑1 (MQ=255)
                                                                                                                       tATGCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGAt                       <  3:111681/150‑1 (MQ=255)
                                                                                                                          gCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAAt                    <  4:452738/150‑1 (MQ=255)
                                                                                                                          gCGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATg                   >  1:975133/1‑151 (MQ=255)
                                                                                                                           cGATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGt                  >  1:545852/1‑151 (MQ=255)
                                                                                                                            gATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTc                 >  1:1434693/1‑151 (MQ=255)
                                                                                                                            gATTTAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTc                 <  4:891655/151‑1 (MQ=255)
                                                                                                                               ttAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCg               <  1:1866877/150‑1 (MQ=255)
                                                                                                                                tAACCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCga              <  4:1325094/150‑1 (MQ=255)
                                                                                                                                  aCCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCgaga            >  1:1115226/1‑150 (MQ=255)
                                                                                                                                  aCCCATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCgaga            >  2:303206/1‑150 (MQ=255)
                                                                                                                                   cccATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCGAGaa           >  4:627648/1‑150 (MQ=255)
                                                                                                                                    ccATTTTACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCGAGAAAt         >  4:1835624/1‑151 (MQ=255)
                                                                                                                                       ttttACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCGAGAAATcc       >  1:909171/1‑150 (MQ=255)
                                                                                                                                        tttACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCGAGAAATCCa      <  1:1452312/150‑1 (MQ=255)
                                                                                                                                         ttACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCGAGAAATCCaa     <  4:1421936/150‑1 (MQ=255)
                                                                                                                                         ttACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAAt    >  2:71794/1‑151 (MQ=255)
                                                                                                                                          tACGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAAt    <  1:293101/150‑1 (MQ=255)
                                                                                                                                           aCGAGTTCGCCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTa  <  2:26831/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
ATTTAGTTTCGCCGGGTGACTGGTTTTCCATACGCATCCTTTCGCACCTGCTTCAAGCAATTTTATTGAATCAATATACGGCATCTTTGCTACAAATGTAATAATATGAGCATCTATCTTATGCGATTTAACCCATTTTACGAGTTCGGCCCCCATCATTCCACCTGCAAATTCTCCTAAAATAAATATTTAGGCTGTAGATTTGTTACACTCCGCTATAGCTTTCTTCAGATCTTTATAGCCCCCAACAACCTCATATAAATCAGGATAATGTCCGAGAAATCCAATTA  >  NC_000913/2994748‑2995037

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