Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 4,997,786 T→C K173E (AAG→GAG)  PA14_RS22840 ← DUF853 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084634,997,7860TC99.8% 1997.5 / NA 507K173E (AAG→GAG) PA14_RS22840DUF853 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (222/284);  minor base A (0/1);  total (222/285)
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 507 total aligned reads displayed.

GCCCCCGCCCGTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTTGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGTTGGC  >  NC_008463/4997639‑4997933
                                                                                                                                                   |                                                                                                                                                   
gcCCCCGCCCGTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGcc                                                                                                                                                  <  2:3391448/151‑1 (MQ=255)
   cccgcccgTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAg                                                                                                                                               <  2:899662/151‑1 (MQ=255)
   cccgcccgTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAg                                                                                                                                               <  2:6693511/151‑1 (MQ=255)
    ccgcccgTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAgg                                                                                                                                              <  2:10430806/151‑1 (MQ=255)
     cgcccgTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGa                                                                                                                                             <  2:2381350/151‑1 (MQ=255)
     cgcccgTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGa                                                                                                                                             >  1:8050387/1‑151 (MQ=255)
      gcccgTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGAt                                                                                                                                            >  1:4896467/1‑151 (MQ=255)
      gcccgTCCGGCTCGGGGCGGAGGATGACTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGACCAGCGCCCGGTCCTCGCCGAgtca                                                                                                                                            >  2:3050732/1‑148 (MQ=255)
           tCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCg                                                                                                                                       <  2:5770158/151‑1 (MQ=255)
            ccGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCgg                                                                                                                                      >  1:4312771/1‑151 (MQ=255)
              ggCGCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGAt                                                                                                                                    >  2:4547089/1‑151 (MQ=255)
               gCTCGGGGCGGAGGATGTCTTCCTGTTGCAGGGCAGGCTCGCCGAAGATGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTAGCCGAGGATCTCCGGAtt                                                                                                                                   <  2:8047650/151‑1 (MQ=255)
               gCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGAtt                                                                                                                                   >  1:9525945/1‑151 (MQ=255)
                cTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTg                                                                                                                                  <  2:10043817/151‑1 (MQ=255)
                cTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTg                                                                                                                                  <  1:2240414/151‑1 (MQ=255)
                 tCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCACGGCGCCCTGCTGCTCCAGCGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGACCAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGt                                                                                                                                 <  2:4064010/151‑1 (MQ=255)
                   ggggCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTcc                                                                                                                               <  2:9100629/151‑1 (MQ=255)
                     ggCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCtt                                                                                                                             <  1:6042398/151‑1 (MQ=255)
                      gCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTc                                                                                                                            <  2:6122766/151‑1 (MQ=255)
                       cGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCa                                                                                                                           >  1:2940115/1‑151 (MQ=255)
                       cGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCa                                                                                                                           >  1:11151674/1‑151 (MQ=255)
                       cGGAGGATGTCGTCCAGTTGCAGGGCAGGCTCGCCGAAGACGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGACCAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCa                                                                                                                           <  2:1609840/151‑1 (MQ=255)
                        ggaggaGGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAg                                                                                                                          >  1:8493318/1‑151 (MQ=255)
                          aggaTGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAggt                                                                                                                        >  1:2429065/1‑151 (MQ=255)
                            gaTGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAggtgg                                                                                                                      >  1:6820517/1‑151 (MQ=255)
                             aTGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAggtggt                                                                                                                     >  1:6964513/1‑151 (MQ=255)
                             aTGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAggtggt                                                                                                                     <  2:249247/151‑1 (MQ=255)
                              tGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGtt                                                                                                                    <  2:6468964/151‑1 (MQ=255)
                                 cTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGag                                                                                                                 >  2:4850877/1‑151 (MQ=255)
                                     cAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGagcagc                                                                                                             <  1:9428918/151‑1 (MQ=255)
                                      aGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAgcg                                                                                                            <  2:6859586/151‑1 (MQ=255)
                                         tGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGcct                                                                                                         <  1:8624857/151‑1 (MQ=255)
                                          gcagggcaggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGGACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGcctt                                                                                                        >  2:10104393/1‑151 (MQ=255)
                                          gcagggcaggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGcctt                                                                                                        <  1:6524360/151‑1 (MQ=255)
                                              ggcaggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcagg                                                                                                    <  2:5668074/151‑1 (MQ=255)
                                               gcaggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggt                                                                                                   <  1:10004002/151‑1 (MQ=255)
                                                 aggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGGAGCGCCGGTGCACAGCGCCCGCTCCTCGACGAGGATCTCCGGATTGACCTACTGGTCGTGGAGCATCGccatcaggtcg                                                                                                 >  1:139092/1‑150 (MQ=255)
                                                  ggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtcct                                                                                                <  2:4112207/151‑1 (MQ=255)
                                                  ggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtcct                                                                                                >  2:8091561/1‑151 (MQ=255)
                                                   gCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtcctt                                                                                               <  2:1712987/151‑1 (MQ=255)
                                                   gCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtcctt                                                                                               <  2:7578085/151‑1 (MQ=255)
                                                    cTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtcctcc                                                                                              <  1:1789086/151‑3 (MQ=255)
                                                      cGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtccttcag                                                                                            <  2:2108788/151‑1 (MQ=255)
                                                      cGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtccttcag                                                                                            <  1:3826991/151‑1 (MQ=255)
                                                       gCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtccttcagg                                                                                           >  2:10270214/1‑151 (MQ=255)
                                                        ccGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtccttcaggt                                                                                          >  2:7740907/1‑151 (MQ=255)
                                                         cGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGGGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCcttcaggtccttcaggtc                                                                                         <  2:3793703/151‑1 (MQ=255)
                                                          gaagaaGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCg                                                                                        >  1:1011786/1‑151 (MQ=255)
                                                          gaagaaGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCg                                                                                        >  1:5402662/1‑151 (MQ=255)
                                                            agaagGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGag                                                                                      <  1:4068663/151‑1 (MQ=255)
                                                            agaagGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGag                                                                                      >  1:5810113/1‑151 (MQ=255)
                                                             gaagGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTGGAGCAGCGCCTTCCGGTCCTTCAGGTCGagc                                                                                     >  2:8893932/1‑151 (MQ=255)
                                                              aagGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGACCAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGagca                                                                                    >  2:4985357/1‑151 (MQ=255)
                                                              aagGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGagca                                                                                    <  1:5214239/151‑1 (MQ=255)
                                                                 gCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGagcagca                                                                                 <  2:9706867/151‑1 (MQ=255)
                                                                 gCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGagcagca                                                                                 >  2:6382339/1‑151 (MQ=255)
                                                                  ccTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGagcagcag                                                                                <  1:2702613/151‑1 (MQ=255)
                                                                   cTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGagcagcagc                                                                               >  1:9049851/1‑151 (MQ=255)
                                                                     cGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGcagcagcag                                                                             <  1:1592857/151‑1 (MQ=255)
                                                                      ggCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAgg                                                                            >  1:5945158/1‑151 (MQ=255)
                                                                       gcgcCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGc                                                                           <  2:7241480/151‑1 (MQ=255)
                                                                         gcCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCt                                                                         >  2:1147977/1‑151 (MQ=255)
                                                                         gcCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCt                                                                         >  1:10706885/1‑151 (MQ=255)
                                                                           ccTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTc                                                                       <  1:9317509/151‑1 (MQ=255)
                                                                                tgctCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGt                                                                  >  2:2665917/1‑151 (MQ=255)
                                                                                  ctCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCg                                                                <  2:10586575/151‑1 (MQ=255)
                                                                                   tCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCgg                                                               >  1:3460588/1‑151 (MQ=255)
                                                                                            gCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGa                                                      >  1:8856417/1‑151 (MQ=255)
                                                                                            gCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGa                                                      >  2:10963455/1‑151 (MQ=255)
                                                                                             ccAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGaa                                                     >  2:4718650/1‑151 (MQ=255)
                                                                                              cAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAg                                                    <  1:9578604/151‑1 (MQ=255)
                                                                                              cAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAg                                                    <  2:4353507/151‑1 (MQ=255)
                                                                                              cAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAg                                                    >  2:7389535/1‑151 (MQ=255)
                                                                                               aGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAgg                                                   <  1:10211833/151‑1 (MQ=255)
                                                                                                 gcggcgCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTTTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGACGGCCACCTGGAAggcg                                                 <  2:11212867/151‑1 (MQ=255)
                                                                                                  cggcgCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAggcgg                                                <  1:4377689/151‑1 (MQ=255)
                                                                                                  cggcgCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAggcgg                                                >  1:10360702/1‑151 (MQ=255)
                                                                                                  cggcgCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAggcgg                                                <  2:5186428/151‑1 (MQ=255)
                                                                                                   ggcgCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAggcggc                                               >  1:2547459/1‑151 (MQ=255)
                                                                                                        agcagcGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACa                                          >  2:5883423/1‑151 (MQ=255)
                                                                                                          cagcGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAgc                                        <  1:7814975/151‑1 (MQ=255)
                                                                                                           agcGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAgcg                                       <  2:8391060/151‑1 (MQ=255)
                                                                                                              gcCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCgg                                    >  2:3962444/1‑151 (MQ=255)
                                                                                                                 tGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCct                                 <  2:2854954/151‑1 (MQ=255)
                                                                                                                  gCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCctg                                >  2:2637401/1‑151 (MQ=255)
                                                                                                                     cGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCctgctg                             <  1:3283926/151‑1 (MQ=255)
                                                                                                                      ggAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTgg                            <  1:8330638/151‑1 (MQ=255)
                                                                                                                        aaGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCt                          >  1:9766169/1‑151 (MQ=255)
                                                                                                                         aGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGACGGCGGCGAACAGCGCGGCCTGCTGGCTa                         >  1:2057478/1‑151 (MQ=255)
                                                                                                                               ggTGAACAGCGCCCGGTCCTCGCCGTGGATCTCCGGAGTGTCCTTCAGGTGGTTGAGCAGCGCCTTCTGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGACCAGCGCGGCCTGCTGGCTATCGGTc                   >  1:10882843/1‑151 (MQ=255)
                                                                                                                               ggTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTc                   >  2:3211449/1‑151 (MQ=255)
                                                                                                                                 tGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCaa                 <  1:3741623/151‑1 (MQ=255)
                                                                                                                                   aaCAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTCTCGGTCAAcg               >  2:10861552/1‑150 (MQ=255)
                                                                                                                                      cgcgcCCGGTCCTCGCCGCCGATCGCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCa            <  2:2593633/150‑1 (MQ=255)
                                                                                                                                         gcCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCagga         >  2:5498757/1‑151 (MQ=255)
                                                                                                                                         gcCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCagga         <  1:1211622/151‑1 (MQ=255)
                                                                                                                                          cccGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCaggag        >  2:8553601/1‑151 (MQ=255)
                                                                                                                                          cccGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCaggag        <  2:9639173/151‑1 (MQ=255)
                                                                                                                                           ccGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCaggagg       <  1:5511172/151‑1 (MQ=255)
                                                                                                                                                ccTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCTGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGTTGGc  >  1:4539764/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCCCCCGCCCGTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTTGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGTTGGC  >  NC_008463/4997639‑4997933

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: ‑