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% 1774.1 / ‑6.0 454K173E (AAG→GAG) PA14_RS22840DUF853 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (226/227);  minor base G (0/1);  total (226/228)
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 454 total aligned reads displayed.

CGCCCCCGCCCGTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTTGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGTTGGC  >  NC_008463/4997638‑4997933
                                                                                                                                                    |                                                                                                                                                   
   ccccGCCCGTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGa                                                                                                                                                >  1:10166010/1‑151 (MQ=255)
     ccgcccgTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGTGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAgg                                                                                                                                              <  2:1756947/151‑1 (MQ=255)
     ccgcccgACCGGCTCGGGGCGGAGGATGACTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCGGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAgg                                                                                                                                              <  2:2659301/151‑1 (MQ=255)
           gTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTcc                                                                                                                                        >  1:2251523/1‑151 (MQ=255)
           gTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTcc                                                                                                                                        <  2:3791509/151‑1 (MQ=255)
            tCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCg                                                                                                                                       <  2:7928704/151‑1 (MQ=255)
             ccGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCgg                                                                                                                                      <  1:6208363/151‑1 (MQ=255)
               ggCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGAt                                                                                                                                    >  1:5665306/1‑151 (MQ=255)
               ggCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGAt                                                                                                                                    <  1:3249970/151‑1 (MQ=255)
                gCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGAtt                                                                                                                                   >  1:6495191/1‑151 (MQ=255)
                 cTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTg                                                                                                                                  >  2:3362537/1‑151 (MQ=255)
                   cGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTc                                                                                                                                >  2:2160570/1‑151 (MQ=255)
                   cGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTc                                                                                                                                <  1:7573055/151‑1 (MQ=255)
                    ggggCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTcc                                                                                                                               >  1:907860/1‑151 (MQ=255)
                    ggggCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTcc                                                                                                                               <  2:9495171/151‑1 (MQ=255)
                     gggCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCt                                                                                                                              <  2:1059018/151‑1 (MQ=255)
                      ggCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCtt                                                                                                                             <  2:7749670/151‑1 (MQ=255)
                       gCGGAGGATGTCTTCCCGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCCGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGACCAGCGCCCGGTCCTCGCCGCGGATCTCCGGATTGTCCGTc                                                                                                                            <  1:4084807/151‑1 (MQ=255)
                       gCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTc                                                                                                                            <  1:6343649/151‑1 (MQ=255)
                         ggaggaTGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAg                                                                                                                          <  2:828179/151‑1 (MQ=255)
                          gaggaTGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAgg                                                                                                                         >  2:8504727/1‑151 (MQ=255)
                              aTGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAggtggt                                                                                                                     <  1:11000580/151‑1 (MQ=255)
                                   ttCCAGTTGCAGGGCAGGCTCGCCGAGGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGagc                                                                                                                >  1:10436222/1‑151 (MQ=255)
                                   ttCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGagc                                                                                                                <  2:9805740/151‑1 (MQ=255)
                                   ttCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGagc                                                                                                                >  1:4824298/1‑151 (MQ=255)
                                    tCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCCGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGagca                                                                                                               <  2:2084830/151‑1 (MQ=255)
                                    tCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGagca                                                                                                               >  1:8036772/1‑151 (MQ=255)
                                     ccAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGagcag                                                                                                              <  2:630588/151‑1 (MQ=255)
                                       aGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAgcg                                                                                                            >  2:1483489/1‑151 (MQ=255)
                                           gcagggcaggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGcctt                                                                                                        <  1:2488601/151‑1 (MQ=255)
                                           gcagggcaggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGcctt                                                                                                        <  1:7010950/151‑1 (MQ=255)
                                            cagggcaggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttc                                                                                                       >  1:3183327/1‑151 (MQ=255)
                                              gggcaggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcag                                                                                                     >  2:4083901/1‑151 (MQ=255)
                                              gggcaggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcag                                                                                                     <  2:5444289/151‑1 (MQ=255)
                                              gggcaggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccgtcag                                                                                                     >  2:3148129/1‑151 (MQ=255)
                                                    gCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACCGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGGTGAGCCGCGccttcaggtcctt                                                                                               >  2:7484890/1‑151 (MQ=255)
                                                      tCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtccttct                                                                                             >  1:1637534/1‑150 (MQ=255)
                                                      tCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtccttca                                                                                             <  2:3564587/151‑1 (MQ=255)
                                                      tCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtccttca                                                                                             >  1:8717735/1‑151 (MQ=255)
                                                         ccGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGGGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGGTGAGCAGCGccgtcaggtccttcaggt                                                                                          >  2:10871602/1‑151 (MQ=255)
                                                         ccGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGGGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtccttcaggt                                                                                          <  1:7674298/151‑1 (MQ=255)
                                                          cGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtccttcaggac                                                                                         >  2:6017872/1‑149 (MQ=255)
                                                           gaagaaGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCg                                                                                        <  1:5921184/151‑1 (MQ=255)
                                                            aagaagGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGa                                                                                       <  2:4990831/151‑1 (MQ=255)
                                                              gaagGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGagc                                                                                     >  2:6822787/1‑151 (MQ=255)
                                                               aagGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGagca                                                                                    >  1:9263094/1‑151 (MQ=255)
                                                                agGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGagcag                                                                                   <  2:6984005/151‑1 (MQ=255)
                                                                 ggCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGagcagc                                                                                  >  1:2417193/1‑151 (MQ=255)
                                                                 ggCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGagcagc                                                                                  <  2:3672846/151‑1 (MQ=255)
                                                                   ccTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGagcagcag                                                                                >  1:2543164/1‑151 (MQ=255)
                                                                          gcCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCt                                                                         >  1:7391079/1‑151 (MQ=255)
                                                                           cccTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCtt                                                                        >  1:439030/1‑151 (MQ=255)
                                                                            ccTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTc                                                                       >  1:7139109/1‑151 (MQ=255)
                                                                            ccTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTc                                                                       >  2:5346691/1‑151 (MQ=255)
                                                                                   ctCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCg                                                                <  1:1872888/151‑1 (MQ=255)
                                                                                    tCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCgg                                                               <  1:8447952/151‑1 (MQ=255)
                                                                                      cAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGcc                                                             >  2:8925510/1‑151 (MQ=255)
                                                                                      cAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGcc                                                             >  2:8754737/1‑151 (MQ=255)
                                                                                      cAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGcc                                                             >  1:2165298/1‑151 (MQ=255)
                                                                                       agagTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCa                                                            <  1:5378017/151‑1 (MQ=255)
                                                                                         agTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCAcc                                                          <  1:8057385/151‑1 (MQ=255)
                                                                                           tGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTg                                                        >  1:2869232/1‑151 (MQ=255)
                                                                                           tGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTGCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTg                                                        >  2:2834257/1‑151 (MQ=255)
                                                                                                  gcggcgCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAggcg                                                 >  2:2526679/1‑151 (MQ=255)
                                                                                                  gcggcgCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAggcg                                                 <  2:7631177/151‑1 (MQ=255)
                                                                                                  gcggcgCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAggcg                                                 <  2:4670371/151‑1 (MQ=255)
                                                                                                    ggcgCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCTGCAGGCCTTCCCGGTCGGCCTCCTGGAAggcggc                                               <  1:1705137/151‑1 (MQ=255)
                                                                                                    ggcgCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAggcggc                                               <  2:9302424/151‑1 (MQ=255)
                                                                                                        cagcagCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCCGATAGTCCTTCAGGTGGTAGAGCAGCGCCTTCAGGTACTTCAGGTCGAGCATCAGCAGGCCATCCCGGTCGGCCACCGGGAAGGCGGCGAAc                                           >  2:3299878/1‑151 (MQ=255)
                                                                                                         agcagcGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCTGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACa                                          >  1:10794884/1‑151 (MQ=255)
                                                                                                         agcagcGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACa                                          >  1:2719743/1‑151 (MQ=255)
                                                                                                         agcagcGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACa                                          <  1:10276052/151‑1 (MQ=255)
                                                                                                          gcagcGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAg                                         >  1:3373294/1‑151 (MQ=255)
                                                                                                           cagcGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCGTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAgc                                        <  2:4838657/151‑1 (MQ=255)
                                                                                                            agcGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAgcg                                       >  2:5154339/1‑151 (MQ=255)
                                                                                                             gcgcCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAgcgc                                      >  1:8280598/1‑151 (MQ=255)
                                                                                                             gcgcCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAgcgc                                      <  1:514742/151‑1 (MQ=255)
                                                                                                              cgcCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGcgcg                                     >  1:10826305/1‑151 (MQ=255)
                                                                                                              cgcCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGcgcg                                     <  2:1690649/151‑1 (MQ=255)
                                                                                                                ccTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGc                                   <  2:6674521/151‑1 (MQ=255)
                                                                                                                         gaGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGGTGAGCAGCGCCTTCAGGTCCTTCAGGTCGTGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGcg                          >  1:4510170/2‑150 (MQ=255)
                                                                                                                          aGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTa                         <  1:6820998/151‑1 (MQ=255)
                                                                                                                            cgcCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATc                       <  1:10008536/151‑1 (MQ=255)
                                                                                                                               cGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGt                    <  2:2429278/151‑1 (MQ=255)
                                                                                                                               cGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGt                    >  1:9948290/1‑151 (MQ=255)
                                                                                                                                ggTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTc                   <  2:1273844/151‑1 (MQ=255)
                                                                                                                                    aaCAGCGCCCGGTCCTCGCCGCGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACt               <  2:2929709/151‑1 (MQ=255)
                                                                                                                                    aaCAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTACAGGTGGTTGAGCCGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACt               <  1:7283757/151‑1 (MQ=255)
                                                                                                                                       aGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCa            >  1:4995191/1‑151 (MQ=255)
                                                                                                                                        gcgcCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCag           <  2:427601/151‑1 (MQ=255)
                                                                                                                                        gcgcCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGGCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCag           <  2:5946657/151‑1 (MQ=255)
                                                                                                                                          gcCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCagga         <  1:2069743/151‑1 (MQ=255)
                                                                                                                                            ccGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCaggagg       <  2:5614781/151‑1 (MQ=255)
                                                                                                                                            ccGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCaggagg       >  2:10160467/1‑151 (MQ=255)
                                                                                                                                              ggTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGtt     >  1:385733/1‑151 (MQ=255)
                                                                                                                                              ggTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGtt     <  1:1440166/151‑1 (MQ=255)
                                                                                                                                              ggTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGtt     >  2:7179956/1‑151 (MQ=255)
                                                                                                                                               gTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGTTg    <  1:1135951/151‑1 (MQ=255)
                                                                                                                                                tCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGTTgg   <  2:10109427/151‑1 (MQ=255)
                                                                                                                                                tCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGTTgg   <  2:10773713/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGCCCCCGCCCGTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTTGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGTTGGC  >  NC_008463/4997638‑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: ‑