Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 2,149,425 T→C Q549Q (CAA→CAG PA14_RS09965 ← M14 family metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084632,149,4250TC99.8% 1612.1 / NA 420Q549Q (CAA→CAGPA14_RS09965M14 family metallopeptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (197/222);  minor base A (0/1);  total (197/223)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 420 total aligned reads displayed.

GCGCAGGCGTTTGTCGAAGGCGGACTTCAGCGCCGGGTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTTTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCGTGCCGCTGACCAGCAGGCGACCTTCGTAGCCCGGAATCCGGTT  >  NC_008463/2149278‑2149571
                                                                                                                                                   |                                                                                                                                                  
gcgcAGGCGTTTGTCGAAGGCGGACTTCAGCGCCGGGTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTg                                                                                                                                                  <  1:1876277/150‑1 (MQ=255)
gcgcAGGCGTTTGTCGAAGGCGGACTTCAGCGCCGGGTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTg                                                                                                                                                  <  2:4786725/150‑1 (MQ=255)
gcgcAGGCGTTTGTCGAAGGCGGACTTCAGCGCCGGGTCTTCCGCCAGCATCGTGAGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTg                                                                                                                                                  <  2:227529/150‑1 (MQ=255)
   cAGGCGTTTGTCGAAGGCGGACTACAGCGCCGGGTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCgg                                                                                                                                               <  2:4925690/150‑1 (MQ=255)
    aGGCGTTTGTCGAAGGCGGACTTCAGCGCCGGGTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGt                                                                                                                                              >  1:5134317/1‑150 (MQ=255)
    aGGCGTTTGTCGAAGGCGGACTTCAGCGCCGGGTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGt                                                                                                                                              >  2:8808399/1‑150 (MQ=255)
    aGGCGTTTGTCGAAGGCGGACTTCAGCGCCGGGACTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTAGTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGt                                                                                                                                              <  2:7054320/150‑1 (MQ=255)
     ggCGTTTGTCGAAGGCGGACTTCAGCGCCGGGTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGtt                                                                                                                                             >  1:2636611/1‑150 (MQ=255)
        gTTTGTCGAAGGCGGACTTCAGCGCCGGGTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTgca                                                                                                                                          >  2:9475345/1‑150 (MQ=255)
             tCGAAGGCGGACTTCAGCGCCGGGTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCCGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGTATCCGGACTCTGCGGTTGCAGCAgg                                                                                                                                     <  2:2651130/150‑1 (MQ=255)
              cgcaGGCGGACTTCAGCGCCGGGTCTTCCGCCAGCATCTGGCGCGCCATCGGCTCCCTCACGAATTCCTCCGGCTCCGCGGCGGCCACCCGGGTGGAATTGAAGAAGCCCCAGGACCAGACCGAATCCGGACTCTGCGGTTGCAGCAGGt                                                                                                                                    <  2:8159402/147‑1 (MQ=255)
                 aGGCGGACTTCAGCGCCGGGTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGa                                                                                                                                 <  2:9769478/150‑1 (MQ=255)
                 aGGCGGACTTCAGCGCCGGGTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGa                                                                                                                                 <  1:9627495/150‑1 (MQ=255)
                 aGGCGGACTTCAGCGCCGGGTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCCGACTCTGCGGTTGCAGCAGGTCGa                                                                                                                                 >  2:3259045/1‑150 (MQ=255)
                   gCGGACTTCAGCGCCGGGTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATg                                                                                                                               <  1:4080366/150‑1 (MQ=255)
                         ttCAGCGCCGGGTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAgg                                                                                                                         >  1:4382662/1‑150 (MQ=255)
                          tCAGCGCCGGGTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGa                                                                                                                        <  2:9351616/150‑1 (MQ=255)
                             gcgcCGGGTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGc                                                                                                                     <  2:421413/150‑1 (MQ=255)
                              cgcCGGGTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGcc                                                                                                                    <  1:7465936/150‑1 (MQ=255)
                               gcCGGGTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGccc                                                                                                                   <  1:3785327/150‑1 (MQ=255)
                               gcCGGGTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTAGTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGccc                                                                                                                   <  1:10530783/150‑1 (MQ=255)
                                ccGGGTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCa                                                                                                                  <  2:1989291/150‑1 (MQ=255)
                                  gggTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGc                                                                                                                <  2:10665848/150‑1 (MQ=255)
                                   ggTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGACTTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCg                                                                                                               >  1:8745533/1‑150 (MQ=255)
                                   ggTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCg                                                                                                               <  1:1549423/150‑1 (MQ=255)
                                   ggTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCg                                                                                                               >  2:2514305/1‑150 (MQ=255)
                                   ggTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCg                                                                                                               >  1:7913112/1‑150 (MQ=255)
                                      cTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCt                                                                                                            >  2:6347064/1‑150 (MQ=255)
                                       ttCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTg                                                                                                           >  2:10420427/1‑150 (MQ=255)
                                        tCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTgg                                                                                                          >  1:2779967/1‑150 (MQ=255)
                                         ccGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTggg                                                                                                         <  1:9320066/150‑1 (MQ=255)
                                            ccAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCgg                                                                                                      <  1:6800747/150‑1 (MQ=255)
                                               gCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGgtgt                                                                                                   >  1:3235029/1‑150 (MQ=255)
                                                 aTCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGtcg                                                                                                 >  1:9073408/1‑150 (MQ=255)
                                                     tGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGtcgatcg                                                                                             >  1:7263226/1‑150 (MQ=255)
                                                       cgcgcCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGAc                                                                                           <  1:8549518/150‑1 (MQ=255)
                                                        gcgcCATCGGCTCCATCACGTCTTCCTCCGGCACCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACGCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACa                                                                                          <  2:8299099/150‑1 (MQ=255)
                                                          gcCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGc                                                                                        <  1:6948642/150‑1 (MQ=255)
                                                          gcCATCGGCTCCATCACGTAGTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGc                                                                                        <  1:2171488/150‑1 (MQ=255)
                                                           ccATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCg                                                                                       <  2:3713373/150‑1 (MQ=255)
                                                            cATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCTCCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGa                                                                                      <  2:5548432/150‑1 (MQ=255)
                                                            cATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGa                                                                                      >  2:567616/1‑150 (MQ=255)
                                                             aTCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAg                                                                                     >  1:3660411/1‑150 (MQ=255)
                                                             aTCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAg                                                                                     >  2:214858/1‑150 (MQ=255)
                                                             aTCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAg                                                                                     >  1:9735261/1‑150 (MQ=255)
                                                              tCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAgc                                                                                    <  2:5156150/150‑1 (MQ=255)
                                                               cGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAgcc                                                                                   >  1:6189004/1‑150 (MQ=255)
                                                                 gCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAgccgg                                                                                 >  2:588005/1‑150 (MQ=255)
                                                                  cTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAgccggc                                                                                >  2:6864205/1‑150 (MQ=255)
                                                                  cTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAgccggc                                                                                <  1:4740240/150‑1 (MQ=255)
                                                                   tCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAgccggcc                                                                               >  1:5628183/1‑150 (MQ=255)
                                                                    ccATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAgccggccg                                                                              <  1:3953114/150‑1 (MQ=255)
                                                                     cATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGccggccgg                                                                             <  2:9039572/150‑1 (MQ=255)
                                                                      aTCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCggg                                                                            >  2:8515679/1‑150 (MQ=255)
                                                                       tCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGa                                                                           >  2:1280944/1‑150 (MQ=255)
                                                                            tATTCCTCCGGCTCCTCGGCGGCCACCCGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTc                                                                      <  2:1409794/150‑1 (MQ=255)
                                                                            tATTCCTCCGGCTCCTCGGCGGCCACCATGGTGGAATTGTAGAAGCCCCAGGACCTGATCGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCCGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGta                                                                      <  2:6484332/150‑2 (MQ=255)
                                                                             agtcctccGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGACTCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCt                                                                     <  2:7947940/148‑1 (MQ=255)
                                                                             aTTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCt                                                                     >  1:4889772/1‑150 (MQ=255)
                                                                               tcctccGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGc                                                                   >  1:2062103/1‑150 (MQ=255)
                                                                                   ccGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAgc                                                               <  2:7597555/150‑1 (MQ=255)
                                                                                    cGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAgcg                                                              >  1:3848824/1‑150 (MQ=255)
                                                                                      gCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCt                                                            <  1:5768942/150‑1 (MQ=255)
                                                                                      gCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCt                                                            <  2:10955819/150‑1 (MQ=255)
                                                                                       ctcctcGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTg                                                           <  1:809564/150‑1 (MQ=255)
                                                                                       ctcctcGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTg                                                           <  2:3860820/150‑1 (MQ=255)
                                                                                        tcctcGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTgg                                                          >  1:2957613/1‑150 (MQ=255)
                                                                                           tcGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAAc                                                       >  2:5098650/1‑150 (MQ=255)
                                                                                            cggcggCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACg                                                      >  1:7597518/1‑150 (MQ=255)
                                                                                            cggcggCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACg                                                      <  2:6539267/150‑1 (MQ=255)
                                                                                             ggcggcCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACgg                                                     <  1:10939754/150‑1 (MQ=255)
                                                                                              gcggcCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGc                                                    >  2:7867685/1‑150 (MQ=255)
                                                                                                  ccaccaGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCg                                                <  2:1718583/150‑1 (MQ=255)
                                                                                                    accaGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCCGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGc                                              >  2:2843345/1‑150 (MQ=255)
                                                                                                     ccaGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCg                                             >  1:4239819/1‑150 (MQ=255)
                                                                                                     ccaGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCg                                             >  1:3477606/1‑150 (MQ=255)
                                                                                                      caGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCGt                                            >  2:6244138/1‑150 (MQ=255)
                                                                                                       aGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGACGGTCTGCAAGCGCTGGAACGGCTTCGGCGTg                                           >  2:8052520/1‑150 (MQ=255)
                                                                                                         ggTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCGTGcc                                         <  1:6341999/150‑1 (MQ=255)
                                                                                                          gTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCGTGCCg                                        >  1:7211448/1‑150 (MQ=255)
                                                                                                               aTTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCGTGCCGCTGAc                                   <  1:10006558/150‑1 (MQ=255)
                                                                                                                ttGAAGAAGCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCGTGCCGCTGAcc                                  <  1:10815490/150‑1 (MQ=255)
                                                                                                                   aagaagCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCGTGCCGCTGACcagc                               <  2:8429474/150‑1 (MQ=255)
                                                                                                                   aagaagCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCGTGCCGCTGACcagc                               <  2:10737881/150‑1 (MQ=255)
                                                                                                                   aagaagCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCGTGCCGCTGACcagc                               <  2:7255418/150‑1 (MQ=255)
                                                                                                                     gaagCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCGTGCCGCTGACcagcag                             <  2:4382662/150‑1 (MQ=255)
                                                                                                                      aagCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCGTGCCGCTGACCAGCAgg                            >  1:11096751/1‑150 (MQ=255)
                                                                                                                       agCCCCAGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCGTGCCGCTGACCAGCAGGc                           >  1:6285405/1‑150 (MQ=255)
                                                                                                                             aGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCGTGCCGCTGACCAGCAGGCGACCtt                     <  1:6011665/150‑1 (MQ=255)
                                                                                                                             aGGACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCGTGCCGCTGACCAGCAGGCGACCtt                     <  2:4591834/150‑1 (MQ=255)
                                                                                                                              ggACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCGTGCCGCTGACCAGCAGGCGACCTTc                    <  1:1604208/150‑1 (MQ=255)
                                                                                                                              ggACCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCGTGCCGCTGACCAGCAGGCGACCTTc                    <  1:8394394/150‑1 (MQ=255)
                                                                                                                                aCCAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCGTGCCGCTGACCAGCAGGCGACCTTCGt                  >  1:9892646/1‑150 (MQ=255)
                                                                                                                                 ccAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCGTGCCGCTGACCAGCAGGCGACCTTCGTa                 >  1:7150612/1‑150 (MQ=255)
                                                                                                                                 ccAGAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCGTGCCGCTGACCAGCAGGCGACCTTCGTa                 >  1:287658/1‑150 (MQ=255)
                                                                                                                                    gAACGAATCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCGTGCCGCTGACCAGCAGGCGACCTTCGTAGcc              <  1:5387204/150‑1 (MQ=255)
                                                                                                                                         aaTCCGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCGTGCCGCTGACCAGCAGGCGACCTTCGTAGCCCGGaa         >  1:886982/1‑150 (MQ=255)
                                                                                                                                             cGGACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGTGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCGTGCCGCTGACCAGCAGGCGACCTTCGTAGCCCGGAATCCg     >  2:6034480/1‑150 (MQ=255)
                                                                                                                                              ggACTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCGTCGGCGTGCCGCTGACCAGCAGGCGACCTTCGTAGCCCGGAATCCgg    >  2:6621136/1‑150 (MQ=255)
                                                                                                                                                aCTCTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCGTGCCGCTGACCAGCAGGCGACCTTCGTAGCCCGGAATCCGGtt  >  2:1145947/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCGCAGGCGTTTGTCGAAGGCGGACTTCAGCGCCGGGTCTTCCGCCAGCATCTTGCGCGCCATCGGCTCCATCACGTATTCCTCCGGCTCCTCGGCGGCCACCAGGGTGGAATTGAAGAAGCCCCAGGACCAGAACGAATCCGGACTTTGCGGTTGCAGCAGGTCGATGGCCAGGACGCCCAGCGGCTGGGCGGTGTCGATCGACAGCGAGCCGGCCGGGAAGGTCTGCAAGCGCTGGAACGGCTTCGGCGTGCCGCTGACCAGCAGGCGACCTTCGTAGCCCGGAATCCGGTT  >  NC_008463/2149278‑2149571

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: