Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,267,804 G→A A231T (GCC→ACC)  tyrB → tyrosine aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,267,8040GA98.4% 208.2 / ‑4.5 64A231T (GCC→ACC) tyrBtyrosine aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (38/25);  total (39/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

CCCCTGGTATGACGAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCAGTGCTGAATGACGAGGCATTGAAAGCCAGCTGGCTGGCGGA  >  NC_000913/4267557‑4268045
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
ccccTGGTATGACGAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCa                                                                                                                                                                                                                                                <  1:593348/251‑1 (MQ=255)
       tATGACGAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCGACGATGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGACGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTCTCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTTGCACCATTGcc                                                                                                                                                                                                                                           >  1:648959/1‑249 (MQ=255)
        aTGACGAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTTTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAgc                                                                                                                                                                                                                                        >  1:178856/1‑251 (MQ=255)
                      aaCGGCGTGCGCTTTAATGACCTGTTGGCGACGCTGAAAACATTCCCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTAccc                                                                                                                                                                                                                            <  1:135388/249‑1 (MQ=255)
                       aCGGCGTGCGCTTTAATGACCGTTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGct                                                                                                                                                                                                                         <  2:461935/251‑1 (MQ=255)
                              gcgcTTTAATGACCTGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTg                                                                                                                                                                                                                   <  1:275179/250‑1 (MQ=255)
                                gcTTTAATGACCTGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGc                                                                                                                                                                                                                >  1:52413/1‑251 (MQ=255)
                                                   gACGCTGAAAACATTACCTGCCCGCAGGATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAAttt                                                                                                                                                                                             >  2:234070/1‑251 (MQ=255)
                                                    aCGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAAtttt                                                                                                                                                                                            >  1:16688/1‑251 (MQ=255)
                                                         gAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTccc                                                                                                                                                                                       <  1:45904/251‑1 (MQ=255)
                                                          aaaaCATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCCTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCt                                                                                                                                                                                      <  1:579011/251‑1 (MQ=255)
                                                            aaCATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCttt                                                                                                                                                                                    >  2:284833/1‑251 (MQ=255)
                                                            aaCATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCttt                                                                                                                                                                                    <  1:77183/251‑1 (MQ=255)
                                                              cATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTAc                                                                                                                                                                                  >  2:499153/1‑251 (MQ=255)
                                                                ttACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAACATTTTCTCCCTTTACgg                                                                                                                                                                                >  1:231956/1‑251 (MQ=255)
                                                                                         tgcATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGtt                                                                                                                                                       >  1:110265/1‑251 (MQ=255)
                                                                                          gcATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATGCCATTCCTCGATATTGCCTAGCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGGTCTCGAAAAGTTTCTCCCTGTACGGCGAGCGCGTCGGCGGACTTTCTGt                                                                                                                                                        <  2:141910/249‑1 (MQ=255)
                                                                                            aTCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTttatgttatg                                                                                                                                                    >  2:174256/1‑251 (MQ=255)
                                                                                              ccATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTAtgtg                                                                                                                                                  <  1:620159/251‑1 (MQ=255)
                                                                                                aTGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCa                                                                                                                                                                                                               >  1:455257/1‑188 (MQ=255)
                                                                                                aTGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCa                                                                                                                                                                                                               <  2:455257/188‑1 (MQ=255)
                                                                                                            cccAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAg                                                                                                                                    <  1:534984/251‑1 (MQ=255)
                                                                                                                 cGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCt                                                                                                                               <  1:229291/251‑1 (MQ=255)
                                                                                                                 cGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCt                                                                                                                               <  1:347570/251‑1 (MQ=255)
                                                                                                                 cGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGGGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATGTGGTGCCGGTATGGAAGAGGATGCCTACGCTCTTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGCCTTTCTGTTATGTGTGAAGATGCCGAAGCCGCt                                                                                                                               >  2:334334/1‑251 (MQ=255)
                                                                                                                  gggTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGCCTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTg                                                                                                                              >  2:399187/1‑251 (MQ=255)
                                                                                                                        cgacctcACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCcgcg                                                                                                                        >  2:157276/5‑251 (MQ=255)
                                                                                                                         gATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGt                                                                                                                       >  2:349138/1‑251 (MQ=255)
                                                                                                                                          gggggATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGGACTGGGGCAATTGaa                                                                                                        >  2:573893/3‑249 (MQ=255)
                                                                                                                                              gATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAAc                                                                                                  >  1:479536/1‑251 (MQ=255)
                                                                                                                                                tGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAAc                                                                                                  >  1:71931/1‑249 (MQ=255)
                                                                                                                                                 gCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAg                                                                                                <  1:143089/250‑1 (MQ=255)
                                                                                                                                                  cGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGtt                                                                                              >  2:5407/1‑251 (MQ=255)
                                                                                                                                                   ggTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTc                                                                                             <  2:9579/251‑1 (MQ=255)
                                                                                                                                                      gATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTcgcc                                                                                          >  2:571830/1‑251 (MQ=255)
                                                                                                                                                       aTTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTcgccg                                                                                         <  1:499153/251‑1 (MQ=255)
                                                                                                                                                                   aaaGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAg                                                                             <  2:110265/251‑1 (MQ=255)
                                                                                                                                                                     aGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGcc                                                                           >  1:41011/1‑251 (MQ=255)
                                                                                                                                                                      gCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGccc                                                                          <  1:82107/251‑1 (MQ=255)
                                                                                                                                                                        ccGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCccgc                                                                        >  1:509080/1‑251 (MQ=255)
                                                                                                                                                                        ccGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCccgc                                                                        >  1:21705/1‑251 (MQ=255)
                                                                                                                                                                         cgcgAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCTTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGCCTTTCTGTTATGTGTGAAGCTGCCGAAGCCGCTGGCCGCGTACTGGGGCAATGGACAGCAACAGTTCGCCGCACCTCCTCGAGAccgcc                                                                       >  2:452979/1‑251 (MQ=255)
                                                                                                                                                                         cgcgAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTAATGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGcccacc                                                                       >  1:481439/1‑248 (MQ=255)
                                                                                                                                                                               cTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGTCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAAGTACTCCAGCCCGCCGAAttt                                                                 >  2:476662/1‑251 (MQ=255)
                                                                                                                                                                                  attccattccTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTcg                                                                                            >  1:528687/1‑221 (MQ=255)
                                                                                                                                                                                  attccattccTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCGTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTcg                                                                                            <  2:528687/221‑1 (MQ=255)
                                                                                                                                                                                  attccattccTCGATATTGCCTATCAAGGATTTGGAGCCGGTATGGAAGAGGATGCATACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGACGATGCCGAAGCCGCTGGCCGCGTACTGGGGCCATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTgg                                                              >  1:396411/1‑251 (MQ=255)
                                                                                                                                                                                     ccattacTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTgg                                                           >  2:281008/1‑250 (MQ=255)
                                                                                                                                                                                          ccTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAg                                                       <  1:284833/250‑1 (MQ=255)
                                                                                                                                                                                           cTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAggt                                                     >  1:397667/1‑251 (MQ=255)
                                                                                                                                                                                           cTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAggt                                                     >  1:494805/1‑251 (MQ=255)
                                                                                                                                                                                            tCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGACGAGGATGCCTACGCTATTCGCACCATTGCAAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGATCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGGgcgaanggg                                                    >  2:597473/1‑245 (MQ=255)
                                                                                                                                                                                                            aGGATTGGGAGCCGGTATGGAAGAGGATGCCTACGCTATTCGAATAATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCCGCCCGCCGCATTTTGGTGCGCAGGTGGTGGCTGCAGTGCTg                                     >  2:16302/1‑250 (MQ=255)
                                                                                                                                                                                                             ggATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCAGTGCTGa                                    >  2:641835/1‑250 (MQ=255)
                                                                                                                                                                                                               aTTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCAGTGCTGAATg                                 <  1:71908/251‑1 (MQ=255)
                                                                                                                                                                                                                 ttGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCAGTGCTGAATGAc                               >  1:500869/1‑251 (MQ=255)
                                                                                                                                                                                                                        cGGTATGGAAGAGGATGCCTACGCTATTCGCACCAGTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTTCGCAGGTGGTGGCTGCAGTGCTGAATGACGAGGcct                        <  2:648959/251‑3 (MQ=255)
                                                                                                                                                                                                                           tATGGAAGAGGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCAGTGCTGAATGACGAGGCATTg                      >  1:106572/1‑250 (MQ=255)
                                                                                                                                                                                                                                aaGAGGATGCCTACGCTATTCGAACCATTGCCAGCGCTGGATTACACGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGGCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGGAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGAAGGAGGTGGCTGCAGAGCTGAATGACGAGGCATTGAAAGcc                >  2:472009/1‑251 (MQ=255)
                                                                                                                                                                                                                                 agagGATGCCTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCAGTGCTGAATGACGAGGCATTGAAAGcc                >  2:155499/1‑250 (MQ=255)
                                                                                                                                                                                                                                         ccTACGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCAGTGCTGAATGACGAGGCATTGAAAGCCAgctggctg       <  1:174256/251‑1 (MQ=255)
                                                                                                                                                                                                                                            aCGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCAGTGCTGAATGACGAGGCATTGAAAGCCAGCtggctggc     <  1:571830/250‑1 (MQ=255)
                                                                                                                                                                                                                                             cGCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCAGTGCTGAATGACGAGGCATTGAAAGCCAGCTGGCTGGCgg   >  1:620998/1‑251 (MQ=255)
                                                                                                                                                                                                                                              gCTATTCGCACCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCAGTGCTGAATGACGAGGCATTGAAAGCCAGCTGGCTGGCGGa  <  2:71931/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
CCCCTGGTATGACGAAGCGACTAACGGCGTGCGCTTTAATGACCTGTTGGCGACGCTGAAAACATTACCTGCCCGCAGTATTGTGTTGCTGCATCCATGTTGCCACAACCCAACGGGTGCCGATCTCACTAATGATCAGTGGGATGCGGTGATTGAAATTCTCAAAGCCCGCGAGCTTATTCCATTCCTCGATATTGCCTATCAAGGATTTGGTGCCGGTATGGAAGAGGATGCCTACGCTATTCGCGCCATTGCCAGCGCTGGATTACCCGCTCTGGTGAGCAATTCGTTCTCGAAAATTTTCTCCCTTTACGGCGAGCGCGTCGGCGGACTTTCTGTTATGTGTGAAGATGCCGAAGCCGCTGGCCGCGTACTGGGGCAATTGAAAGCAACAGTTCGCCGCAACTACTCCAGCCCGCCGAATTTTGGTGCGCAGGTGGTGGCTGCAGTGCTGAATGACGAGGCATTGAAAGCCAGCTGGCTGGCGGA  >  NC_000913/4267557‑4268045

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: