Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,879,197 T→G E192D (GAA→GAC ygbT ← multifunctional endonuclease Cas1, CRISPR adaptation protein; DNA repair enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,879,1970TG97.5% 296.9 / ‑0.7 79E192D (GAA→GACygbTmultifunctional endonuclease Cas1, CRISPR adaptation protein; DNA repair enzyme
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base G (43/34);  total (45/34)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

ACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGAAGTGCGTAGGTTGCCCGCACGC  >  NC_000913/2879050‑2879342
                                                                                                                                                   |                                                                                                                                                 
aCCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCa                                                                                                                                                <  2:645809/151‑1 (MQ=255)
  cAGGGTTACGACGCGCTCTCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCa                                                                                                                                                <  2:952233/149‑1 (MQ=255)
    gggTTACGACGCGCTAGCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAACCCCAATAGCTGGTGCATAACCAGCTGCAGGTACCGCCGCGTCAGTTa                                                                                                                                            <  2:250902/151‑1 (MQ=255)
     ggTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTAc                                                                                                                                           <  2:826014/151‑1 (MQ=255)
       ttACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGc                                                                                                                                         >  1:84396/1‑151 (MQ=255)
         acgacgCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGcc                                                                                                                                        <  2:337620/150‑1 (MQ=255)
               cgcTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATaaa                                                                                                                                 <  2:1111196/151‑1 (MQ=255)
                   aTCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAgg                                                                                                                             >  1:583305/1‑151 (MQ=255)
                    tctcAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAgg                                                                                                                             >  1:728909/1‑150 (MQ=255)
                        aaaaGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGtt                                                                                                                        <  2:958624/151‑1 (MQ=255)
                             cTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAgc                                                                                                                   >  1:747573/1‑151 (MQ=255)
                             cTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAgc                                                                                                                   >  2:532285/1‑151 (MQ=255)
                                   gTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAAt                                                                                                             <  2:1166940/151‑1 (MQ=255)
                                    tACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATg                                                                                                            >  2:185866/1‑151 (MQ=255)
                                      caacaGTGACAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGc                                                                                                           >  2:253656/1‑150 (MQ=255)
                                        acaGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCAtt                                                                                                        <  2:224947/151‑1 (MQ=255)
                                               cAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGAt                                                                                                 >  1:319445/1‑151 (MQ=255)
                                                aaaTTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATc                                                                                                >  1:781036/1‑151 (MQ=255)
                                                aaaTTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATc                                                                                                >  1:282048/1‑151 (MQ=255)
                                                  aTTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGt                                                                                              <  2:984151/151‑1 (MQ=255)
                                                   tttAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGt                                                                                              >  1:913571/1‑150 (MQ=255)
                                                    ttAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTa                                                                                             <  1:409196/150‑1 (MQ=255)
                                                    ttAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTAt                                                                                            >  1:625044/1‑151 (MQ=255)
                                                     tAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCGAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATc                                                                                           >  1:842754/1‑151 (MQ=255)
                                                     tAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTCCGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTa                                                                                             <  2:732521/149‑1 (MQ=255)
                                                     tAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATCCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGCAGTTGCAGCGCTAATGCATTGGTTGATCGTAt                                                                                            >  2:983451/1‑150 (MQ=255)
                                                      aaTGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCg                                                                                          <  1:553279/151‑1 (MQ=255)
                                                      aaTGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCg                                                                                          >  2:717505/1‑151 (MQ=255)
                                                       atgatgTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGc                                                                                         >  1:117633/1‑151 (MQ=255)
                                                              cTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTc                                                                                  >  2:952683/1‑151 (MQ=255)
                                                               tGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTcc                                                                                 >  2:335698/1‑151 (MQ=255)
                                                                gCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTccc                                                                                >  2:336026/1‑151 (MQ=255)
                                                                     atCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTc                                                                            <  1:382500/150‑1 (MQ=255)
                                                                     atCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTc                                                                            <  1:937565/150‑1 (MQ=255)
                                                                        gTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTc                                                                        <  2:817698/151‑1 (MQ=255)
                                                                          aaaCAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCgg                                                                      >  2:444907/1‑151 (MQ=255)
                                                                          aaaCAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCg                                                                       <  1:1042663/150‑1 (MQ=255)
                                                                             cAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATc                                                                   >  2:656294/1‑151 (MQ=255)
                                                                             cAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATc                                                                   >  2:69506/1‑151 (MQ=255)
                                                                              aaaGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCg                                                                  >  2:3692/1‑151 (MQ=255)
                                                                               aaGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGt                                                                 <  1:1149808/151‑1 (MQ=255)
                                                                               aaGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCCCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCTCAGTCTTTCGGATCg                                                                  >  2:1192361/1‑150 (MQ=255)
                                                                                  gAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGc                                                              >  1:1100621/1‑151 (MQ=255)
                                                                                  gAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGc                                                              >  2:374492/1‑151 (MQ=255)
                                                                                  gAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTAACGCCCTTTTCCCAGTCTTTCGGATCGTAGc                                                              >  1:677342/1‑151 (MQ=255)
                                                                                   aaaGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCg                                                             <  2:435445/151‑1 (MQ=255)
                                                                                     agagGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGAc                                                           <  1:91072/151‑1 (MQ=255)
                                                                                     agagGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGAc                                                           <  2:921874/151‑1 (MQ=255)
                                                                                     agagGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGAc                                                           >  2:894405/1‑151 (MQ=255)
                                                                                     agagGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGAc                                                           >  1:539410/1‑151 (MQ=255)
                                                                                       agGCTTGCCTGTATGCACAAACCCAATAGCTGATGCATAACTAGCTGCAAGTATCGCCACGTCAGTTACGCCGTATAAACAGGATGTTGCAGCGGTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGAGCGTAGCTAc                                                           <  2:741762/149‑1 (MQ=255)
                                                                                           tttCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGtccat                                                     >  1:277335/1‑151 (MQ=255)
                                                                                              ccTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGtccattc                                                   <  2:891608/150‑1 (MQ=255)
                                                                                               cTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGtccattcc                                                  >  1:1088465/1‑150 (MQ=255)
                                                                                                tGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGtccattccat                                                >  1:270181/1‑151 (MQ=255)
                                                                                                 gTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGtccattccat                                                <  1:532310/150‑1 (MQ=255)
                                                                                                  tATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGTCGTCCATTCCATg                                               <  2:1057970/150‑1 (MQ=255)
                                                                                                      cacaAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACg                                          <  1:309267/151‑1 (MQ=255)
                                                                                                      cacaAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACg                                          <  1:1076915/151‑1 (MQ=255)
                                                                                                        caAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGcc                                        <  1:1005177/151‑1 (MQ=255)
                                                                                                                  aGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCtt                               <  1:836231/150‑1 (MQ=255)
                                                                                                                  aGCTGGTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCtt                               <  2:953760/150‑1 (MQ=255)
                                                                                                                       gTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGGTGATCGTATCGTCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCa                          >  2:1110513/1‑150 (MQ=255)
                                                                                                                       gTGCATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGATGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGcc                           >  2:728835/1‑149 (MQ=255)
                                                                                                                          cATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGa                        >  2:905810/1‑149 (MQ=255)
                                                                                                                          cATAACCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGa                        >  2:231454/1‑149 (MQ=255)
                                                                                                                             aaCCAGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGAAGTGc                   >  2:128021/1‑151 (MQ=255)
                                                                                                                                 aGCTGCAAGTATCGCCTCGTCAGTTAGGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGTTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGAAGTGCGt                 >  2:947335/1‑149 (MQ=255)
                                                                                                                                 aGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGAAGTGCGt                 >  1:167100/1‑149 (MQ=255)
                                                                                                                                 aGCTGCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGAAGTGCGt                 >  1:1211776/1‑149 (MQ=255)
                                                                                                                                  gCTGCAAGTATCGCCGCGTCAGTTACGCCGCATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGAAGTGCGTAgg              <  1:680414/151‑1 (MQ=255)
                                                                                                                                    tGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGAAGTGCGTAGGtt            >  2:124008/1‑151 (MQ=255)
                                                                                                                                     gCAAGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGAAGTGCGTAGGTTg           <  2:500470/151‑1 (MQ=255)
                                                                                                                                       aaGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGAAGTGCGTAGGTTGcc         >  2:152098/1‑151 (MQ=255)
                                                                                                                                       aaGTATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGAAGTGCGTAGGTTGcc         >  2:676478/1‑151 (MQ=255)
                                                                                                                                          tATCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGAAGTGCGTAGGTTGCCCGc      >  1:841330/1‑151 (MQ=255)
                                                                                                                                            tCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGAAGTGCGTAGGTTGCCCGCa     <  1:529693/150‑1 (MQ=255)
                                                                                                                                            tCGCCGCGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTAATCGTATCGCCCTTTTCCTAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGCATTGCTTCGCCAGAAGTGCGTAGGTTGCCCGCAc    >  2:306409/1‑151 (MQ=255)
                                                                                                                                               ccgcGTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGAAGTGCGTAGGTTGCCCGCAcgc  <  1:656325/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ACCAGGGTTACGACGCGCTATCTCAAAAGCTTTCGGTACAACAGTGTCAAATTTAATGATGTCTGCAATATCGTAAACAAAGGAAAGAGGCTTTCCTGTATGCACAAACCCAATAGCTGGTGCATAACCAGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGAAGTGCGTAGGTTGCCCGCACGC  >  NC_000913/2879050‑2879342

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