Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,722,928 Δ1 bp coding (467/1527 nt) PA14_RS16920 → inactive transglutaminase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,722,9280A.89.9% 1976.5 / 243.1 487coding (467/1527 nt)PA14_RS16920inactive transglutaminase family protein
Reads supporting (aligned to +/- strand):  ref base A (24/25);  new base . (191/247);  total (215/272)
Fisher's exact test for biased strand distribution p-value = 5.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 487 total aligned reads displayed.

CAGACCCTGTACTACCGCCTGGTGCTGACCAAGCGCTACACCGGCGAGAAGACCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGACCCTGCTGTCCATCGCCCACG  >  NC_008463/3722783‑3723074
                                                                                                                                                 |                                                                                                                                                  
 aGACCCTGTACTACCGCCTGGTGCTGACCAAGCGCTACACCGGCGAGAAGACCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCgg                                                                                                                                               <  1:2356942/149‑1 (MQ=255)
  gACCCTGTACTACCGCCTGGTGCTGACCAAGCGCTACACCGGCGAGAAGACCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGcc                                                                                                                                             >  1:3089137/1‑150 (MQ=255)
  gACCCTGTACTACCGCCTGGTGCTGACCAAGCGCTACACCGGCGAGAAGACCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGcc                                                                                                                                             >  2:1071708/1‑150 (MQ=255)
    cccTGTACTACCGCCTGGTGCTGACCAAGCGCTACACCGGCGAGAAGACCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCACTCGGCCg                                                                                                                                            <  1:7963964/150‑1 (MQ=255)
       tGTACTACCGCCTGGTGCTGACCAAGCGCTACACCGGCGAGAAGACCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGt                                                                                                                                        >  2:10950116/1‑150 (MQ=255)
        gTACTACCGCCTGGTGCTGACCAAGCGCTACACCGGCGAGAAGACCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTc                                                                                                                                       >  2:8825698/1‑150 (MQ=255)
        gTACTACCGCCTGGTGCTGACCAAGCGCTACACCGGCGAGAAGACCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTc                                                                                                                                       <  1:8344409/150‑1 (MQ=255)
        gTACGACCGCCTGGTGCTGACCAAGCGCTACACCGGCGAGAAGACCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTc                                                                                                                                       <  2:7571060/150‑1 (MQ=255)
          actacCCCCCGGTGCTGACCAAGCGCGACCCCGGCGCGAAGCCCAAGGAGAAAGGTCCGAACGGCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCgc                                                                                                                                     <  2:4392946/150‑2 (MQ=255)
              ccGCCTGGTGCTGACCAAGCGCTACACCGGCGAGAAGACCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGAcc                                                                                                                                 <  2:10645383/150‑1 (MQ=255)
                gCCTGGAGCTGACCAAGCGCTACACCGGCGAGAAGACCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCACTCGGCCGACGTCGAGACCt                                                                                                                                <  2:9798403/150‑1 (MQ=255)
                      tGCTGACCAAGCGCTACACCGGCGAGAAGACCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAg                                                                                                                         <  2:3235830/150‑1 (MQ=255)
                       gCTGACCAAGCGCTACACCGGCGAGAAGACCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGc                                                                                                                        <  1:3602281/150‑1 (MQ=255)
                         tGACCAAGCGCTACACCGGCGAGAAGACCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGa                                                                                                                      <  2:8432820/150‑1 (MQ=255)
                             cAAGCGCTACACCGGCGAGAAGACCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGc                                                                                                                   >  2:1476230/1‑150 (MQ=255)
                                 cgcTACCCCCGCGAGAAGACCAAGGAGACAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCAGTCGCCAGCCTCGGCCGCCGTCGAGACCTTCGTCAGCGAAGCCATCa                                                                                                              <  1:817096/150‑1 (MQ=255)
                                  gcTACACCGGCGAGAAGACCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCaa                                                                                                             <  1:7366185/150‑1 (MQ=255)
                                   cTACACCGGCGAGAAGACCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAg                                                                                                            <  1:7255341/150‑1 (MQ=255)
                                    tACACCGGCGAGAAGACCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAgc                                                                                                           >  1:5599049/1‑150 (MQ=255)
                                        ccGGCGAGAAGACCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAgcgcgc                                                                                                       >  2:6499805/1‑150 (MQ=255)
                                        ccGGCGAGAAGACCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAgcgcgc                                                                                                       >  1:6043831/1‑150 (MQ=255)
                                          ggCGAGAAGACCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCca                                                                                                     >  1:4354198/1‑150 (MQ=255)
                                              agaagaCCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCcaacaa                                                                                                 >  2:6654240/1‑150 (MQ=255)
                                              agaagaCCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCcaacaa                                                                                                 <  2:8452677/150‑1 (MQ=255)
                                               gaagaCCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCaacaac                                                                                                >  2:10797882/1‑150 (MQ=255)
                                               gaagaCCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCaacaac                                                                                                <  2:3851039/150‑1 (MQ=255)
                                                 atacaaaGGATAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGTAAGGCCCGGAGAAGAGAGCCGCAGAAGAACTGCTGGCCCCCATTCGACATCCTCGGCCGACGTCGAGACCAACGTAAGCGAAGCCAGCAACCGCGCCAACAACgg                                                                                              >  2:1777771/6‑150 (MQ=255)
                                                 agaCCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACgg                                                                                              >  1:8983822/1‑150 (MQ=255)
                                                 agaCCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACgg                                                                                              <  1:9966180/150‑1 (MQ=255)
                                                    ccAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCaa                                                                                           <  1:9516572/150‑1 (MQ=255)
                                                    ccAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCACCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCaa                                                                                           >  1:10896868/1‑150 (MQ=255)
                                                     cAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAac                                                                                          <  1:7056069/150‑1 (MQ=255)
                                                       aGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAacga                                                                                        >  2:3042709/1‑150 (MQ=255)
                                                        ggAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAacgac                                                                                       >  2:548457/1‑150 (MQ=255)
                                                        ggAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAacgac                                                                                       >  2:9450989/1‑150 (MQ=255)
                                                          agaAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAAcgacga                                                                                     >  2:8163342/1‑150 (MQ=255)
                                                          agaAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAAcgacga                                                                                     >  1:527144/1‑150 (MQ=255)
                                                           gaAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTGGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACgacgac                                                                                    <  1:9617336/150‑1 (MQ=255)
                                                           gaAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACgacgac                                                                                    <  1:3395325/150‑1 (MQ=255)
                                                             aaGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACaa                                                                                  >  1:2037800/1‑150 (MQ=255)
                                                             aaGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACaa                                                                                  <  1:10502759/150‑1 (MQ=255)
                                                              aGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAAc                                                                                 <  2:3997400/150‑1 (MQ=255)
                                                               ggTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACg                                                                                >  2:5262762/1‑150 (MQ=255)
                                                                gTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGt                                                                               <  1:7714931/150‑1 (MQ=255)
                                                                 tCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGt                                                                               <  1:10147703/150‑1 (MQ=255)
                                                                    gATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAg                                                                           <  1:5726759/150‑1 (MQ=255)
                                                                     aTCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAg                                                                           <  1:4564065/150‑1 (MQ=255)
                                                                     aTCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAgc                                                                          >  1:2965643/1‑150 (MQ=255)
                                                                     aTCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAgc                                                                          >  2:1248339/1‑150 (MQ=255)
                                                                     aTCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAgc                                                                          <  2:501188/150‑1 (MQ=255)
                                                                     aTCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAgc                                                                          <  1:6458064/150‑1 (MQ=255)
                                                                       cTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAgct                                                                         <  1:6938995/150‑1 (MQ=255)
                                                                       cTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAgctg                                                                        <  1:8593851/150‑1 (MQ=255)
                                                                        ttCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAgctgc                                                                       >  2:6291989/1‑150 (MQ=255)
                                                                           cgcgACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGa                                                                    <  2:347879/150‑1 (MQ=255)
                                                                            gcgACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGAt                                                                   >  2:3770461/1‑150 (MQ=255)
                                                                             cgACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATg                                                                  >  2:4568103/1‑150 (MQ=255)
                                                                              gACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATgg                                                                 >  2:4872989/1‑150 (MQ=255)
                                                                               aCAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGc                                                                >  2:903573/1‑150 (MQ=255)
                                                                               aCAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGc                                                                <  2:417781/150‑1 (MQ=255)
                                                                                cAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGcc                                                               <  1:8702540/150‑1 (MQ=255)
                                                                                  gCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCgg                                                             <  1:2677341/150‑1 (MQ=255)
                                                                                  gCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCgg                                                             <  2:2331190/150‑1 (MQ=255)
                                                                                  gCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCgg                                                             >  2:2530382/1‑150 (MQ=255)
                                                                                       ccccTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGaca                                                        <  1:3812962/150‑1 (MQ=255)
                                                                                         ccTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACAcc                                                      >  2:2131162/1‑150 (MQ=255)
                                                                                           tGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAg                                                    <  1:328475/150‑1 (MQ=255)
                                                                                                 gCCCGGAGAAGATCGCCGCAGACGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACACCGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGa                                              <  2:10284764/150‑1 (MQ=255)
                                                                                                  cccGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGaa                                             <  2:973573/150‑1 (MQ=255)
                                                                                                     ggAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAaaag                                          <  2:10424853/150‑1 (MQ=255)
                                                                                                      gagaAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAaaaga                                         <  1:170872/150‑1 (MQ=255)
                                                                                                      gagaAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCGGGacacga                                         >  2:8012339/1‑145 (MQ=255)
                                                                                                       agaagaTCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAaaagaa                                        <  1:655838/150‑1 (MQ=255)
                                                                                                         aagaTCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAaaagaaag                                      <  2:9241376/150‑1 (MQ=255)
                                                                                                         aagaTCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAaaagaaag                                      <  1:6330763/150‑1 (MQ=255)
                                                                                                               gccgcAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCata                                <  2:2707839/150‑1 (MQ=255)
                                                                                                               gccgcAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCata                                >  1:1546665/1‑150 (MQ=255)
                                                                                                               gccgcAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCata                                <  1:5839488/150‑1 (MQ=255)
                                                                                                               gccgcAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCata                                >  2:5931242/1‑150 (MQ=255)
                                                                                                                ccgcAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCatat                               >  1:6013353/1‑150 (MQ=255)
                                                                                                                ccgcAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCatat                               <  1:9292450/150‑1 (MQ=255)
                                                                                                                   cAGAAGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATAtcgt                            <  1:3986004/150‑1 (MQ=255)
                                                                                                                    aGAAGCCCTGCTGGCCCCCATTCGCCAGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATAtcgt                            >  1:4993534/1‑150 (MQ=255)
                                                                                                                       aGCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCgag                        <  2:7188302/150‑1 (MQ=255)
                                                                                                                        gCCCTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCgaga                       >  1:4821449/1‑150 (MQ=255)
                                                                                                                          ccTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGAcc                     <  2:4211818/150‑1 (MQ=255)
                                                                                                                          ccTGCTGGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGAcc                     >  1:3764336/1‑150 (MQ=255)
                                                                                                                            tgctgGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGACCct                   <  2:6794468/150‑1 (MQ=255)
                                                                                                                             gctgGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGACCctg                  >  2:10132471/1‑150 (MQ=255)
                                                                                                                              ctgGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGACCctgc                 <  2:2941161/150‑1 (MQ=255)
                                                                                                                               tgGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGGCGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGACCctgct                <  2:1175866/150‑1 (MQ=255)
                                                                                                                               tgGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGACCctgct                <  2:7066348/150‑1 (MQ=255)
                                                                                                                               tgGCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGACCctgct                <  1:2449349/150‑1 (MQ=255)
                                                                                                                                ggCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGACCctgctg               <  2:4714321/150‑1 (MQ=255)
                                                                                                                                 gCCCCCATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGACCCTGCTGt              <  1:619516/150‑1 (MQ=255)
                                                                                                                                   ccccATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCAGATCGTCGAGACCCTGCTGTcc            >  2:608093/1‑150 (MQ=255)
                                                                                                                                      cATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGACCCTGCTGTCCATc         <  1:7888226/150‑1 (MQ=255)
                                                                                                                                      cATTCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGACCCTGCTGTCCATc         >  2:7729758/1‑150 (MQ=255)
                                                                                                                                        ttCGCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGACCCTGCTGTCCATCGc       <  1:7459913/150‑1 (MQ=255)
                                                                                                                                           gCCAGCCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGACCCTGCTGTCCATCGCCCa    <  1:10331405/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CAGACCCTGTACTACCGCCTGGTGCTGACCAAGCGCTACACCGGCGAGAAGACCAAGGAGAAAGGTCCGATCTTCCGCGACAGCCTGCCCCTGGAAGGCCCGGAGAAGATCGCCGCAGAAGCCCTGCTGGCCCCCATTCGCCAGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGACCCTGCTGTCCATCGCCCACG  >  NC_008463/3722783‑3723074

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