Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,980,727 T→C 6.4% intergenic (+31/‑160) wecG → / → yifK UDP‑N‑acetyl‑D‑mannosaminuronic acid transferase/putative APC family amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,980,7270TC 6.4% 298.2 / 8.4 93intergenic (+31/‑160)wecG/yifKUDP‑N‑acetyl‑D‑mannosaminuronic acid transferase/putative APC family amino acid transporter
Reads supporting (aligned to +/- strand):  ref base T (48/39);  new base C (4/2);  total (52/41)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

GCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAA  >  NC_000913/3980585‑3980875
                                                                                                                                              |                                                                                                                                                    
gCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTcc                                                                                                                                                >  3:344399/1‑149 (MQ=255)
 cAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCg                                                                                                                                               >  1:428945/1‑149 (MQ=255)
   aaCGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGACCACTCCGTc                                                                                                                                             <  4:463227/149‑1 (MQ=255)
    aCGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCg                                                                                                                                            >  1:146070/1‑149 (MQ=255)
         ggggCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGt                                                                                                                                       <  1:606814/149‑1 (MQ=255)
               ggAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCtt                                                                                                                                 <  1:126038/149‑1 (MQ=255)
               ggAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCtt                                                                                                                                 >  4:194710/1‑149 (MQ=255)
                 aGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCtttt                                                                                                                               >  1:189951/1‑149 (MQ=255)
                  gTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCttttt                                                                                                                              <  1:37319/149‑1 (MQ=255)
                    ggCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTtta                                                                                                                            <  2:57987/149‑1 (MQ=255)
                     gCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTttat                                                                                                                           >  3:659946/1‑149 (MQ=255)
                       tctACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTttatta                                                                                                                         >  4:131099/1‑149 (MQ=255)
                       tctACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTttatta                                                                                                                         >  3:213766/1‑149 (MQ=255)
                       tctACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTttatta                                                                                                                         >  4:719606/1‑149 (MQ=255)
                         tACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTaca                                                                                                                       <  3:508579/149‑1 (MQ=255)
                            cGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACaa                                                                                                                    >  2:117481/1‑149 (MQ=255)
                             gCCTGCTTTCGGAGCCGAGCCGCATTAAGCATCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAGAGCGACGTAGTGGTCACTACGTCGCCTGTCTCCCTTTTTATTACACaaa                                                                                                                   <  4:468859/149‑1 (MQ=255)
                             gCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACaaa                                                                                                                   >  2:604540/1‑149 (MQ=255)
                              ccTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTTTCTGCTTTTTTATTACACAAAg                                                                                                                  >  3:465486/1‑149 (MQ=255)
                              ccTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAg                                                                                                                  >  4:592740/1‑149 (MQ=255)
                                tGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCa                                                                                                                >  2:266107/1‑149 (MQ=255)
                                  cTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCAtt                                                                                                              <  2:279858/149‑1 (MQ=255)
                                  cTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCAtt                                                                                                              >  2:497242/1‑149 (MQ=255)
                                       gCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAAt                                                                                                         >  4:655842/1‑149 (MQ=255)
                                         agccgagccgCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAAttt                                                                                                       >  1:355140/1‑149 (MQ=255)
                                          gccgagccgCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAAtttt                                                                                                      >  2:65545/1‑149 (MQ=255)
                                           ccgagccgCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAAttttt                                                                                                     >  3:735947/1‑149 (MQ=255)
                                           ccgagccgCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAAttttt                                                                                                     >  1:428107/1‑149 (MQ=255)
                                                ccgCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGc                                                                                                <  3:93285/149‑1 (MQ=255)
                                                ccgCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGc                                                                                                <  1:266107/149‑1 (MQ=255)
                                                    aTTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTa                                                                                            <  2:414814/149‑1 (MQ=255)
                                                    aTTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTATTAATGCTTTa                                                                                            >  2:350434/1‑149 (MQ=255)
                                                     ttAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTAt                                                                                           >  4:668812/1‑149 (MQ=255)
                                                      tAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTTACTCCGTTGACTGTTTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTAtt                                                                                          <  2:253617/149‑1 (MQ=255)
                                                       aaGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTAttt                                                                                         <  1:12420/149‑1 (MQ=255)
                                                          cGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGcc                                                                                      <  1:457507/149‑1 (MQ=255)
                                                                cTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTtct                                                                                >  2:163889/1‑149 (MQ=255)
                                                                 ttCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTtctc                                                                               >  3:441500/1‑149 (MQ=255)
                                                                 ttCGTTTGCTGCGTTATTTACGCTGGCACTACACAGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTtctc                                                                               >  2:496596/1‑149 (MQ=255)
                                                                    gTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTg                                                                            >  2:627440/1‑149 (MQ=255)
                                                                      ttGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGaa                                                                          <  1:37210/149‑1 (MQ=255)
                                                                         ctgcGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTa                                                                       <  3:572733/149‑1 (MQ=255)
                                                                         ctgcGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTa                                                                       <  1:568668/149‑1 (MQ=255)
                                                                         ctgcGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTa                                                                       <  3:259666/149‑1 (MQ=255)
                                                                                    aCGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAACTTTTTAATGCTTTACTTCCCATTTCTCCTGAATTACGAAAACAttt                                                            <  4:465486/149‑1 (MQ=255)
                                                                                    aCGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACAttt                                                            <  3:5752/149‑1 (MQ=255)
                                                                                        tGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCaa                                                        <  1:627440/149‑1 (MQ=255)
                                                                                            aCTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACt                                                    <  1:542417/149‑1 (MQ=255)
                                                                                              tACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCg                                                  >  4:679172/1‑149 (MQ=255)
                                                                                              tACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCg                                                  <  2:189951/149‑1 (MQ=255)
                                                                                              tACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCg                                                  >  4:551062/1‑149 (MQ=255)
                                                                                                cacCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACa                                                               <  3:189986/134‑1 (MQ=255)
                                                                                                cacCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACa                                                               >  4:189986/1‑134 (MQ=255)
                                                                                                          cTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATaa                                      >  1:129494/1‑149 (MQ=255)
                                                                                                           tATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGCATTACGAAAACATTTGCAACACTCGATGTACCCATAAc                                     >  1:406292/1‑149 (MQ=255)
                                                                                                           tATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAAc                                     >  1:141552/1‑149 (MQ=255)
                                                                                                           tATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAAc                                     >  2:37400/1‑149 (MQ=255)
                                                                                                            aTGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACg                                    >  2:106654/1‑149 (MQ=255)
                                                                                                             tGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGa                                   <  3:679172/149‑1 (MQ=255)
                                                                                                             tGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGa                                   <  1:285518/149‑1 (MQ=255)
                                                                                                               aTTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATa                                 >  3:247569/1‑149 (MQ=255)
                                                                                                                 tttCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAAc                               <  4:213766/149‑1 (MQ=255)
                                                                                                                 tttCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAAc                               <  4:441500/149‑1 (MQ=255)
                                                                                                                   tCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCg                             <  4:83679/149‑1 (MQ=255)
                                                                                                                     ctctctTTGTAAAGCGACGGAGTGGTCACTCCGTGGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTAtt                                                                                          >  3:506223/1‑86 (MQ=255)
                                                                                                                     ctctctTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTAtt                                                                                          <  4:506223/86‑1 (MQ=255)
                                                                                                                     ctctctTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGt                           <  2:406292/149‑1 (MQ=255)
                                                                                                                     ctctctTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGt                           >  1:426769/1‑149 (MQ=255)
                                                                                                                      tctctTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTa                          >  1:567902/1‑149 (MQ=255)
                                                                                                                      tctctTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTa                          >  1:244195/1‑149 (MQ=255)
                                                                                                                      tctctTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTa                          >  1:441018/1‑149 (MQ=255)
                                                                                                                      tctctTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTAAGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTa                          >  4:499188/1‑149 (MQ=255)
                                                                                                                       ctctTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTaa                         >  2:226643/1‑149 (MQ=255)
                                                                                                                          tttGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAacac                      >  4:560838/1‑149 (MQ=255)
                                                                                                                          tttGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAacac                      >  2:442775/1‑149 (MQ=255)
                                                                                                                            agacaggCGACGGAGTGACCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCg                    >  2:332905/7‑149 (MQ=255)
                                                                                                                                aggCGACGGAGTGACCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGc                                                                                       >  2:397290/3‑78 (MQ=38)
                                                                                                                                aggCGACGGAGTGACCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGc                                                                                       <  1:397290/76‑1 (MQ=38)
                                                                                                                                aggCGACGGAGTGACCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGaaa                >  2:526893/3‑149 (MQ=255)
                                                                                                                                aaGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGaaa                <  2:98239/149‑1 (MQ=255)
                                                                                                                                aaGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGaaa                <  3:668812/149‑1 (MQ=255)
                                                                                                                                 aGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAg               >  2:170994/1‑149 (MQ=255)
                                                                                                                                 aGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCAATGTACCCATAACGATAACCGGTAACACCGGAAAg               <  4:603390/149‑1 (MQ=255)
                                                                                                                                  gCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGc              <  4:720801/149‑1 (MQ=255)
                                                                                                                                  gCGACGGAGTGACCACTCCGTCGCTTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGc              >  1:414594/1‑149 (MQ=255)
                                                                                                                                    gACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCAt            <  2:146070/149‑1 (MQ=255)
                                                                                                                                    gACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCAt            <  1:562551/149‑1 (MQ=255)
                                                                                                                                         aGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCaaa       <  2:129494/149‑1 (MQ=255)
                                                                                                                                          gTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAaac      <  3:131099/149‑1 (MQ=255)
                                                                                                                                          gTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAaac      <  3:719606/149‑1 (MQ=255)
                                                                                                                                           tGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAaaca     >  1:60364/1‑149 (MQ=255)
                                                                                                                                             gTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAaacaca   >  2:419321/1‑149 (MQ=255)
                                                                                                                                             gTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAaacaca   >  1:54045/1‑149 (MQ=255)
                                                                                                                                              tCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAaacacaa  >  4:213596/1‑149 (MQ=255)
                                                                                                                                              tCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAaacacaa  <  3:551062/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GCAAACGCTGGGGCTGGAGTGGCTCTACCGCCTGCTTTCGCAGCCGAGCCGCATTAAGCGTCAGCTTCGTTTGCTGCGTTATTTACGCTGGCACTACACCGGCAACCTATGATTTTCCTCTCTTTGTAAAGCGACGGAGTGGTCACTCCGTCGCCTGTCTGCTTTTTTATTACACAAAGCATTCAAATTTTTAATGCTTTATTTGCCATTTCTCCTGAATTACGAAAACATTTGCAACACTCGATGTACCCATAACGATAACCGGTAACACCGGAAAGCATGCAAACACAA  >  NC_000913/3980585‑3980875

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: