Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,430,345 A→G 13.9% E161G (GAG→GGG)  cycA → D‑alanine/D‑serine/glycine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,430,3450AG13.9% 174.5 / 2.5 76E161G (GAG→GGG) cycAD‑alanine/D‑serine/glycine transporter
Reads supporting (aligned to +/- strand):  ref base A (25/38);  new base G (10/2);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 9.46e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAA  >  NC_000913/4430118‑4430579
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
aTCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTTCTGACGCTCCATCTCGCCCCCCGGAAAATGTTCGGGGAGATGGAGTTCTGGtttt                                                                                                                                                                                                                                        >  1:909911/1‑238 (MQ=255)
   tttttgtttcTTTGCTTCCGTTTTTCTTGGGCTGTGGGGATGTTATTTCCCCGGCTGGATTTCATGGTTGTGGTGGGTTGTAACCGGTATGGCAGACGGGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCga                                                                                                                                                                                                                                     <  2:390545/230‑1 (MQ=255)
    ttttGTGGCTTCGGTTTCTGTTTTCTTGGGCCGGTGGTAGGATTTTTCACCGGTGGGTCTTTCTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCgat                                                                                                                                                                                                                                    <  2:216166/239‑1 (MQ=255)
          gACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCGTCGCTGGCGGTGTTAGGGCGTCTGCCGCCGCTCCAATTCCGCCCCGGGGAAAATTGCGGTGGGGGGGAGGTTTGGTTTGGCATGAACaa                                                                                                                                                                                                                             >  2:535596/1‑239 (MQ=255)
            cTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGCGGCGATCACGGCTTATGCTCAGTTCTGGTTACCCGATCTCTCGGCCGGGGTGGCCTCGCTGGCGGGGAGTGTGCTCGTGATGACGCCCAATCTCGGCCCCGGGAAAAAGTTCGGTGGGGTGGGGGCTGGGTTTGGGTTGAACaaaa                                                                                                                                                                                                                           >  2:163608/1‑239 (MQ=255)
             ttCGCTTCCGATTTGCTCGGGCCGTGGGCAGGATATTTCCCCGGCTGGGCTTACTGGTTCTGCTGGGTTGTAACCGGTAGGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAAt                                                                                                                                                                                                                          <  1:81013/239‑1 (MQ=255)
                  ttCCGGTTTCCTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGGTCTGCTGGGCTGTAACCGGTATGGCAAACGTGGTGGCGATCAAGGCTTAAGCTACGGTCTGGTTCCCCGAACTCTCCGGCTGGGTCGCCTCGCAGGCGGGGGTATGGCCGCTGACGACGCCCCACATCGGCCCCGCGAAAATGTTCCGTTCGATGGGGTTCCTGTTTGGGATGATCaacaccagac                                                                                                                                                                                                                     >  2:569050/1‑231 (MQ=255)
                   tCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGc                                                                                                                                                                                                                    <  1:12137/239‑1 (MQ=255)
                     cgttttACCCGGGCCTTGGGCAGGATATTTCCCCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCa                                                                                                                                                                                                                  <  1:201144/236‑1 (MQ=255)
                       ttttACCGGGGCCGTGGGCATGATATTTCACCGGGTGTACTTACTGTTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCAtt                                                                                                                                                                                                                <  1:594930/238‑1 (MQ=255)
                                 gCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGCTCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGGTAGGGCTGCCGCTTACGCCCAATTCCGCCCCCCTGAAAAAGTTTGGGGGGATGGGGGTTTGGTTTTCAATGTAACAAAACGTCGGCATTGGGTTGcggg                                                                                                                                                                                                      >  2:874192/1‑236 (MQ=255)
                                 gCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGGGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGGGGTGGGGGTCTGGGTTTCGATGATCAAAATACTCCCCAtttggtgcctgg                                                                                                                                                                                                      >  1:651863/1‑229 (MQ=255)
                                       ggcaggtgttttCCCCGGTTGGATTTATTGGTTCGTGTTGGTTTTAACCTGTATGGCAGGCGGGGGGGCGATCACGGCTTATGCCCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTgtcg                                                                                                                                                                                                <  2:383214/230‑1 (MQ=255)
                                         gaGGATTTTTCCCCGGGTGGACTTCGTGGTTCTGTTTGGTTGTAACCGTGATGGCAGACGTGGTGGCGTTCACGGCTTATGCCTAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTgtcgt                                                                                                                                                                                               <  2:245993/237‑1 (MQ=255)
                                                      cGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCg                                                                                                                                                                                                                                      <  2:93113/187‑1 (MQ=255)
                                                      cGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCg                                                                                                                                                                                                                                      >  1:93113/1‑187 (MQ=255)
                                                                        ctgctgGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGc                                                                                                                                                                                           <  2:671153/211‑1 (MQ=255)
                                                                        ctgctgGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGc                                                                                                                                                                                           >  1:671153/1‑211 (MQ=255)
                                                                               gTTGTAACCGGTAGGGCAGACGGGGGGGTGATCGCGGGTTGTGCTCAGTTCGGGTTCCCCGAGCTCCCCGGGTGGGTCGCCTCGTGGGGGGGGATAGTGCTGGTGCTGACGCTCAATCGCGCACCCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTCTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGAc                                                                                                                                                       <  2:545094/239‑1 (MQ=255)
                                                                                ttttAACCGGTTTGGCGGGCGTGGTGGGCGTCACGGTTTTTTCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGGTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGAc                                                                                                                                                       <  2:801941/235‑1 (MQ=255)
                                                                                tttgaaCCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATTCTCAGTTCTGGTTCCCCGATCCCGCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGAc                                                                                                                                                       <  1:200142/234‑1 (MQ=255)
                                                                                ttGTAACCGGTATGGCAGGCGTGGTTGCCATCACGGCTTATGCCCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGAc                                                                                                                                                       <  1:733578/238‑1 (MQ=255)
                                                                                      ccGGTATGGCAGACGTTGTTGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGtgtg                                                                                                                                                <  1:489440/239‑1 (MQ=255)
                                                                                        ggTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGGGATGGGGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGGGGGGATGCACTTTCAGTCCCCGACTGGTGTGGa                                                                                                                                              >  1:525181/1‑239 (MQ=255)
                                                                                                     gtggggctgtCACGTTTTATTCTCAGTTTGGGTCCCCCGCTCTCCCCGGCTGGGGCGCCTGGGTGGCGGTGATGTGGCTGCTGCTGACCGTCAATCCCGCCACCGTGAAAGTGTTCGGTGAGATGGAGTTCTGGTTTGGGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCCTTTTAAGTCCCCGACTGGTGTGGAAGCGTCATTcgcg                                                                                                                                <  2:620452/230‑1 (MQ=255)
                                                                                                      ggtggCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGCCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTcgcg                                                                                                                                <  1:6366/238‑1 (MQ=255)
                                                                                                         ggtgATCCCGGCTTTTGCGCAGTTCTGGTTCCCCGATCTCGCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCAt                                                                                                                            <  2:909911/235‑1 (MQ=255)
                                                                                                         ggCGATAACGGTTTATGCTCAGTTCTGGGTCCCCGAGCCCTCCGGCTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAAGGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCAt                                                                                                                            <  1:662232/238‑1 (MQ=255)
                                                                                                              tCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGtt                                                                                                                                                                                                                                                >  1:50776/1‑120 (MQ=255)
                                                                                                              tCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGtt                                                                                                                                                                                                                                                <  2:50776/120‑1 (MQ=255)
                                                                                                                       aTGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGGGGGATGGGGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGGGGCGGTGCACTTTCAGTCCACGACTGGTGTGGAAGCGGCATTTGCGGATTTTGGGAAGGAcggg                                                                                                            >  1:368015/1‑238 (MQ=255)
                                                                                                                               ttCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCCCTTTCCGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGAcggcggtcggt                                                                                                    >  1:306748/1‑233 (MQ=255)
                                                                                                                               ttCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCCCTTTCAGTGACCGACTGGTGTGGAAGCGTCATTCACGCATTTTTGGAATAAcggcgggcggt                                                                                                    >  1:142398/1‑233 (MQ=255)
                                                                                                                               ttCTGGTTCCCCGACCTGTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCACGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGa                                                                                                                 <  1:586292/227‑1 (MQ=255)
                                                                                                                                    gTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATACTGCTGCTGCTGACCCTCAATCTCGCCCCCGTGAAAAAGTTTGGGTTGATGGGGGTCTTGGTTGCGATGATCAAAAACGCGGCCCGTTGGTCGCTGATTGTCTGCGGCCTCGGCAAGGGTGCGCTTGCCCTCCCATCCCCTACTGGGGTGGAAGGGCCTTTCGCGCATTTTGGGGAATGCTGCCGGCGGGTcccca                                                                                              >  2:833624/1‑238 (MQ=255)
                                                                                                                                     ttCCCCGTTCTCCCCGAGTGGCGCGCCGCGCTGGCGGTGGTAGGGGTTCGGCTGACGCTCAATCCCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGGCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGa                                                                                             <  2:525181/238‑1 (MQ=255)
                                                                                                                                       ccccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGGGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGGCATGGTGGCGATGCACTTTCATTCACCGACTGGTGTGGAAACGTCATTTGGGCCTTTTTGGAATGACGGCGGCCGGGTTCCGGAAg                                                                                          >  1:149246/1‑239 (MQ=255)
                                                                                                                                         ccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGAGGGGGGTCTGGTTTGCGATGATCCAAATCGTCGCCATTGTGTCGCTGGTTTTTGTCGGCCTGTTCATGGGGGCGGTGGACTTTCAGTGCCCCGCTGGTGGGGGAGCCGTCTTCCCGGCCTTTTGGAAGGACGGGCGGCGGGTTCCCGAAggg                                                                                         >  2:4639/1‑239 (MQ=255)
                                                                                                                                         ccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGGGCGATGGGGGTCTGGTTTGCGATGATCAAAATCGTCGCCCTTGGTTCCCTGATTGTCGTCCGCCTGGTCAAGGTGGCGATGCACTTTCAGTCACCCACCTGTTGGGGAGCGCCATTCGCGCACTTTTGGGATGAAGGCGGGCGGGTCCCCAAAAgg                                                                                         >  2:52021/1‑239 (MQ=255)
                                                                                                                                         ccGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCCCAATCTCGCCCCCGTGAAAATGTTCGGTGAGATGGGGTTCTGGTTTGCGATGAGCAAAATCTGCGCCATTGTGTCGCTGAGTGTCGTCGGCCTGGTCCATGTGGCGGTGGACTTTCAGGCCCCACCCGGGTTGGAAGGGTCCATCGCGCCTTTTTGGTATGaagctcccacc                                                                                                       >  2:586292/1‑217 (MQ=255)
                                                                                                                                                           tccccTCGCTTGCGGTTATAGTTGTTTTTCTGGCGCTCAATCTGGCCACCGTTAAAAGGTGCGGTGGAGTGGAGTTTGGGTTTGCGATGATCAAAATCGTCGTCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGcc                                                                      <  2:396152/236‑1 (MQ=255)
                                                                                                                                                           tCGCCTCGCTGGCGGTGTTAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTTCGCTGATCAAAATCGTCGCCATTGTGTCGCTGATTGTGGTCGGCCTGGTCATAGTGGGGATGCACTTTCGGTCCACGACGGCGGTGGAAGAGGCAGTCTCGCGTTTGTGGACAGAAGGCCGCGGGGACCCGCAAGGTTTTAGGGGGTTTTTTTcc                                                                      >  2:190071/1‑239 (MQ=255)
                                                                                                                                                           tCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCCCTTTCAGTCACCGACTGGGGTGGAAGCGTCATTTGCGCATTTTGGGAATGACGGCGGCTGGTTCCCGGAAGGTTTAAATGGCTTTTTTGcc                                                                      >  1:188319/1‑239 (MQ=255)
                                                                                                                                                            cGCCTCGCTGGCGGTGATAGTTCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGcc                                                                      <  1:87261/238‑1 (MQ=255)
                                                                                                                                                               cTCGCTGGCGGTGATAGTGCTGCTGCTGAAGCTCAATCTCGCCACCGTGAAAGGGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAACGGTCATTCGCGCATTTTTGGGATGACGGCGGCTGGGTTCCCGAAGGGTTAAATGGCGTCTTTTCCgggt                                                                  >  1:518668/1‑237 (MQ=255)
                                                                                                                                                               cTCGCGGGCGGTGATAGTGCTGCTGCTTACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGGCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGa                                                                   <  1:190071/238‑1 (MQ=255)
                                                                                                                                                                      gCGGTGATAGTGCTGCTTCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACttt                                                                                                                                                                  <  1:606276/141‑1 (MQ=255)
                                                                                                                                                                        ggTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTTCCGGGTTCCAGATAg                                                         >  1:257741/1‑239 (MQ=255)
                                                                                                                                                                               gtgctggtgctgGCTCTTAATTTCGCTCCCGGGAAAATTGTCGGTTAGGTGGAGTTCTGGTTTGCGATTATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGtttt                                                  <  2:306748/239‑1 (MQ=255)
                                                                                                                                                                                tgctgctgcTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTAATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGTTTCTTTCCCGGATTCCAGATAGCGGattt                                                  >  1:122327/1‑238 (MQ=255)
                                                                                                                                                                                tgctgctgcTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGATTGCGATGATCAAAATCGTCGCCATTGTGCCGCTGATTGTCGTCGGCCTGGTCATGAGGGCGGTGCACTTTCAGTCACCGa                                                                                                                                                        >  2:606276/1‑141 (MQ=255)
                                                                                                                                                                                     ctgctgACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGCGGGCGGCTGGTTCCCGAAAGGTTTAATGGGCTTTTTTGCCGGATTCCAAATAGCGGTTTTTGGttt                                            >  2:60541/1‑239 (MQ=255)
                                                                                                                                                                                      tgctgACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAGGGGGGTTCTTTTCCGGGTTCCAAGTAGGggttttggttttt                                           >  2:776650/1‑238 (MQ=255)
                                                                                                                                                                                          gACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCCGATTGCGGTTTTCGCTTTCGTgg                                       >  1:825773/1‑239 (MQ=255)
                                                                                                                                                                                          gACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTgg                                       >  1:125101/1‑239 (MQ=255)
                                                                                                                                                                                           aCGCCCAATCGCGCCACCGTGAAAATGGTCGGTTAGGTGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGt                                         <  1:569050/236‑1 (MQ=255)
                                                                                                                                                                                              gtccAATCCCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGa                                    <  1:684039/236‑1 (MQ=255)
                                                                                                                                                                                               tCAATGTCGCCACCGTGAATATGTTCGGTGAGGTGGGGTTCTGGTTTGGGATTATCAAAAGCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCTCCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGTGAtt                                  <  2:758451/239‑1 (MQ=255)
                                                                                                                                                                                                 cAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGCCTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGGTTCCAGGTAGCGGGTTTCGCTTTCGGGGGGGttt                                 >  2:703481/1‑238 (MQ=255)
                                                                                                                                                                                                   aTCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGGGCTACCGTGATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAGTGGGCTTATTTGCCGGATTCCAAGTAGCGGTTTTCGCTTTCGGGGGGAtgtag                               >  2:84836/1‑235 (MQ=255)
                                                                                                                                                                                                        gacACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggt                          <  1:193069/237‑1 (MQ=255)
                                                                                                                                                                                                        gCCACCGTGAAAATGTTCTGGGGGTGGGGGTTCGGTGTTGCGATGATCAAAATCGTCGCCATTGGGTCGCTGATTGGCGTTGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGGCTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggt                          <  2:651863/239‑1 (MQ=255)
                                                                                                                                                                                                        gCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGCAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATACCGGTTTTCGCTTTCGTGGGGATTGAGCTggt                          >  1:424384/1‑239 (MQ=255)
                                                                                                                                                                                                        gCCACCGTGAAAATGGTCGGTGATATGGAGTTCCGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCCGTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggt                          <  1:72406/239‑1 (MQ=255)
                                                                                                                                                                                                        gCCACCGTGAAAAGGGTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggt                          <  1:4639/239‑1 (MQ=255)
                                                                                                                                                                                                         ccACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggta                         <  1:703481/239‑1 (MQ=255)
                                                                                                                                                                                                           cggttGAAAATTTTCGGTGGGATGGAGTTCTGGTTTGCGATGATCAAAATCGGCGACATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTggta                         <  2:485323/233‑1 (MQ=255)
                                                                                                                                                                                                            ccGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCTAGATAGCGGTTTTCGCTTTCTTGGGGTTTGAGCTgggaggt                      >  1:81541/1‑239 (MQ=255)
                                                                                                                                                                                                            ccGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTTGTGGGGGTTGAGCGggtgggg                      >  2:489796/1‑238 (MQ=255)
                                                                                                                                                                                                                         cGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTTCCGGATTCCAGATAGCGGGTTTCGCTTTCGTGGGGATTGGGCTGGTAGGTACAAAAg               >  1:485323/1‑233 (MQ=255)
                                                                                                                                                                                                                          ggTGGGATTGTGTTTTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGa        <  2:257741/239‑1 (MQ=255)
                                                                                                                                                                                                                          ggTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAg                               <  1:311001/216‑1 (MQ=255)
                                                                                                                                                                                                                                aTGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAg                               <  2:670764/210‑1 (MQ=255)
                                                                                                                                                                                                                                aTGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAg                               >  1:670764/1‑210 (MQ=255)
                                                                                                                                                                                                                                aTGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCa   <  1:535596/238‑1 (MQ=255)
                                                                                                                                                                                                                                 tGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCaa  <  1:60541/238‑1 (MQ=255)
                                                                                                                                                                                                                                    aGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGGTTAAGGGGCTTTTTTTCCGGATTCCAGATAACGGGTTTTGCTTTCGTGGGGATTGGGCTggtaggta                     >  2:311001/1‑216 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
ATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTTCTGGTTTGCGATGATCAAAATCGTCGCCATTGTGTCGCTGATTGTCGTCGGCCTGGTCATGGTGGCGATGCACTTTCAGTCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTACAACAGCTGCGGAAACCAA  >  NC_000913/4430118‑4430579

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