Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 65,814 T→A 7.7% intergenic (‑34/+41) polB ← / ← araD DNA polymerase II/L‑ribulose‑5‑phosphate 4‑epimerase AraD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091365,8140TA 7.7% 318.8 / 12.1 104intergenic (‑34/+41)polB/araDDNA polymerase II/L‑ribulose‑5‑phosphate 4‑epimerase AraD
Reads supporting (aligned to +/- strand):  ref base T (47/49);  new base A (3/5);  total (50/54)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 109 total aligned reads displayed.

TGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCC  >  NC_000913/65665‑65962
                                                                                                                                                     |                                                                                                                                                    
tGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGtat                                                                                                                                                      <  2:1451339/150‑1 (MQ=255)
 gCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGtata                                                                                                                                                     >  2:68491/1‑150 (MQ=255)
    gTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGcc                                                                                                                                                  <  2:710472/150‑1 (MQ=255)
      gCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTg                                                                                                                                                >  2:1501785/1‑150 (MQ=255)
         aGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGtt                                                                                                                                             <  1:345787/150‑1 (MQ=255)
          gCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGttt                                                                                                                                            >  1:1269893/1‑150 (MQ=255)
          gCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGttt                                                                                                                                            >  2:1461573/1‑150 (MQ=255)
          gCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTATACTCAAGCCTGGttt                                                                                                                                            <  1:1225774/150‑1 (MQ=255)
            gTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCg                                                                                                                                          <  1:135254/150‑1 (MQ=255)
            gTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTATACTCAAGCCTGGttttt                                                                                                                                          <  1:155457/150‑3 (MQ=255)
             tAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTATACTCAAGCCTGGtttttt                                                                                                                                         <  1:227172/150‑4 (MQ=255)
              aaCCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGtt                                                                                                                                        <  1:1403386/150‑1 (MQ=255)
              aaCCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGtt                                                                                                                                        <  2:316576/150‑1 (MQ=255)
              aaCCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGtt                                                                                                                                        <  1:283666/150‑1 (MQ=255)
              aaCCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTATACTCAAGCCTGGttttttg                                                                                                                                        <  1:576601/150‑5 (MQ=255)
               aCCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGttt                                                                                                                                       <  1:216208/150‑1 (MQ=255)
                ccTGCAACGGCTCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTg                                                                                                                                      >  2:126038/1‑150 (MQ=255)
                ccTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTg                                                                                                                                      <  1:1490223/150‑1 (MQ=255)
                ccTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTg                                                                                                                                      >  1:1510728/1‑150 (MQ=255)
                ccTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTg                                                                                                                                      <  1:1035495/150‑1 (MQ=255)
                ccTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCACTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTg                                                                                                                                      >  1:163826/1‑150 (MQ=255)
                  tGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGCGTTAAGATAAAACCTGCCTGCGCCACGCTCAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGAt                                                                                                                                    >  1:1396030/1‑150 (MQ=255)
                   gCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGAtt                                                                                                                                   <  2:66722/150‑1 (MQ=255)
                    cAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTg                                                                                                                                  <  2:612480/150‑1 (MQ=255)
                    cAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCCTTTGATTg                                                                                                                                  >  2:382467/1‑150 (MQ=255)
                               tGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGtttt                                                                                                                       <  1:515895/150‑1 (MQ=255)
                                gTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTta                                                                                                                      <  2:348100/150‑1 (MQ=255)
                                 tCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCTGGTTGTCCCGCCAGTGTCTGGTTAAGATAAAACCTTCCTGCTCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTtat                                                                                                                     <  1:757708/150‑1 (MQ=255)
                                   cGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTTTAc                                                                                                                   <  1:876891/150‑1 (MQ=255)
                                   cGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATAc                                                                                                                   <  2:1314315/150‑1 (MQ=255)
                                      cGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGt                                                                                                                >  2:1313599/1‑150 (MQ=255)
                                       gccagccaGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAGGATAAAACCCGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTc                                                                                                               >  1:591695/1‑150 (MQ=255)
                                       gccagccaGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTc                                                                                                               >  1:207717/1‑150 (MQ=255)
                                        ccagccagAAGGACACTTCTGTCCCTTGCGCGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACAAGCCACCGCCACGCTGAAAATCCATCAAAAAACCACGCTTGAGAATAGCCTGCTTTAGGTTGATTGGATGTGGTTTTATACAGTCa                                                                                                              >  2:411868/1‑150 (MQ=255)
                                         cagccagAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGc                                                                                                                                <  1:1435550/131‑1 (MQ=255)
                                         cagccagAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGc                                                                                                                                >  2:1435550/1‑131 (MQ=255)
                                          agccagAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCAtt                                                                                                            >  2:331654/1‑150 (MQ=255)
                                             cagAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACt                                                                                                         <  2:1269893/150‑1 (MQ=255)
                                                aaGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGcc                                                                                                      >  1:1301704/1‑150 (MQ=255)
                                                       cTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAata                                                                                               <  1:747298/150‑1 (MQ=255)
                                                         tCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATg                                                                                             <  1:1343386/150‑1 (MQ=255)
                                                         tCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATg                                                                                             >  1:814208/1‑150 (MQ=255)
                                                          cTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGc                                                                                            >  1:1170285/1‑150 (MQ=255)
                                                            gTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCt                                                                                          >  1:1165656/1‑150 (MQ=255)
                                                              cccTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTc                                                                                        >  1:249975/1‑150 (MQ=255)
                                                              cccTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCCGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTc                                                                                        >  1:1565191/1‑150 (MQ=255)
                                                               ccTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTcg                                                                                       <  2:923597/150‑1 (MQ=255)
                                                                cTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTcgc                                                                                      >  1:662304/1‑150 (MQ=255)
                                                                cTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTcgc                                                                                      >  2:24535/1‑150 (MQ=255)
                                                                   gCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGcc                                                                                   <  2:640602/150‑1 (MQ=255)
                                                                    cGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCa                                                                                  <  2:1325110/150‑1 (MQ=255)
                                                                     ggggTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCAt                                                                                 <  2:207717/150‑1 (MQ=255)
                                                                      gggTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATg                                                                                <  1:845904/150‑1 (MQ=255)
                                                                        gtgtCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCt                                                                              <  2:1592555/150‑1 (MQ=255)
                                                                            cccGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACg                                                                          >  2:1514687/1‑150 (MQ=255)
                                                                             ccGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGc                                                                         <  1:234129/150‑1 (MQ=255)
                                                                             ccGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGc                                                                         >  1:51288/1‑150 (MQ=255)
                                                                             ccGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGc                                                                         <  1:677532/150‑1 (MQ=255)
                                                                               gCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAg                                                                       >  2:712316/1‑150 (MQ=255)
                                                                                ccAGTGTCGGGGTAAGATAAAACCTGCCTGCGCCACGCTGAAACTCCATCAAAAACCCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGa                                                                      <  2:1594841/150‑1 (MQ=255)
                                                                                 cAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGAt                                                                     <  1:1313599/150‑1 (MQ=255)
                                                                                  aGTGTCGGGTTAGGATAAAACCCGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATa                                                                    <  2:591695/150‑1 (MQ=255)
                                                                                      tCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTGGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGATTACGCAGATAgtgt                                                                >  2:425860/1‑150 (MQ=255)
                                                                                          gTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTAt                                                            >  2:139965/1‑150 (MQ=255)
                                                                                            tAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATcc                                                          >  1:1105141/1‑150 (MQ=255)
                                                                                            tAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATcc                                                          >  2:1281002/1‑150 (MQ=255)
                                                                                             aaGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCca                                                         >  2:1177006/1‑150 (MQ=255)
                                                                                              aGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCTATGCTTACGCAGATAGTGTTTATCcag                                                        <  1:24535/150‑1 (MQ=255)
                                                                                              aGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCcag                                                        >  1:636805/1‑150 (MQ=255)
                                                                                              aGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATACCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCcag                                                        <  1:544859/150‑1 (MQ=255)
                                                                                                 tAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCcagcag                                                     >  1:519300/1‑150 (MQ=255)
                                                                                                        tgcctgcGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTTATTACTGCCCGTAATATGCCTTCGCGCCATACTTACGCAGATAGTGTTTATCCAGCAGCGTTtgc                                              >  2:1015760/1‑150 (MQ=255)
                                                                                                        tgcctgcGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTtgc                                              <  1:382467/150‑1 (MQ=255)
                                                                                                        tgcctgcGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTtgc                                              <  1:126038/150‑1 (MQ=255)
                                                                                                          cctgcGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTtgctg                                            <  1:1281002/150‑1 (MQ=255)
                                                                                                          cctgcGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTtgctg                                            <  2:366949/150‑1 (MQ=255)
                                                                                                           ctgcGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTtgctgc                                           <  2:1308052/150‑1 (MQ=255)
                                                                                                                ccACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATc                                      >  1:576633/1‑150 (MQ=255)
                                                                                                                     cTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTa                                 <  2:814208/150‑1 (MQ=255)
                                                                                                                       gAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAAc                               <  2:1396030/150‑1 (MQ=255)
                                                                                                                        aaaaTCCATCGAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTACCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACt                              >  2:1219307/1‑150 (MQ=255)
                                                                                                                         aaaTCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTg                             >  1:1202996/1‑150 (MQ=255)
                                                                                                                          aaTTCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTGGTTTGATTGGCCGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATTCTTACTCAGATAGTGTTTATCCAGCAGCGTTTTCTGCATATCCGGTAACTgc                            <  2:1565191/150‑1 (MQ=255)
                                                                                                                              cATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTgcggcg                        >  1:1329983/1‑150 (MQ=255)
                                                                                                                               aTCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGgcgc                       <  1:712316/150‑1 (MQ=255)
                                                                                                                                 cAAAAAACCAGGCTTGAGTATAGCCTGGTTTCTTTTGATTGGCTGTGGTTTTATAAAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTa                     >  2:1369332/1‑150 (MQ=255)
                                                                                                                                 cAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTa                     >  1:727916/1‑150 (MQ=255)
                                                                                                                                 aacgaaaCCAGGCTATACTCAAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTa                     >  2:888072/5‑150 (MQ=255)
                                                                                                                                  acgaaaCCAGGCTATACTCAAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTaa                    >  2:358798/4‑150 (MQ=255)
                                                                                                                                  acgaaaCCAGGCTATACTCAAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTaa                    >  1:944207/4‑150 (MQ=255)
                                                                                                                                   aaaaaCCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAAc                   <  2:51288/150‑1 (MQ=255)
                                                                                                                                      aaCCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGa                <  2:1052319/150‑1 (MQ=255)
                                                                                                                                        ccAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACg              >  2:470606/1‑150 (MQ=255)
                                                                                                                                          aGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGc            <  2:1105141/150‑1 (MQ=255)
                                                                                                                                          aGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGc            <  1:1502076/150‑1 (MQ=255)
                                                                                                                                            gCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAg          >  1:1298916/1‑150 (MQ=255)
                                                                                                                                               tGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAat       <  1:1177006/150‑1 (MQ=255)
                                                                                                                                                gAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAata      >  2:1312615/1‑150 (MQ=255)
                                                                                                                                                   tataGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATAACCGGTAACTGCGGCGCTAACTGACGGCAGAATATcc   <  1:1501785/150‑1 (MQ=255)
                                                                                                                                                    ataGCCTGGTTTCGTTTGTTTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATccc  <  2:558246/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGCGGTGCAAGCGTAACCTGCAACGGCCCGTTGTCCGTCGCCAGCCAGAAGGAGACTTCTGTCCCTTGCGGGGTGTCCCGCCAGTGTCGGGTTAAGATAAAACCTGCCTGCGCCACGCTGAAAATCCATCAAAAAACCAGGCTTGAGTATAGCCTGGTTTCGTTTGATTGGCTGTGGTTTTATACAGTCATTACTGCCCGTAATATGCCTTCGCGCCATGCTTACGCAGATAGTGTTTATCCAGCAGCGTTTGCTGCATATCCGGTAACTGCGGCGCTAACTGACGGCAGAATATCCC  >  NC_000913/65665‑65962

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