Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,961,248 A→T 5.8% L394Q (CTG→CAG)  recC ← exodeoxyribonuclease V subunit RecC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,961,2480AT 5.8% 272.1 / 4.2 85L394Q (CTG→CAG) recCexodeoxyribonuclease V subunit RecC
Reads supporting (aligned to +/- strand):  ref base A (36/43);  new base T (2/3);  total (39/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTTTCCAGT  >  NC_000913/2961099‑2961397
                                                                                                                                                     |                                                                                                                                                     
cACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGAGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTAGATAACAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGAATCGCca                                                                                                                                                       >  1:1304494/1‑150 (MQ=255)
    acgGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCcagcag                                                                                                                                                   >  1:1665765/1‑150 (MQ=255)
            aaaTGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCCCCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCTGCTGGCGAtgct                                                                                                                                           <  1:630091/150‑4 (MQ=255)
            aaaTGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCTGCTGGCGAtgct                                                                                                                                           <  1:866999/150‑4 (MQ=255)
              aTGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCTGCTGGCGAtgctgg                                                                                                                                         <  1:1647230/150‑6 (MQ=255)
                 gCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTTGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTaaa                                                                                                                                      >  1:612379/1‑150 (MQ=255)
                    tAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAaact                                                                                                                                   >  1:1521530/1‑150 (MQ=255)
                      aggtaggtAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAaacttc                                                                                                                                 <  1:1027104/150‑1 (MQ=255)
                      aggtaggtAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAaacttc                                                                                                                                 >  1:1448081/1‑150 (MQ=255)
                       ggtaggtaACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAaacttca                                                                                                                                <  1:891641/150‑1 (MQ=255)
                         taggtaACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAaacttcaac                                                                                                                              <  2:388255/150‑1 (MQ=255)
                           ggtaACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAaacttcaactt                                                                                                                            <  1:526971/150‑1 (MQ=255)
                            gtaACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGGCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAaacttcaacttc                                                                                                                           <  2:1613960/150‑1 (MQ=255)
                               aCGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACg                                                                                                                        >  1:1472151/1‑150 (MQ=255)
                                   tCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAaacttcaactt                                                                                                                            <  2:1532227/142‑1 (MQ=255)
                                   tCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAaacttcaactt                                                                                                                            >  1:1532227/1‑142 (MQ=255)
                                    ccGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCg                                                                                                                   >  2:809161/1‑150 (MQ=255)
                                    ccGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGAATGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCg                                                                                                                   >  2:815079/1‑150 (MQ=255)
                                     cGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCCGCAGGCGATCGTGAAAACCATCAACTTCACGCTGCgg                                                                                                                  >  2:254403/1‑150 (MQ=255)
                                       cAGGTGCGCTACCAAACACAGCCTTAATAAACGGACTGTATCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGc                                                                                                                <  2:1243119/150‑1 (MQ=255)
                                           tGCGCTACCAAACACAGCCTTAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTGTTCCTCCAGCATCGCGAGCAGGCGATCGTTTAAAACTTGAACTTCACGCTTCGGGCTATg                                                                                                            <  2:725547/150‑1 (MQ=255)
                                           tGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATg                                                                                                            >  1:1530633/1‑150 (MQ=255)
                                            gcgcTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATgg                                                                                                           >  1:1552464/1‑150 (MQ=255)
                                              gcTACCAAACACAGCCTGAATAACCGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCa                                                                                                         <  1:733550/150‑1 (MQ=255)
                                                 aCCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAaacttca                                                                                                                                <  1:566623/124‑1 (MQ=255)
                                                 aCCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAaacttca                                                                                                                                >  2:566623/1‑124 (MQ=255)
                                                 aCCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAAc                                                                                                      <  1:1125415/150‑1 (MQ=255)
                                                   cAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGt                                                                                                    <  1:1143675/150‑1 (MQ=255)
                                                    aaaCACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTg                                                                                                   <  1:339391/150‑1 (MQ=255)
                                                     aaCACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTgg                                                                                                  <  1:542636/150‑1 (MQ=255)
                                                       cacaGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGaa                                                                                                >  1:714215/1‑150 (MQ=255)
                                                       cacaGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGaa                                                                                                >  1:730650/1‑150 (MQ=255)
                                                           gCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTg                                                                                            <  2:1328810/150‑1 (MQ=255)
                                                           gCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTg                                                                                            >  1:1498503/1‑150 (MQ=255)
                                                           gCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTg                                                                                            <  2:242417/150‑1 (MQ=255)
                                                           gCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTg                                                                                            <  2:94773/150‑1 (MQ=255)
                                                             cTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGAt                                                                                          <  1:609532/150‑1 (MQ=255)
                                                             cTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGAt                                                                                          <  1:615396/150‑1 (MQ=255)
                                                               gAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATAc                                                                                        >  1:1022123/1‑150 (MQ=255)
                                                               gAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATAc                                                                                        >  2:995494/1‑150 (MQ=255)
                                                                  tAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGc                                                                                     <  1:683151/150‑1 (MQ=255)
                                                                  tAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGc                                                                                     <  1:1656501/150‑1 (MQ=255)
                                                                       ggAGTGTAGCTGTCTAAATCAGCGACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATcc                                                                                <  2:925672/150‑1 (MQ=255)
                                                                          cTGTAGCTGTCGATATGAGCCACGATCACGATGATGGCTCGCGGAGTATGTGTGGGGGCTTCCTCCAGCATGGCCAGGAGGCGATCGTGTAAAACTTGAACTTCACGATGCGGGCTATGGCAAACGTGGAAGGTGATACTGGTATCCAGt                                                                             <  2:360180/150‑1 (MQ=255)
                                                                          cTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATAc                                                                                        <  2:1613883/139‑1 (MQ=255)
                                                                          cTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATAc                                                                                        >  1:1613883/1‑139 (MQ=255)
                                                                            ttAGCTGTCGATATCAGCGACTATCTCGATGATGTCGGGCGGAGTAAGTGTGTGGTCTTCCTCCAGCATCGCGAGCAGGCGATCGTGTAACACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTgg                                                                           <  1:703822/149‑1 (MQ=255)
                                                                            gTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTgg                                                                           <  1:809161/150‑1 (MQ=255)
                                                                            gTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTgg                                                                           <  1:815079/150‑1 (MQ=255)
                                                                            gTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCACGCGATCGTGTAAAACTTCAACTTAAAGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTgg                                                                           >  1:253281/1‑150 (MQ=255)
                                                                             tAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCTGAGTAAGTGGCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGa                                                                          <  2:1619645/150‑1 (MQ=255)
                                                                                  gTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAg                                                                     <  1:476905/150‑1 (MQ=255)
                                                                                   tCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGc                                                                    <  2:1530633/150‑1 (MQ=255)
                                                                                    cGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATGGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCg                                                                   >  1:194384/1‑150 (MQ=255)
                                                                                    cGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCg                                                                   >  1:646376/1‑150 (MQ=255)
                                                                                      atatCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATAGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCggg                                                                 >  2:128503/1‑150 (MQ=255)
                                                                                       tatCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGc                                                                <  1:995494/150‑1 (MQ=255)
                                                                                           tgCCACCATCACGATGATGTCGCGCGTAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCATGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGttt                                                            <  2:381654/149‑1 (MQ=255)
                                                                                            gCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTg                                                           <  2:714215/150‑1 (MQ=255)
                                                                                            gCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTg                                                           <  2:814982/150‑1 (MQ=255)
                                                                                              caccaTCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGtt                                                         >  2:1195492/1‑150 (MQ=255)
                                                                                                 caTCACGATGATGTTGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATc                                                      <  2:612379/150‑1 (MQ=255)
                                                                                                 caTCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATc                                                      >  2:746700/1‑150 (MQ=255)
                                                                                                  aTCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCg                                                     <  2:1665765/150‑1 (MQ=255)
                                                                                                   tCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGc                                                    >  1:1265811/1‑150 (MQ=255)
                                                                                                    cACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCt                                                   >  2:1225640/1‑150 (MQ=255)
                                                                                                    cACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCt                                                   >  1:647927/1‑150 (MQ=255)
                                                                                                         tgatgTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACggg                                              >  1:1155797/1‑150 (MQ=255)
                                                                                                         tgatgTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACggg                                              <  1:1087331/150‑1 (MQ=255)
                                                                                                             gTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGaaaa                                          <  2:1498503/150‑1 (MQ=255)
                                                                                                               cgcgcgGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAAct                                        >  1:1026137/1‑150 (MQ=255)
                                                                                                                    ggAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCg                                   >  2:313101/1‑150 (MQ=255)
                                                                                                                     gAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGa                                  <  2:1448081/150‑1 (MQ=255)
                                                                                                                          aGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTc                             >  1:581063/1‑150 (MQ=255)
                                                                                                                               cGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACAcc                        <  1:1290769/150‑1 (MQ=255)
                                                                                                                                  gTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGc                     <  1:1372819/150‑1 (MQ=255)
                                                                                                                                  gTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGc                     >  2:425298/1‑150 (MQ=255)
                                                                                                                                  gTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGc                     >  2:120460/1‑150 (MQ=255)
                                                                                                                                       cctccAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAATCGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACgg                <  2:1090256/150‑1 (MQ=255)
                                                                                                                                       cctccAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACgg                >  1:1617970/1‑150 (MQ=255)
                                                                                                                                        ctccAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAAc                                         >  1:1315363/1‑124 (MQ=255)
                                                                                                                                        ctccAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAAc                                         <  2:1315363/124‑1 (MQ=255)
                                                                                                                                        ctccAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGgc               >  1:1095015/1‑150 (MQ=255)
                                                                                                                                         tacacgaTCGCCTGCTGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGgcg              >  2:114475/7‑150 (MQ=255)
                                                                                                                                           cacgaTCGCCTGCTGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGcgcg            >  2:569102/5‑150 (MQ=255)
                                                                                                                                                 cGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTttc      >  2:1007483/1‑150 (MQ=255)
                                                                                                                                                 cGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTttc      >  2:1013632/1‑150 (MQ=255)
                                                                                                                                                    cagcagGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTttccag   <  1:1195492/150‑1 (MQ=255)
                                                                                                                                                     agcagGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTttccagt  >  1:1600118/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTTTCCAGT  >  NC_000913/2961099‑2961397

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: