Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 559,051 +T intergenic (+619/+70) PA14_RS02545 → / ← rarD glycine cleavage system protein R/EamA family transporter RarD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_008463559,0501.T98.8% 1427.3 / 5.7 342intergenic (+618/+71)PA14_RS02545/rarDglycine cleavage system protein R/EamA family transporter RarD
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base T (172/165);  total (176/166)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 372 total aligned reads displayed.

ACAGCGGTGTGCTCAAAGACTGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACGGCCAGGGCGAGCCAGATGAA  >  NC_008463/558900‑559195
                                                                                                                                                        |                                                                                                                                                 
   gCGGTGTGCTCAAAGACTGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTAttg                                                                                                                                                >  2:9643622/1‑149 (MQ=255)
    cGGTGTGCTCAAAGACTGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTAttgg                                                                                                                                               <  2:6165779/151‑4 (MQ=255)
       tgtgCTCAAAGACTGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCgg                                                                                                                                             <  1:7767059/151‑1 (MQ=255)
        gtgCTCAAAGACTGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCggg                                                                                                                                            <  1:6381159/151‑1 (MQ=255)
        gtgCTCAAAGACTGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCggg                                                                                                                                            <  2:5089126/151‑1 (MQ=255)
          gCTCAAAGACTGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCggggg                                                                                                                                          >  1:1166330/1‑151 (MQ=255)
             cAAAGACTGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATa                                                                                                                                       >  1:3678474/1‑151 (MQ=255)
             cAAAGACTGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATa                                                                                                                                       <  2:1482360/151‑1 (MQ=255)
              aaaGACTGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATag                                                                                                                                      <  2:4131607/151‑1 (MQ=255)
                aGACTGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATagga                                                                                                                                    >  1:2175810/1‑151 (MQ=255)
                     gATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGAAGGTATATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATAGCTGAGCGGCAGTATTGGCGGGGGAGAGGAAGGAg                                                                                                                               >  2:2223211/1‑151 (MQ=255)
                       tCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTg                                                                                                                             <  2:10634973/151‑1 (MQ=255)
                       tCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTg                                                                                                                             <  1:7557191/151‑1 (MQ=255)
                        cAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGc                                                                                                                            >  2:11156641/1‑151 (MQ=255)
                           gcgcTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTCTTGGCGGGGGATAGGAAGGAGTGCgtg                                                                                                                         >  2:6560789/1‑151 (MQ=255)
                              cTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGttt                                                                                                                      >  1:8035322/1‑151 (MQ=255)
                                   cAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCg                                                                                                                 >  2:183097/1‑151 (MQ=255)
                                     gggCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGcc                                                                                                               >  1:290317/1‑151 (MQ=255)
                                     gggCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGcc                                                                                                               >  2:4629849/1‑151 (MQ=255)
                                      ggCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCTGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCa                                                                                                              >  1:1325671/1‑151 (MQ=255)
                                      ggCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGGCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCa                                                                                                              >  2:8124672/1‑151 (MQ=255)
                                        cttcttCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAgg                                                                                                            <  2:306171/151‑1 (MQ=255)
                                         ttcttcCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGc                                                                                                           >  1:4092860/1‑151 (MQ=255)
                                          tcttcCGACGAGTTATCCACGCGCCTGCGGGGATGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCt                                                                                                          <  1:751175/151‑1 (MQ=255)
                                          tcttcCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCt                                                                                                          >  1:5407485/1‑151 (MQ=255)
                                             tcCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGc                                                                                                       <  1:7631925/151‑1 (MQ=255)
                                               cgacgaGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCgg                                                                                                     <  1:4952807/151‑1 (MQ=255)
                                                gacgaGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCGTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGa                                                                                                    >  1:8462332/1‑151 (MQ=255)
                                                   gagaaaTCCACGCTACTGCTGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGTGGGATAGGAAGGACTGCGTGTTTTTTCGCCAGGCTTGCGGAGtt                                                                                                 >  2:1361694/6‑151 (MQ=255)
                                                    aGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTc                                                                                                >  2:264236/1‑151 (MQ=255)
                                                     gTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGGATTGGCGGGGGATAGGAAGGAGTGCGTGTTGTTTCGCCAGGCTTGCGGAGTTcc                                                                                               >  2:8225220/1‑151 (MQ=255)
                                                      ttATCCACGCGCCTGCGGGGAGTAGAAGGCAGGTCTAAGCACAAATTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGATTTTTCGCCAGGCTTGCGGAGTTCCg                                                                                              >  1:6444694/1‑151 (MQ=255)
                                                      ttATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCg                                                                                              <  1:9856974/151‑1 (MQ=255)
                                                      ttATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTAGTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCg                                                                                              >  2:9154833/1‑151 (MQ=255)
                                                        aTCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCa                                                                                            <  2:3571301/151‑1 (MQ=255)
                                                           cACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGc                                                                                         <  1:7124351/151‑1 (MQ=255)
                                                            aCGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGcc                                                                                        <  1:1500581/151‑1 (MQ=255)
                                                            aCGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGcc                                                                                        >  2:7164783/1‑151 (MQ=255)
                                                              gcgcCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTg                                                                                      >  1:863920/1‑151 (MQ=255)
                                                               cgcCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTgg                                                                                     >  2:4876207/1‑151 (MQ=255)
                                                                 ccTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGTATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCATgggg                                                                                   >  1:4229947/1‑151 (MQ=255)
                                                                 ccTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTgggg                                                                                   <  2:2002238/151‑1 (MQ=255)
                                                                  cTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGa                                                                                  >  2:10239620/1‑151 (MQ=255)
                                                                  cTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGa                                                                                  >  1:3110197/1‑151 (MQ=255)
                                                                     cGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATaa                                                                               >  1:2442344/1‑151 (MQ=255)
                                                                       gggAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACg                                                                             <  1:6869290/151‑1 (MQ=255)
                                                                        ggaggaGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGt                                                                            >  1:10790758/1‑151 (MQ=255)
                                                                         gaggagAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTc                                                                           <  2:3986591/151‑1 (MQ=255)
                                                                           ggagAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAg                                                                         <  1:9601293/151‑1 (MQ=255)
                                                                            gagaAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGGATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGc                                                                        >  1:6695368/1‑151 (MQ=255)
                                                                             agaAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGcc                                                                       <  1:3357442/151‑1 (MQ=255)
                                                                             agaAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGcc                                                                       >  2:4624806/1‑151 (MQ=255)
                                                                              gaAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGccc                                                                      <  2:6992301/151‑1 (MQ=255)
                                                                               aaGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCg                                                                     >  1:3939754/1‑151 (MQ=255)
                                                                                 ggCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGc                                                                   <  1:1460358/151‑1 (MQ=255)
                                                                                  gCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCt                                                                  <  1:2971904/151‑1 (MQ=255)
                                                                                   cAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTg                                                                 <  2:4285203/151‑1 (MQ=255)
                                                                                   cAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTg                                                                 >  1:19221/1‑151 (MQ=255)
                                                                                    aGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTgg                                                                >  2:1754320/1‑151 (MQ=255)
                                                                                     ggTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTggg                                                               >  1:8913004/1‑151 (MQ=255)
                                                                                      gTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTgggg                                                              >  2:10132227/1‑151 (MQ=255)
                                                                                      gTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGAGAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTgggg                                                              >  2:882486/1‑151 (MQ=255)
                                                                                       tCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGa                                                             >  2:5839010/1‑151 (MQ=255)
                                                                                             cacaAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGAGTGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAcggcgg                                                       >  2:5449358/1‑151 (MQ=255)
                                                                                              acaAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAggcggc                                                      <  2:6406421/151‑1 (MQ=255)
                                                                                              acaAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAggcggc                                                      >  1:10152687/1‑151 (MQ=255)
                                                                                                aaaGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGgcgc                                                    >  2:830260/1‑151 (MQ=255)
                                                                                                 aaGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCa                                                   >  1:7971954/1‑151 (MQ=255)
                                                                                                  aGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAg                                                  <  1:10962822/151‑1 (MQ=255)
                                                                                                     aCGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCt                                               <  2:444912/151‑1 (MQ=255)
                                                                                                     aCGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCt                                               <  1:7220853/151‑1 (MQ=255)
                                                                                                      cGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTg                                              >  1:5917646/1‑151 (MQ=255)
                                                                                                           aTGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCa                                         >  2:10843703/1‑151 (MQ=255)
                                                                                                           aTGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCa                                         >  1:8861238/1‑151 (MQ=255)
                                                                                                            tGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAg                                        <  1:2011902/151‑1 (MQ=255)
                                                                                                             gAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCGTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAgg                                       <  2:7983674/151‑1 (MQ=255)
                                                                                                              aaCCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGc                                      >  2:947458/1‑151 (MQ=255)
                                                                                                               aCCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCg                                     <  1:9524948/151‑1 (MQ=255)
                                                                                                                  tgtgGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCg                                  >  2:365349/1‑151 (MQ=255)
                                                                                                                  tgtgGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCg                                  <  2:6779059/151‑1 (MQ=255)
                                                                                                                   gtgGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGa                                 >  2:1155832/1‑151 (MQ=255)
                                                                                                                     ggAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAg                               >  1:9275564/1‑151 (MQ=255)
                                                                                                                       aaaaGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCt                             >  1:1840796/1‑151 (MQ=255)
                                                                                                                       aaaaGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCt                             <  1:4674414/151‑1 (MQ=255)
                                                                                                                       aaaaGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCt                             >  1:4815800/1‑151 (MQ=255)
                                                                                                                         aaGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGt                           <  2:1102899/151‑1 (MQ=255)
                                                                                                                             cTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGAc                       >  1:10404661/1‑151 (MQ=255)
                                                                                                                             cTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGAc                       <  1:5719889/151‑1 (MQ=255)
                                                                                                                              tGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACg                      <  1:6823454/151‑1 (MQ=255)
                                                                                                                               ggggTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACgg                     <  1:5147584/151‑1 (MQ=255)
                                                                                                                               ggggTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACgg                     >  2:2483219/1‑151 (MQ=255)
                                                                                                                                 ggTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACGGcc                   <  1:10685125/151‑1 (MQ=255)
                                                                                                                                   tGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACGGCCAg                 <  2:10398495/151‑1 (MQ=255)
                                                                                                                                   tGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACGGCCAg                 <  2:1066139/151‑1 (MQ=255)
                                                                                                                                    ggATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACGGCCAgg                <  1:3257421/151‑1 (MQ=255)
                                                                                                                                      aTCGCTGAGCGTCAGTATTGGCGGGGGATAGGACGGAGTGCGAGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGAAAACGTCAGCCCGGGGGGGGAAGTCGGCGCAGGCTGCGCCAGGCGACGAAGCTGTAGACGGCCAGGGc              >  2:9284209/1‑151 (MQ=255)
                                                                                                                                        cGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACGGCCAGGGCGa            <  1:6084357/151‑1 (MQ=255)
                                                                                                                                        cGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACGGCCAGGGCGa            <  2:3810000/151‑1 (MQ=255)
                                                                                                                                              gCGTCAGTATTGGCGGGGGATCGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACGGCCAGGGCGAGCCAGa      >  2:7460457/1‑151 (MQ=255)
                                                                                                                                               cGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACGGCCAGGGCGAGCCAGAt     <  1:6583497/151‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                 
ACAGCGGTGTGCTCAAAGACTGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACGGCCAGGGCGAGCCAGATGAA  >  NC_008463/558900‑559195

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: