Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 22,773 A→T 8.1% E128V (GAA→GTA)  ileS → isoleucyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091322,7730AT 8.1% 391.7 / 12.0 117E128V (GAA→GTA) ileSisoleucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (52/55);  new base T (3/7);  total (55/62)
Fisher's exact test for biased strand distribution p-value = 3.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

TGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAA  >  NC_000913/22626‑22914
                                                                                                                                                   |                                                                                                                                             
tGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCATCGCCGCCGAGTTCCGCGCCAAGTGCCGCGaa                                                                                                                                              <  4:662328/149‑1 (MQ=255)
tGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGaa                                                                                                                                              <  1:317960/149‑1 (MQ=255)
tGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGaa                                                                                                                                              >  4:521487/1‑149 (MQ=255)
    gTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATAcg                                                                                                                                          >  1:1047722/1‑149 (MQ=255)
     tCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATAcgc                                                                                                                                         >  3:896736/1‑149 (MQ=255)
       cAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACgcgg                                                                                                                                       <  2:420925/149‑1 (MQ=255)
       cAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACgcgg                                                                                                                                       <  2:599492/149‑1 (MQ=255)
        aaaGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACgcggc                                                                                                                                      >  1:399238/1‑149 (MQ=255)
          aGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGTATTCgcggcac                                                                                                                                    <  2:430692/149‑3 (MQ=255)
          aGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGa                                                                                                                                    <  1:445012/149‑1 (MQ=255)
          aGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGa                                                                                                                                    <  3:309917/149‑1 (MQ=255)
           gggCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGAc                                                                                                                                   >  4:672506/1‑149 (MQ=255)
            ggCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGTATTCgcggcactt                                                                                                                                  <  3:223973/149‑5 (MQ=255)
            ggCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGAcc                                                                                                                                  >  2:62062/1‑149 (MQ=255)
             gCTTTCCGGTTATGACTCGCCGTATGTGCCTGGGTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGAccc                                                                                                                                 <  4:714254/149‑1 (MQ=255)
             gCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCTCGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGAccc                                                                                                                                 <  3:924352/149‑1 (MQ=255)
             gCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGTATTCgcggcacttg                                                                                                                                 <  2:581626/149‑6 (MQ=255)
             gCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGAccc                                                                                                                                 >  4:813163/1‑149 (MQ=255)
               tttCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAg                                                                                                                               >  3:336018/1‑149 (MQ=255)
                    ggTTATGACTCGCCGTATGTGCCTGGCTTGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGa                                                                                                                          <  2:136155/149‑1 (MQ=255)
                      ttATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACg                                                                                                                        >  4:281630/1‑149 (MQ=255)
                          gACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCa                                                                                                                    >  1:1079919/1‑149 (MQ=255)
                          gACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCa                                                                                                                    >  1:65242/1‑149 (MQ=255)
                            cTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGTATTCgcggcacttgg                                                                                                                                <  4:793241/135‑7 (MQ=255)
                            cTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGTATTCgcggcacttgg                                                                                                                                >  3:793241/1‑129 (MQ=255)
                             tCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACg                                                                                                                 >  4:437282/1‑149 (MQ=255)
                               gCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGGGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCTCGAATACGCGGCGACCCAGGTTGACGGTCAACGCa                                                                                                               >  4:704691/1‑149 (MQ=255)
                                ccGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCaa                                                                                                              <  1:995071/149‑1 (MQ=255)
                                 cGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCaaa                                                                                                             <  4:278429/149‑1 (MQ=255)
                                      gtgCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACtt                                                                                                        >  3:41909/1‑149 (MQ=255)
                                       tgCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACttt                                                                                                       >  4:1055757/1‑149 (MQ=255)
                                          ctggctggGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGGGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCCGGTTGACGGTCAACGCAAAGACTTTATc                                                                                                    >  1:176294/1‑149 (MQ=255)
                                          ctggctggGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGGGAGAAATTCACCGCCGCCGAGGTCAGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATc                                                                                                    >  2:679975/1‑149 (MQ=255)
                                           tggctggGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATcc                                                                                                   >  2:272130/1‑149 (MQ=255)
                                              ctggGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTc                                                                                                >  1:262756/1‑149 (MQ=255)
                                              ctggGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTc                                                                                                >  3:278732/1‑149 (MQ=255)
                                                gggACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGATGGTCAACGCAAAGACTTTATCCGTCTg                                                                                              <  4:788449/149‑1 (MQ=255)
                                                   aCTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGc                                                                                           <  2:109321/149‑1 (MQ=255)
                                                      gCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTg                                                                                        <  1:425548/149‑1 (MQ=255)
                                                         aCGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTTACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTg                                                                                     >  4:160185/1‑149 (MQ=255)
                                                          cGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTgg                                                                                    <  4:517930/149‑1 (MQ=255)
                                                          cGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTgg                                                                                    <  1:465887/149‑1 (MQ=255)
                                                            gTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGc                                                                                  >  1:251250/1‑149 (MQ=255)
                                                               tGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGAc                                                                               >  4:425263/1‑149 (MQ=255)
                                                               tGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGAc                                                                               <  4:336018/149‑1 (MQ=255)
                                                                  cGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTgg                                                                            <  1:822191/149‑1 (MQ=255)
                                                                  cGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTgg                                                                            >  2:568953/1‑149 (MQ=255)
                                                                   gATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGt                                                                           >  2:83210/1‑149 (MQ=255)
                                                                        aGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCAc                                                                      <  1:336343/149‑1 (MQ=255)
                                                                        aGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGGTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCAc                                                                      <  4:798210/149‑1 (MQ=255)
                                                                          cTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCAccc                                                                    >  1:584623/1‑149 (MQ=255)
                                                                            gAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTa                                                                                                      <  1:692634/113‑1 (MQ=255)
                                                                              aaGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTAc                                                                <  1:272130/149‑1 (MQ=255)
                                                                               aGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTAcc                                                               >  2:399830/1‑149 (MQ=255)
                                                                                 tCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTg                                                             >  4:809132/1‑149 (MQ=255)
                                                                                   gAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGAc                                                           >  2:1000686/1‑149 (MQ=255)
                                                                                   gAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGAc                                                           >  1:680096/1‑149 (MQ=255)
                                                                                    aGCAAGAATACGGTAAGCCGGGTGAGAAATTCACTGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGAcc                                                          >  3:50546/1‑149 (MQ=255)
                                                                                    aGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGa                                                                                                           >  4:536686/1‑100 (MQ=255)
                                                                                    aGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGa                                                                                                           <  3:536686/100‑1 (MQ=255)
                                                                                    aGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGAcc                                                          >  2:822998/1‑149 (MQ=255)
                                                                                     gCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCa                                                         >  2:224742/1‑149 (MQ=255)
                                                                                     gCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCa                                                         >  4:16763/1‑149 (MQ=255)
                                                                                       aaGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATg                                                       >  2:272510/1‑149 (MQ=255)
                                                                                              cGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCa                                                <  2:65242/149‑1 (MQ=255)
                                                                                               ggTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCaa                                               >  1:579876/1‑149 (MQ=255)
                                                                                                 tAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCaaaa                                             <  1:568953/149‑1 (MQ=255)
                                                                                                 tAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCaaaa                                             <  3:821243/149‑1 (MQ=255)
                                                                                                      cGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGaa                                        <  2:842185/149‑1 (MQ=255)
                                                                                                        ggTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGc                                      <  3:336172/149‑1 (MQ=255)
                                                                                                          tGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCa                                    >  4:1080222/1‑149 (MQ=255)
                                                                                                           gagaAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCaa                                   <  1:793460/149‑1 (MQ=255)
                                                                                                            agaAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCACAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACGGAAGCCAAc                                  <  2:176294/149‑1 (MQ=255)
                                                                                                            agaAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAAc                                  <  1:628248/149‑1 (MQ=255)
                                                                                                             gaAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAAca                                 >  3:985811/1‑149 (MQ=255)
                                                                                                              aaaTTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAAcat                                >  4:484832/1‑149 (MQ=255)
                                                                                                               aaTTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAAcatc                               >  4:736366/1‑149 (MQ=255)
                                                                                                               aaTTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAAcatc                               >  4:551438/1‑149 (MQ=255)
                                                                                                                 ttCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAAcatcat                             >  1:177768/1‑149 (MQ=255)
                                                                                                                 ttCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAAcatcat                             >  3:390391/1‑149 (MQ=255)
                                                                                                                  tCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACatcatc                            <  2:86203/149‑1 (MQ=255)
                                                                                                                   cACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACttt                                                                                                       >  1:456255/1‑73 (MQ=255)
                                                                                                                   cACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACttt                                                                                                       <  2:456255/73‑1 (MQ=255)
                                                                                                                         cgccgAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCgcgcgc                     <  2:253067/149‑1 (MQ=255)
                                                                                                                         cgccgAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTAAACGCAAAGACTTTATCCGGCTGGGCGTGCTGCGCGTCTTGTCGCACCCGTACGTGACCATGGACTTCAAAACTGAAGCCAACATCATCCgcgcgc                     <  2:293247/149‑1 (MQ=255)
                                                                                                                           ccgAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTg                   <  1:922349/149‑1 (MQ=255)
                                                                                                                            cgAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTgg                  <  1:495526/149‑1 (MQ=255)
                                                                                                                            cgAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTgg                  <  2:115926/149‑1 (MQ=255)
                                                                                                                            cgAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTgg                  <  3:76459/149‑1 (MQ=255)
                                                                                                                             gAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTggg                 <  1:224742/149‑1 (MQ=255)
                                                                                                                                 tCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCaaa             <  1:399830/149‑1 (MQ=255)
                                                                                                                                  ccGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCaaaa            <  3:281630/149‑1 (MQ=255)
                                                                                                                                      gcCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAatcat        <  3:494867/149‑1 (MQ=255)
                                                                                                                                      gcCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAatcat        <  3:218627/149‑1 (MQ=255)
                                                                                                                                       ccAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGCCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAatcatc       <  4:1067178/149‑1 (MQ=255)
                                                                                                                                         aaGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCgg     <  3:869460/149‑1 (MQ=255)
                                                                                                                                         aaGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCgg     <  2:257488/149‑1 (MQ=255)
                                                                                                                                            cgccgccTATTCGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCaa  >  2:779221/2‑149 (MQ=255)
                                                                                                                                            cgCCGTGTATTCGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAatcat        <  2:924876/142‑1 (MQ=255)
                                                                                                                                            cgCCGCGTATTCGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAatcat        >  1:924876/2‑143 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
TGAAGTCCAAAGGGCTTTCCGGTTATGACTCGCCGTATGTGCCTGGCTGGGACTGCCACGGTCTGCCGATCGAGCTGAAAGTCGAGCAAGAATACGGTAAGCCGGGTGAGAAATTCACCGCCGCCGAGTTCCGCGCCAAGTGCCGCGAATACGCGGCGACCCAGGTTGACGGTCAACGCAAAGACTTTATCCGTCTGGGCGTGCTGGGCGACTGGTCGCACCCGTACCTGACCATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAA  >  NC_000913/22626‑22914

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