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.T97.9% 1592.0 / 19.5 388intergenic (+618/+71)PA14_RS02545/rarDglycine cleavage system protein R/EamA family transporter RarD
Reads supporting (aligned to +/- strand):  ref base . (7/1);  new base T (198/178);  total (207/181)
Fisher's exact test for biased strand distribution p-value = 7.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 406 total aligned reads displayed.

ACAGCGGTGTGCTCAAAGACTGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACGGCCAGGGCGAGCCAGATGAA  >  NC_008463/558900‑559195
                                                                                                                                                         |                                                                                                                                                  
aCAGCGGTGTGCTCAAAGACTGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTa                                                                                                                                                     <  1:6075145/151‑1 (MQ=255)
   gCGGTGTGCTCAAAGACTGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTAttg                                                                                                                                                 <  1:2839854/151‑3 (MQ=255)
   gCGGTGTGCTCAAAGACTGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTAttg                                                                                                                                                 <  2:7085773/151‑3 (MQ=255)
        gtgCTCAAAGACTGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCggg                                                                                                                                            >  2:5113854/1‑151 (MQ=255)
            tCAAAGACTGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGAt                                                                                                                                        >  1:2661540/1‑151 (MQ=255)
            tCAAAGACTGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGGCAGTCTTGGCGGGGGAt                                                                                                                                        >  2:7952556/1‑151 (MQ=255)
             cAAAGACTGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATa                                                                                                                                       >  2:918080/1‑151 (MQ=255)
               aaGACTGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAGAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATagg                                                                                                                                     >  2:1808829/1‑151 (MQ=255)
                  aCTGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATaggaag                                                                                                                                  >  1:11379177/1‑151 (MQ=255)
                   cTGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATaggaagg                                                                                                                                 <  2:5002464/151‑1 (MQ=255)
                   cTGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATaggaagg                                                                                                                                 >  1:612115/1‑151 (MQ=255)
                    tGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATaggaagga                                                                                                                                >  2:3175936/1‑151 (MQ=255)
                    tGATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATaggaagga                                                                                                                                <  1:8649472/151‑1 (MQ=255)
                     gATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAg                                                                                                                               <  1:1480784/151‑1 (MQ=255)
                     gATCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAg                                                                                                                               <  2:10796514/151‑1 (MQ=255)
                      aTCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGt                                                                                                                              <  2:10722401/151‑1 (MQ=255)
                       tCAGGCGCTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTg                                                                                                                             >  2:7197254/1‑151 (MQ=255)
                              cTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGttt                                                                                                                      >  1:7826392/1‑151 (MQ=255)
                              cTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGttt                                                                                                                      >  1:8148594/1‑151 (MQ=255)
                              cTTTGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGttt                                                                                                                      >  1:9928179/1‑151 (MQ=255)
                                ttGCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGttttt                                                                                                                    <  2:2535238/151‑1 (MQ=255)
                                  gCAGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTc                                                                                                                  >  2:9575667/1‑151 (MQ=255)
                                    aGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGc                                                                                                                >  2:10667818/1‑151 (MQ=255)
                                    aGGGCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGc                                                                                                                >  1:10778430/1‑151 (MQ=255)
                                     gggCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGcc                                                                                                               <  1:8569611/151‑1 (MQ=255)
                                     gggCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGcc                                                                                                               <  1:8818146/151‑1 (MQ=255)
                                       gCTTCTTCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAg                                                                                                             <  1:6700142/151‑1 (MQ=255)
                                        cttcttCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAgg                                                                                                            <  2:2152761/151‑1 (MQ=255)
                                        cttcttCCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAgg                                                                                                            <  1:7519644/151‑1 (MQ=255)
                                          tcttcCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCt                                                                                                          <  1:7418082/151‑1 (MQ=255)
                                           cttcCGACGAGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTAGTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCtt                                                                                                         >  1:2526101/1‑151 (MQ=255)
                                               cgacgaGTTATCCACGCGCCTGCGGGGCGGAGAAGGCAGGTCTATGCACAAAGTACGTGCCTGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCgg                                                                                                     >  2:8258226/1‑151 (MQ=255)
                                                gacgaGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGGGCGTGTTTTTTCGCCAGGCTTGCGGa                                                                                                    >  2:624014/1‑151 (MQ=255)
                                                gacgaGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGGATGGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGa                                                                                                    >  2:5241283/1‑151 (MQ=255)
                                                    aGTTATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTc                                                                                                >  2:4312296/1‑151 (MQ=255)
                                                      ttATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCg                                                                                              <  2:9086195/151‑1 (MQ=255)
                                                       tATCCACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGc                                                                                             <  2:11433651/151‑1 (MQ=255)
                                                           cACGCGCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGc                                                                                         <  1:1985124/151‑1 (MQ=255)
                                                             cgcgCCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCt                                                                                       <  2:3025325/151‑1 (MQ=255)
                                                              gcgcCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTg                                                                                      <  2:6953911/151‑1 (MQ=255)
                                                              gcgcCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGGATTGGCGGGGGATAGGAGGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTg                                                                                      >  1:4471497/1‑151 (MQ=255)
                                                               cgcCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTgg                                                                                     >  1:5544180/1‑151 (MQ=255)
                                                               cgcCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGGGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTng                                                                                     <  2:3899834/151‑1 (MQ=255)
                                                                gcCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTggg                                                                                    >  1:1251880/1‑151 (MQ=255)
                                                                gcCTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTggg                                                                                    <  2:3118618/151‑1 (MQ=255)
                                                                 ccTGCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTgggg                                                                                   >  1:10300720/1‑151 (MQ=255)
                                                                    gCGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATa                                                                                <  1:9228654/151‑1 (MQ=255)
                                                                     cGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATaa                                                                               >  2:7662242/1‑151 (MQ=255)
                                                                     cGGGGAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATaa                                                                               <  2:10137600/151‑1 (MQ=255)
                                                                      ggggAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGCACCTGTGGAATAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGATGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAAc                                                                              >  2:11196909/1‑151 (MQ=255)
                                                                      ggggAGGAGAAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAAc                                                                              >  2:4068305/1‑151 (MQ=255)
                                                                            gagaAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGc                                                                        >  2:7475920/1‑151 (MQ=255)
                                                                            gagaAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCGGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGAGAGGTAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATCACGTCAGc                                                                        >  2:4635213/1‑151 (MQ=255)
                                                                             agaAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGcc                                                                       >  1:5876400/1‑151 (MQ=255)
                                                                              gaAGGCAGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGccc                                                                      <  1:4071960/151‑1 (MQ=255)
                                                                                    aGGTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGCGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTgg                                                                <  2:7820313/151‑1 (MQ=255)
                                                                                     ggTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTggg                                                               >  2:5860263/1‑151 (MQ=255)
                                                                                      gTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTgggg                                                              >  1:164337/1‑151 (MQ=255)
                                                                                      gTCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTgggg                                                              >  1:8020882/1‑151 (MQ=255)
                                                                                       tCTATGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGa                                                             <  1:962743/151‑1 (MQ=255)
                                                                                          aTGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAgg                                                          <  2:10906997/151‑1 (MQ=255)
                                                                                           tGCACACAGTACGTGCATGAACCTGTGGAAACGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCAGGCTGGGGAAggc                                                         <  2:2680228/151‑1 (MQ=255)
                                                                                           tGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAggc                                                         <  2:2919306/151‑1 (MQ=255)
                                                                                           tGCACAAAGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAggc                                                         >  1:2185019/1‑151 (MQ=255)
                                                                                                 aaGTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCa                                                   >  2:4381240/1‑151 (MQ=255)
                                                                                                  aGTACGTGCAAGAACCTGTGGAAAAGACTGGGGTGGATCGCTGAGCGTCACGATTGGCGGGGGATCGCAAGCAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGATCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGAAg                                                  >  1:11276745/1‑151 (MQ=255)
                                                                                                   gTACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAgg                                                 >  1:10949877/1‑151 (MQ=255)
                                                                                                    tACGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGc                                                <  1:6310095/151‑1 (MQ=255)
                                                                                                     aCGTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCt                                               >  2:3709390/1‑151 (MQ=255)
                                                                                                       gTGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTgc                                             <  2:9278730/151‑1 (MQ=255)
                                                                                                        tGCATGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTgcg                                            <  1:10454476/151‑1 (MQ=255)
                                                                                                           aTGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGGATTGGCGGGGGATCGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCa                                         <  2:5589196/151‑1 (MQ=255)
                                                                                                            tGAACCTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAg                                        >  1:898732/1‑151 (MQ=255)
                                                                                                              aaCCTGTGGAAACGCCTGGGGTGGATCGCTGAGCGTCAGGATTGGCGGGGGATAGGTCGGAGGGCGGGTTTTTTCGCCAGGCTTGCGGAGTGCCGCAAGCCTGGGGGTGACGTCAGCCCGGCGGGGGAATGCGGCGCAGGCTGCGCCAgtc                                      >  2:6768529/1‑149 (MQ=255)
                                                                                                                ccTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGt                                    >  1:3625309/1‑151 (MQ=255)
                                                                                                                 cTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTc                                   <  1:10150175/151‑1 (MQ=255)
                                                                                                                 cTGTGGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTc                                   >  1:1114060/1‑151 (MQ=255)
                                                                                                                  tgtgGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCg                                  <  1:3426996/151‑1 (MQ=255)
                                                                                                                   gtgGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCCGGCGTCGa                                 >  1:7148014/1‑151 (MQ=255)
                                                                                                                    tgGAAAAGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGaa                                >  1:181912/1‑151 (MQ=255)
                                                                                                                       aaaaGCCTGGGGTGGATCGCTGAGCGGCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCACGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCt                             >  2:4855144/1‑151 (MQ=255)
                                                                                                                        aaaGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTg                            >  1:5423100/1‑151 (MQ=255)
                                                                                                                        aaaGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTg                            >  2:10515395/1‑151 (MQ=255)
                                                                                                                          aGCCTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGATTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTa                          <  2:4740840/151‑1 (MQ=255)
                                                                                                                            ccTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGa                        >  1:84099/1‑151 (MQ=255)
                                                                                                                            ccTGGGGTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGa                        >  1:4280206/1‑151 (MQ=255)
                                                                                                                               ggggTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACgg                     <  1:9303765/151‑1 (MQ=255)
                                                                                                                                 ggTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACGGcc                   >  1:5093503/1‑151 (MQ=255)
                                                                                                                                  gTGGATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACGGCCa                  >  2:5386541/1‑151 (MQ=255)
                                                                                                                                    ggATCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACGGCCAgg                >  2:3389247/1‑151 (MQ=255)
                                                                                                                                      aTCGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACGGCCAGGGc              <  2:94244/151‑1 (MQ=255)
                                                                                                                                        cGCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACGGCCAGGGCGa            >  2:8734705/1‑151 (MQ=255)
                                                                                                                                        cGCTGAGCGTCAGGATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCGTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGCAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACGGCCAGGGCGa            <  2:9013922/151‑1 (MQ=255)
                                                                                                                                         gCTGAGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACGGCCAGGGCGAg           <  1:3107485/151‑1 (MQ=255)
                                                                                                                                             aGCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACGGCCAGGGCGAGCCAg       >  2:7339511/1‑151 (MQ=255)
                                                                                                                                              gCGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACGGCCAGGGCGAGCCAGa      >  1:2041170/1‑151 (MQ=255)
                                                                                                                                               cGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGTGGATAACGTCAGCACGGCTGGGGAATGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACGGCCAGGGCGAGCCAGAt     >  2:1422415/1‑151 (MQ=255)
                                                                                                                                               cGTCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACGGCCAGGGCGAGCCAGAt     <  2:4198692/151‑1 (MQ=255)
                                                                                                                                                 tCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACGGCCAGGGCGAGCCAGATGa   <  1:9473766/151‑1 (MQ=255)
                                                                                                                                                 tCAGTATTGGCGGGGGATAGGAAGGAGTGCGTGTTTTTTCGCCAGGCTTGCGGAGTTCCGCAAGCCTGGGGATAACGTCAGCCCGGCTGGGGAAGGCGGCGCAGGCTGCGCCAGGCGTCGAAGCTGTAGACGGCCAGGGCGAGCCAGATGa   >  1:6938497/1‑151 (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: