Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,643,723 T→C K194K (AAA→AAG ydfU ← Qin prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,643,7230TC98.8% 322.0 / ‑4.6 84K194K (AAA→AAGydfUQin prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (35/48);  total (36/48)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

GATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATC  >  NC_000913/1643576‑1643870
                                                                                                                                                   |                                                                                                                                                   
gATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGAtctt                                                                                                                                                   <  2:112727/150‑1 (MQ=255)
gATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGAtctt                                                                                                                                                   <  2:137818/150‑1 (MQ=255)
   gCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTc                                                                                                                                                <  2:741029/150‑1 (MQ=255)
    cTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCg                                                                                                                                              >  2:809727/1‑151 (MQ=255)
       accaccGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCt                                                                                                                                            <  1:1038891/150‑1 (MQ=255)
       accaccGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCt                                                                                                                                            <  2:451314/150‑1 (MQ=255)
       accaccGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTg                                                                                                                                           <  1:974033/151‑1 (MQ=255)
       accaccGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTg                                                                                                                                           <  1:366151/151‑1 (MQ=255)
        ccaccGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTgg                                                                                                                                          <  2:119771/151‑1 (MQ=255)
           ccGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTa                                                                                                                                        <  2:628510/150‑1 (MQ=255)
             gTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATc                                                                                                                                     <  1:514839/151‑1 (MQ=255)
             gTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATc                                                                                                                                     <  1:3461/151‑1 (MQ=255)
                ttttCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATccc                                                                                                                                   >  2:1138298/1‑150 (MQ=255)
                 tttCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATcc                                                                                                                                    <  2:304637/148‑1 (MQ=255)
                  ttCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAg                                                                                                                                 <  2:420474/150‑1 (MQ=255)
                  ttCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGc                                                                                                                                >  2:187622/1‑151 (MQ=255)
                     cTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGa                                                                                                                              >  2:275278/1‑150 (MQ=255)
                       gtggtgATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAac                                                                                                                            <  2:848987/150‑1 (MQ=255)
                        ttgtcatatttGTTGCTTGTGGGCGTAAGGCGACGGTGCAGTATTTCTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCCCTTTCGCGGTACACCGAGCAGATCTGTTCCGCTGGTAATAACAGc                                                                                                                                <  1:975164/137‑1 (MQ=255)
                             aTTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCAAGCGAACGAcgcg                                                                                                                      >  1:550765/1‑150 (MQ=255)
                                 ttgttgCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATAc                                                                                                                  >  1:517760/1‑150 (MQ=255)
                                   gttgCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAg                                                                                                                <  2:241120/150‑1 (MQ=255)
                                     tgCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTc                                                                                                            >  2:1204662/1‑151 (MQ=255)
                                          gtggGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGt                                                                                                        <  1:1033671/151‑1 (MQ=255)
                                                  aaGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCAt                                                                                                 <  1:124931/150‑1 (MQ=255)
                                                    ggCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATcc                                                                                               >  2:15634/1‑150 (MQ=255)
                                                      cAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATccacc                                                                                            >  1:75740/1‑151 (MQ=255)
                                                       aaCGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATccacc                                                                                            >  1:263397/1‑150 (MQ=255)
                                                       aaCGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATccac                                                                                             >  2:1015091/1‑149 (MQ=255)
                                                           gTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTg                                                                                       <  2:937612/151‑1 (MQ=255)
                                                              cAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTTGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGAtt                                                                                    >  1:989073/1‑151 (MQ=255)
                                                              cAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGAtt                                                                                    <  2:163652/151‑1 (MQ=255)
                                                              cAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGAtt                                                                                    >  2:577883/1‑151 (MQ=255)
                                                               aaGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTg                                                                                   >  1:51063/1‑151 (MQ=255)
                                                                aGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGc                                                                                  <  1:423375/151‑1 (MQ=255)
                                                                aGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGc                                                                                  <  1:28032/151‑1 (MQ=255)
                                                                 gATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTGCGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGc                                                                                  <  2:565037/150‑1 (MQ=255)
                                                                 gATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGGCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGc                                                                                  >  1:819474/1‑150 (MQ=255)
                                                                     ttttGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGAc                                                                              <  1:1007955/150‑1 (MQ=255)
                                                                      tttGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGAcc                                                                             <  1:251125/150‑1 (MQ=255)
                                                                        tgtgCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGAcc                                                                             <  1:393165/148‑1 (MQ=255)
                                                                        tgtgCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGAGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGcc                                                                          >  1:475905/1‑151 (MQ=255)
                                                                         gtgCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGccc                                                                         >  2:803011/1‑151 (MQ=255)
                                                                            cGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCccacc                                                                      <  2:115648/151‑1 (MQ=255)
                                                                             gttgttTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCccacc                                                                      >  1:345351/1‑150 (MQ=255)
                                                                               tgttTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAgg                                                                   >  2:48203/1‑151 (MQ=255)
                                                                                   tCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATa                                                                <  1:486642/150‑1 (MQ=255)
                                                                                     aaTATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCGGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGAATAAtt                                                             <  1:856348/151‑1 (MQ=255)
                                                                                     aaTATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAAtt                                                             <  2:340633/151‑1 (MQ=255)
                                                                                     aaTATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAAtt                                                             <  1:618676/151‑1 (MQ=255)
                                                                                      atatGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCTGGATAAt                                                              >  1:590995/1‑149 (MQ=255)
                                                                                      atatGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTc                                                            <  2:989995/151‑1 (MQ=255)
                                                                                           cTGGTGGCGGTCAGCTCTCCCGGTACGATGGCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGcc                                                        >  1:744790/1‑150 (MQ=255)
                                                                                            tggtggCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGCACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCCCCCGGATAATTCAGcc                                                        >  1:364339/1‑149 (MQ=255)
                                                                                             ggtggCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCaa                                                      >  2:192568/1‑150 (MQ=255)
                                                                                             ggtggCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCaa                                                      <  2:691342/150‑1 (MQ=255)
                                                                                              gtggCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGGAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCaa                                                      >  2:454337/1‑149 (MQ=255)
                                                                                                   ggTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAAt                                               >  1:664146/1‑151 (MQ=255)
                                                                                                      cAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATCCAGCTTCATGTAGTGCATCCACCACCTGATTGTAGACCGCCGACCATGATAATTCAGCCAAAGATAAtt                                              <  1:89080/149‑1 (MQ=255)
                                                                                                      cAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTcc                                            <  2:1105305/151‑1 (MQ=255)
                                                                                                          tctcCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCg                                          >  2:796080/1‑149 (MQ=255)
                                                                                                             cccGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGACCGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTcc                                      <  1:72764/150‑1 (MQ=255)
                                                                                                             cccGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCt                                     >  2:414087/1‑151 (MQ=255)
                                                                                                                ggTACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCg                                  <  1:895077/151‑1 (MQ=255)
                                                                                                                  tACGATGTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGt                                 >  1:431024/1‑150 (MQ=255)
                                                                                                                        gTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCtt                          <  1:50314/151‑1 (MQ=255)
                                                                                                                        gTCGCTTTCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCt                           <  2:1092367/150‑1 (MQ=255)
                                                                                                                              ttCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCg                    >  1:266544/1‑151 (MQ=255)
                                                                                                                               tCGCGGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCg                    <  2:423923/150‑1 (MQ=255)
                                                                                                                                  cgGTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGAc                >  1:209514/1‑151 (MQ=255)
                                                                                                                                   ggTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACg               <  1:1018805/151‑1 (MQ=255)
                                                                                                                                    gTACACCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACgg              <  2:684679/151‑1 (MQ=255)
                                                                                                                                     tACACCGAGCAGATCTTTTCCGCTGGTAACCCCAGCGAACGACGCGATACATCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATACCCGCTCCTGCGTACCGCTTATTGCGTGACgg              <  2:875859/150‑1 (MQ=255)
                                                                                                                                       cacCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATg           <  2:704547/151‑1 (MQ=255)
                                                                                                                                       cacCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATg           >  1:16310/1‑151 (MQ=255)
                                                                                                                                       cacCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATg           >  2:946462/1‑151 (MQ=255)
                                                                                                                                        acCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGATCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGa          <  2:308148/151‑1 (MQ=255)
                                                                                                                                        acCGAGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGa          <  2:974733/151‑1 (MQ=255)
                                                                                                                                            aGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTc      <  2:263403/151‑1 (MQ=255)
                                                                                                                                            aGCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTc      <  2:412594/151‑1 (MQ=255)
                                                                                                                                             gCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCa     <  2:640550/151‑1 (MQ=255)
                                                                                                                                             gCAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCa     <  2:251982/151‑1 (MQ=255)
                                                                                                                                              cAGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCa     >  1:716967/1‑150 (MQ=255)
                                                                                                                                               aGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAAt   >  1:563315/1‑151 (MQ=255)
                                                                                                                                               aGATCTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAAt   >  2:429356/1‑151 (MQ=255)
                                                                                                                                                 aTTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATc  >  2:38246/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGCAGACCGCCCACCAGGATAATTCAGCCAAAGATAATTCCCGCTCCTGCGTACCGCTTATTGCGTGACGGATGACGTCAATC  >  NC_000913/1643576‑1643870

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