Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,087,396 T→G 8.0% F231V (TTC→GTC)  metK → S‑adenosylmethionine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,087,3960TG 8.0% 1349.2 / 23.5 528F231V (TTC→GTC) metKS‑adenosylmethionine synthetase
Reads supporting (aligned to +/- strand):  ref base T (183/290);  new base G (45/10);  total (228/300)
Fisher's exact test for biased strand distribution p-value = 7.56e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 569 total aligned reads displayed.

TCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGT  >  NC_000913/3087247‑3087545
                                                                                                                                                     |                                                                                                                                                     
   aTGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttcg                                                                                                                                                    <  2:5936683/150‑1 (MQ=255)
   aTGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttcg                                                                                                                                                    >  2:659475/1‑150 (MQ=255)
   aTGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttcg                                                                                                                                                    <  2:386079/150‑1 (MQ=255)
   aTGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttcg                                                                                                                                                    <  1:2723193/150‑1 (MQ=255)
       tGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTAt                                                                                                                                                <  2:8690518/150‑1 (MQ=255)
           gTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGt                                                                                                                                            <  1:5278541/150‑1 (MQ=255)
            tGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTg                                                                                                                                           <  2:4558354/150‑1 (MQ=255)
            tGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTg                                                                                                                                           <  1:8648551/150‑1 (MQ=255)
                ttCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGccc                                                                                                                                       >  2:4712326/1‑150 (MQ=255)
                ttCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGccc                                                                                                                                       <  2:7749412/150‑1 (MQ=255)
                   cACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAAt                                                                                                                                    >  1:2834904/1‑150 (MQ=255)
                   cACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATAGGTGGCCCAAt                                                                                                                                    >  2:6899900/1‑150 (MQ=255)
                    aCTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATg                                                                                                                                   >  1:4950702/1‑150 (MQ=255)
                     cTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGtcgtttc                                                                                                                                                     >  2:1386948/1‑131 (MQ=255)
                     cTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATgg                                                                                                                                  <  1:8955876/150‑1 (MQ=255)
                     cTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATgg                                                                                                                                  <  1:693511/150‑1 (MQ=255)
                      tCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATggg                                                                                                                                 <  1:515543/150‑1 (MQ=255)
                      tCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATggg                                                                                                                                 >  2:2875710/1‑150 (MQ=255)
                       cAGCACTCTTAATAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATAAATCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGt                                                                                                                                <  1:1686839/150‑1 (MQ=255)
                       cAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGt                                                                                                                                <  2:7069983/150‑1 (MQ=255)
                       cAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGt                                                                                                                                <  1:8458561/150‑1 (MQ=255)
                       cAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGt                                                                                                                                <  1:5574867/150‑1 (MQ=255)
                        aGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTg                                                                                                                               >  2:1060993/1‑150 (MQ=255)
                        aGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTg                                                                                                                               >  1:1817381/1‑150 (MQ=255)
                          cACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGGCCAATTGGTGAc                                                                                                                             >  1:4105905/1‑150 (MQ=255)
                          cACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGAc                                                                                                                             >  2:5320821/1‑150 (MQ=255)
                           aCTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACt                                                                                                                            <  1:4622145/150‑1 (MQ=255)
                            ctctGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAAGTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTAGCGGTGGCCCAATGGGTGACTg                                                                                                                           <  1:792283/150‑1 (MQ=255)
                               tGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCgg                                                                                                                        >  1:8005533/1‑150 (MQ=255)
                                gAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGt                                                                                                                       >  1:101092/1‑150 (MQ=255)
                                 aaGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTc                                                                                                                      <  1:3588043/150‑1 (MQ=255)
                                    agaTCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGa                                                                                                                   <  2:4186415/150‑1 (MQ=255)
                                       tCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATTGGTGACTGCGGTCTGACTg                                                                                                                <  1:9158314/150‑1 (MQ=255)
                                       tCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTg                                                                                                                <  1:6669247/150‑1 (MQ=255)
                                          aCCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGTCCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTc                                                                                                             <  1:3115262/150‑1 (MQ=255)
                                          aCCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTc                                                                                                             <  1:3747153/150‑1 (MQ=255)
                                             aGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTa                                                                                                          >  2:256220/1‑150 (MQ=255)
                                             aGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTa                                                                                                          >  2:8077470/1‑150 (MQ=255)
                                              gAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTaa                                                                                                         >  2:793472/1‑150 (MQ=255)
                                               aaaTCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGTTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTaaa                                                                                                        <  1:1252227/150‑1 (MQ=255)
                                                aaTCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTaaaa                                                                                                       <  1:6396553/150‑1 (MQ=255)
                                                    gCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTAt                                                                                                   <  1:6903999/150‑1 (MQ=255)
                                                    gCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCACGCTGAAAGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGAACGGTCGTTTCGTTATCGGTTGCCCAATGGGTGAATGAGGTCTGACTGGTCGTAAAATTAt                                                                                                   >  2:6760361/1‑150 (MQ=255)
                                                        cAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGtt                                                                                               <  2:5556286/150‑1 (MQ=255)
                                                         aagaagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAAGGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTg                                                                                              >  2:3611302/1‑150 (MQ=255)
                                                           gaagCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGAt                                                                                            <  2:5156235/150‑1 (MQ=255)
                                                                ggTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTa                                                                                       <  1:7351620/150‑1 (MQ=255)
                                                                 gTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTAc                                                                                      <  1:7626616/150‑1 (MQ=255)
                                                                  tAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTaac                                                                                     >  1:3374908/1‑148 (MQ=255)
                                                                   aaTGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTAcgg                                                                                    >  2:7643203/1‑150 (MQ=255)
                                                                       gAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCCAAGGGTGACTGCGGGCTGACTGGTCGTAAAATTATCGGTGATACCAACggcggc                                                                                >  1:2413831/1‑150 (MQ=255)
                                                                       gAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACggcggc                                                                                >  1:6056392/1‑150 (MQ=255)
                                                                       gAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACggcggc                                                                                <  1:7656977/150‑1 (MQ=255)
                                                                       gAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACggcggc                                                                                >  1:7449706/1‑150 (MQ=255)
                                                                       gAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCCAGGGGTGACGGCGGGCTGACGGGTCGGAAAATTATCGTTGATACCGACggcggc                                                                                >  1:6337316/1‑150 (MQ=255)
                                                                        aaGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCa                                                                               <  1:4727944/150‑1 (MQ=255)
                                                                          gagaTCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATAGGTGGCCCAAAGGGTGACTGCGGTCTGACTGGTCGTAAAATTATAGTTGATACCTACGGCGGCATg                                                                             >  1:4226234/1‑150 (MQ=255)
                                                                           agaTCATCAAGCCAATTCTGCCCGCTGAATGGCTGTCTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGTTACCTACGGCGGCATgg                                                                            <  2:2760323/150‑1 (MQ=255)
                                                                           agaTCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATgg                                                                            <  2:3256910/150‑1 (MQ=255)
                                                                            gaTCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGgc                                                                           <  2:8235535/150‑1 (MQ=255)
                                                                            gaTCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGgc                                                                           >  2:9148820/1‑150 (MQ=255)
                                                                             atcatcAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGgcg                                                                          >  2:1536441/1‑150 (MQ=255)
                                                                             atcatcAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGgcg                                                                          <  2:8834261/150‑1 (MQ=255)
                                                                                    aGCCAATTCGGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGGTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACg                                                                   <  1:7794983/150‑1 (MQ=255)
                                                                                      ccAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGt                                                                 >  2:7257197/1‑150 (MQ=255)
                                                                                             ggCCCGCTGAAGGGCTGACGTCGGCCCCCACCGTCGTCATCACCCCGCCCGGGCGGGTCGTTATCGGTGGCCCCCTGGGGGACTGCGGTCTGACTGGTCGGCACATTATCGTGGCTACCTCCGGCGGCATGGCGCGTCACGGGGGCGGTg                                                          <  2:737897/149‑1 (MQ=255)
                                                                                               cccGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCa                                                        >  2:9049784/1‑150 (MQ=255)
                                                                                               cccGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCa                                                        <  2:403482/150‑1 (MQ=255)
                                                                                               cccGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTGGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATAGTTGATAACTACCGCGGCAAGGCGCGTCACGGTGGCGGTGCa                                                        >  1:444322/1‑150 (MQ=255)
                                                                                               cccGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCAACGGCGGCATGGCGCGTCACGGTGGCGGTGCa                                                        >  2:7830826/1‑150 (MQ=255)
                                                                                                  gCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATtc                                                     <  1:7270878/150‑1 (MQ=255)
                                                                                                    tGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATtctc                                                   <  1:1556716/150‑1 (MQ=255)
                                                                                                    tGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATtctc                                                   <  1:8211520/150‑1 (MQ=255)
                                                                                                     gAATGGCTGACTTCTGACACCAAATTCTTCATCAACCCGACCGGTCGTGTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTctct                                                  <  1:6264554/150‑1 (MQ=255)
                                                                                                      aaTGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAACTTATCGTTGATACCTCCGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTg                                                 <  2:223002/150‑1 (MQ=255)
                                                                                                      aaTGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTg                                                 <  2:8426841/150‑1 (MQ=255)
                                                                                                           cTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTaaa                                            <  2:6521234/150‑1 (MQ=255)
                                                                                                            tGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATGGGTGGGCCAATGGGTGTCTGGGGTGTGAGTGGTGGTAAAATTATGGTTGTTACCTACGGCGGCATGGGGCGTCACGGTGGCGGTGGATTGTGTTGTaaat                                           >  2:4337300/1‑149 (MQ=255)
                                                                                                            tGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAg                                           >  2:4898791/1‑150 (MQ=255)
                                                                                                             gACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGa                                          <  2:1181058/150‑1 (MQ=255)
                                                                                                             gACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGa                                          <  1:2052826/150‑1 (MQ=255)
                                                                                                                       accaAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaaa                                <  1:7120446/150‑1 (MQ=255)
                                                                                                                       accaAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATACAACaaaa                                >  2:1760514/1‑150 (MQ=255)
                                                                                                                         caAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGt                              <  1:4334842/150‑1 (MQ=255)
                                                                                                                         caAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAAAGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATAACTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGt                              >  2:8909461/1‑150 (MQ=255)
                                                                                                                            aTTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGa                           <  2:5762044/150‑1 (MQ=255)
                                                                                                                             ttcttcATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGAc                          >  2:2029021/1‑150 (MQ=255)
                                                                                                                             ttcttcATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGAc                          <  1:3451614/150‑1 (MQ=255)
                                                                                                                                   atccACCCGACCGGTCGGGTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTcc                    <  1:1088235/150‑1 (MQ=255)
                                                                                                                                   atcaACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAACGTGGACCGTTcc                    <  1:590853/150‑1 (MQ=255)
                                                                                                                                    tcaACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCg                   >  2:2117230/1‑150 (MQ=255)
                                                                                                                                         ccGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGcc              <  2:6345134/150‑1 (MQ=255)
                                                                                                                                           gACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGATTGGTCGTAAAATTATCGTTGATTCCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTa            <  1:2177094/150‑1 (MQ=255)
                                                                                                                                             ccGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACg          <  1:5676469/150‑1 (MQ=255)
                                                                                                                                             ccGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACg          <  1:8102386/150‑1 (MQ=255)
                                                                                                                                               ggTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACgca        <  2:6165884/150‑1 (MQ=255)
                                                                                                                                                gTCGTTTCGTTATCGGTTGCCCAATGGGTTACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATTGCGCGTCACGGTGGCGGTGCATTATATGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACgcag       >  2:345406/1‑150 (MQ=255)
                                                                                                                                                    tttcgttATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACg   >  1:3919591/1‑150 (MQ=255)
                                                                                                                                                    tttcgttATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACg   >  2:3117882/1‑150 (MQ=255)
                                                                                                                                                     ttcgttATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGt  >  1:5385652/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGATGCTGTCGTGCTTTCCACTCAGCACTCTGAAGAGATCGACCAGAAATCGCTGCAAGAAGCGGTAATGGAAGAGATCATCAAGCCAATTCTGCCCGCTGAATGGCTGACTTCTGCCACCAAATTCTTCATCAACCCGACCGGTCGTTTCGTTATCGGTGGCCCAATGGGTGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGT  >  NC_000913/3087247‑3087545

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