Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,981,102 T→A 6.7% V72V (GTT→GTA yifK → putative APC family amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,981,1020TA 6.7% 469.2 / 20.0 135V72V (GTT→GTAyifKputative APC family amino acid transporter
Reads supporting (aligned to +/- strand):  ref base T (64/62);  new base A (4/5);  total (68/67)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.89e-04

Only 100 of 145 total aligned reads displayed.

ACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTATGTCCAGTTCTGGTTCCCGGAGATGGCGCAGTGGATACCCGCATTGATCGCAGTGGCGCTGGTGGCGTTGGCGAATCTGGCGGCGGTGCGGTTGTACGGCGAAATCGAGTTCTGGTTCGCGATGATCAAAGTCAC  >  NC_000913/3980839‑3981362
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                    
    tAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAACTGTTGTTCCTCGAACCGGTTGCCGGTTCGTTCGcc                                                                                                                                                                                                                                                       >  2:150631/1‑275 (MQ=255)
    tAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGcc                                                                                                                                                                                                                                                       >  2:40932/1‑275 (MQ=255)
                              gtataagagacaGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGcc                                                                                                                                                                                                                              <  1:150848/263‑1 (MQ=255)
                                                      gATAAAAAACCAGGGCTACAGCGGGGGCTGGAAGCTCGACATATCGAACTCATGGCCCTGGGGGGAACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCt                                                                                                                                                                                                      <  4:40931/274‑1 (MQ=255)
                                                            aaaCCGGGGCTACGGCGGGGTCTGGAAGCNNNCCGTATCGAACTCACCGCCCTGGGGGGCACCATTGGGGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTa                                                                                                                                                                                               <  3:151264/275‑1 (MQ=255)
                                                                   agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATCGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCCTCCTGAGTTCAATTGGCGACATGTTTTTCCTCGACCCGGTTACCGGTTCGTTCGTCGTTTATGCGCATCGTTATTTGAGCCCGATCTTTTGCTTTCCCCCCGCCCTGCCTTACTTGttt                                                                                                                                                                                        >  4:214393/1‑275 (MQ=255)
                                                                             tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGTCGTGTATGTGCCTCGTTTTATGAGCCCGTTCTTTGGCTATCTTACCACCTGGTCTTACTGGgttgggtctgtc                                                                                                                                                                               >  1:18411/1‑265 (MQ=255)
                                                                             tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTAGGCGCATCGGTATATGAGCCCGTTCTTTGGCTTTCCCCCCGCCTGGTCTTACCGGttttttgtggtgg                                                                                                                                                                              >  4:469923/1‑265 (MQ=255)
                                                                                    gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTggg                                                                                                                                                                       <  2:170318/275‑1 (MQ=255)
                                                                                     aaGCTCGACATTTCAAACTCATGGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCCTCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTggg                                                                                                                                                                       <  2:469924/274‑1 (MQ=255)
                                                                                                    ccgtatgcGCCGTGGGGGGACCTTTGGCGGAGGCCGTTTTGTGGGGGCCGCCAGTACCCTCAGATGGGCCGTGTCTTCCGTTTTGTTGGCCTATATCGTCGCNGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGttt                                                                                                                                                                                        <  3:111935/234‑1 (MQ=255)
                                                                                                                      tacccttgggCGTCGCCTTGTGTACGGGTGCGCCCAGTACCCGCAAGTGGGCCGGGCCATCCGTATGGTTGGCCTATATCTTCGCCGGGCTTTTGGTCGTTTTCATCATGCGTTCAATGGGCGAAATGTGGTTCCTCGAACCGGTTACCGGTTCGTTCTCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTATGTc                                                                                                                                     <  4:111133/267‑1 (MQ=255)
                                                                                                                      gCACCATGGGCGTCGGCCTGTTTAGGGGGGCCTCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGACCCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTATGTc                                                                                                                                     <  1:41690/275‑1 (MQ=255)
                                                                                                                       cACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggtaacc                                                                                                                                                                                                                                                                   >  5:777305/1‑148 (MQ=255)
                                                                                                                            ttGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTATGTCCAGTTc                                                                                                                               >  2:233314/1‑275 (MQ=255)
                                                                                                                             tGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTAACCGgttcg                                                                                                                                                                                                                                                             <  6:777305/148‑1 (MQ=255)
                                                                                                                                              tGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCa                                                                                                                                                                                                                                            >  8:824341/1‑148 (MQ=255)
                                                                                                                                               gggggCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCCCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTATTTCCCGTTCTGGTTCCCGGAGATGgcc                                                                                                             >  4:143072/1‑273 (MQ=255)
                                                                                                                                               gggggCCGCCAGTAACCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTGCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTATGTCCAGTTCTGGTTCCCGGAGATAgcg                                                                                                             <  2:407737/274‑1 (MQ=255)
                                                                                                                                                ggggCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGt                                                                                                                                                                                                                                                     >  8:1005753/1‑137 (MQ=255)
                                                                                                                                                ggggCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGt                                                                                                                                                                                                                                                     <  7:1005753/137‑1 (MQ=255)
                                                                                                                                                  ggCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTATGTCCAGTTCTGGTTCCCGGAGATGGCGCAg                                                                                                          <  2:143073/274‑1 (MQ=255)
                                                                                                                                                   gccgccAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTGTATGCGCATCGtt                                                                                                                                                                                                                                       <  5:844720/148‑1 (MQ=255)
                                                                                                                                                    ccgccAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGNGCATCGTTATATTAGCCCCGTCTTTTGCTCTCTCCCCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGCTCTCTGNAACCACCGACATTGGCGGTTTTGTGNAGTGTGGGTTTCCGGAGTTgggcgcgggg                                                                                                       >  3:84754/1‑267 (MQ=255)
                                                                                                                                                    ccgccAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATg                                                                                                                                                                                                                                  >  4:315066/1‑152 (MQ=255)
                                                                                                                                                    ccgccAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCa                                                                                                                                                   >  4:376304/1‑231 (MQ=255)
                                                                                                                                                    ccgccAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTATGTCCAGTTCTGGTTCCCGGAGATGGCGCAGTg                                                                                                        >  2:104112/1‑274 (MQ=255)
                                                                                                                                                     cgccAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATgcgc                                                                                                                                                                                                                                             >  6:231939/1‑140 (MQ=255)
                                                                                                                                                     cgccAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATgcgc                                                                                                                                                                                                                                             <  5:231939/140‑1 (MQ=255)
                                                                                                                                                        cAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCTGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTgg                                                                                                                                               <  2:376305/231‑1 (MQ=255)
                                                                                                                                                        cAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGcc                                                                                                                                                                                                                              <  2:315067/152‑1 (MQ=255)
                                                                                                                                                               cTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGt                                                                                                                                                                                                                                        <  6:110945/135‑1 (MQ=255)
                                                                                                                                                                tGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAg                                                                                                                                                                                                                                <  5:1114791/142‑1 (MQ=255)
                                                                                                                                                                 gAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccggtt                                                                                                                                                                                                                                                               <  8:847516/110‑1 (MQ=255)
                                                                                                                                                                 gAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccggtt                                                                                                                                                                                                                                                               >  7:847516/1‑110 (MQ=255)
                                                                                                                                                                 gAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTc                                                                                                                                                                                                                         >  6:804158/1‑148 (MQ=255)
                                                                                                                                                                    aTGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTtata                                                                                                                                                                                                                                    <  7:654069/134‑1 (MQ=255)
                                                                                                                                                                    aTGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTtata                                                                                                                                                                                                                                    >  8:654069/1‑134 (MQ=255)
                                                                                                                                                                       ggccgggccATCCGTATTGTTGGCCTATATCATCGCGGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAGATGTTGTTCCTAGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCCCACAGCCTGGCCTTCCTGGTTTATGTGGCTGGCGGTGGGGATCTCTGAAATCACCGGCATTGGCGGTTATGTCCAGTCTTGGTTCACGGAGATGGCGCAGTGGCTACCCACATTCATCGCAg                                                                                    >  4:176636/1‑275 (MQ=255)
                                                                                                                                                                                 tCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTAACCGgttcg                                                                                                                                                                                                                                                             <  7:1028202/96‑1 (MQ=255)
                                                                                                                                                                                 tCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTAACCGgttcg                                                                                                                                                                                                                                                             >  5:788223/1‑96 (MQ=255)
                                                                                                                                                                                            ggCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGgttcgtt                                                                                                                                                                                                                                                           >  8:351420/1‑87 (MQ=255)
                                                                                                                                                                                              ccTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccggtt                                                                                                                                                                                                                                                               <  8:132214/81‑1 (MQ=255)
                                                                                                                                                                                              ccTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggttaccggtt                                                                                                                                                                                                                                                               >  7:132214/1‑81 (MQ=255)
                                                                                                                                                                                              ccTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTg                                                                                                                                                                                            <  8:695232/148‑1 (MQ=255)
                                                                                                                                                                                                tataTCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGtt                                                                                                                                                                                                                          <  6:345064/116‑1 (MQ=255)
                                                                                                                                                                                                  tatCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGgttc                                                                                                                                                                                                                                                              >  6:779211/1‑78 (MQ=255)
                                                                                                                                                                                                  tatCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCgg                                                                                                                                                                           >  4:249331/1‑161 (MQ=255)
                                                                                                                                                                                                   atcatcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGtt                                                                                                                                                                                                                                                    >  5:818938/1‑87 (MQ=255)
                                                                                                                                                                                                   atcatcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGtt                                                                                                                                                                                                                                                    <  6:818938/87‑1 (MQ=255)
                                                                                                                                                                                                    tcatcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAaccggtta                                                                                                                                                                                                                                                                     >  8:573410/1‑69 (MQ=255)
                                                                                                                                                                                                     catcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCt                                                                                                                                                                                                 >  7:950741/1‑136 (MQ=255)
                                                                                                                                                                                                      atcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTtata                                                                                                                                                                                                                                    >  8:127178/1‑100 (MQ=255)
                                                                                                                                                                                                      atcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTggg                                                                                                                                                                       <  2:249332/161‑1 (MQ=255)
                                                                                                                                                                                                       tcGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTgg                                                                                                                                                                                           >  7:719650/1‑140 (MQ=255)
                                                                                                                                                                                                        cGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTATTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTATGTCCAGTTCTGGTTCCCGGAGATGGCGCAGTGGATACCCGCATTGATCGCAGTGGAGCTGGTGGCGTTGGCGAATCTGGcggcgg                                                   >  2:348204/1‑275 (MQ=255)
                                                                                                                                                                                                          ccGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTgg                                                                                                                                                                                                                    >  5:1006837/1‑112 (MQ=255)
                                                                                                                                                                                                          ccGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTgg                                                                                                                                                                                                                    >  7:677743/1‑112 (MQ=255)
                                                                                                                                                                                                            gggCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTATGTCCAGTTCTGGTTCCCGGAGATGGCGCAGTGGATACCCGCATTGATCGCAGTGGCGCTGGTGGCGTTGGCGAATCTGGCGgcggtccg                                               >  2:497291/1‑275 (MQ=255)
                                                                                                                                                                                                                   tCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTtat                                                                                                                                                                                                                                     <  6:265353/86‑1 (MQ=255)
                                                                                                                                                                                                                   tCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTtat                                                                                                                                                                                                                                     >  5:265353/1‑86 (MQ=255)
                                                                                                                                                                                                                   tCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCttt                                                                                                                                                                                                                      >  7:30797/1‑101 (MQ=255)
                                                                                                                                                                                                                   tCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCttt                                                                                                                                                                                                                      <  8:30797/101‑1 (MQ=255)
                                                                                                                                                                                                                     gTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTa                                                                                                                                                                                       <  6:212226/130‑1 (MQ=255)
                                                                                                                                                                                                                       cTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGtt                                                                                                                                                                                                                                                    <  6:791856/67‑1 (MQ=255)
                                                                                                                                                                                                                       cTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGtt                                                                                                                                                                                                                                                    >  5:791856/1‑67 (MQ=255)
                                                                                                                                                                                                                       cTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGt                                                                                                                                                                                                                                        >  8:947255/1‑79 (MQ=255)
                                                                                                                                                                                                                       cTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGt                                                                                                                                                                                                                                        <  7:947255/79‑1 (MQ=255)
                                                                                                                                                                                                                         ttttCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGtt                                                                                                                                                                                                                                                    <  7:328351/65‑1 (MQ=255)
                                                                                                                                                                                                                         ttttCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGtt                                                                                                                                                                                         >  8:913119/1‑124 (MQ=255)
                                                                                                                                                                                                                         ttttCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTAtg                                                                                                                                                                                     <  8:1052551/128‑1 (MQ=255)
                                                                                                                                                                                                                         ttttCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGt                                                                                                                                                                          >  7:124596/1‑139 (MQ=255)
                                                                                                                                                                                                                         ttttCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGt                                                                                                                                                                          <  8:124596/139‑1 (MQ=255)
                                                                                                                                                                                                                             catcatGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTAc                                                                                                                                                                                              <  7:891712/115‑1 (MQ=255)
                                                                                                                                                                                                                                  tGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTgg                                                                                                                                                                                           <  5:65063/113‑1 (MQ=255)
                                                                                                                                                                                                                                  tGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTgg                                                                                                                                                                                           >  6:65063/1‑113 (MQ=255)
                                                                                                                                                                                                                                          aTGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCtt                                                                                                                                                                                                >  6:602575/1‑100 (MQ=255)
                                                                                                                                                                                                                                           tGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGt                                                                                                                                                                          >  7:125707/1‑121 (MQ=255)
                                                                                                                                                                                                                                           tGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGt                                                                                                                                                                          >  7:246990/1‑121 (MQ=255)
                                                                                                                                                                                                                                            gggCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTATGTCCAGTTCTGGTTCCCGGAGATGGCGCAGTGGATACCCGCATTGATCGCAGTGGCGCTGGTGGCGTTGGCGAATCTGGCGGCGGTGCGGTTGTACGGCGAAATCGAGTTCTGGTTCGCga               <  2:214394/275‑1 (MQ=255)
                                                                                                                                                                                                                                               cGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGt                                                                                                                                                                                          <  6:1040234/101‑1 (MQ=255)
                                                                                                                                                                                                                                               cGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGt                                                                                                                                                                                          >  5:1040234/1‑101 (MQ=255)
                                                                                                                                                                                                                                                 aaaTGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTg                                                                                                                                                                                                     >  6:248445/1‑88 (MQ=255)
                                                                                                                                                                                                                                                     gttgttCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGcc                                                                                                                                                                                                       >  5:243396/1‑82 (MQ=255)
                                                                                                                                                                                                                                                         ttCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCtt                                                                                                                                                                                                                       <  7:458655/62‑1 (MQ=255)
                                                                                                                                                                                                                                                         ttCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTATGTCCAGTTCTGGTTCCCGGAGATGGCGCAGTGGATACCCGCATTGATCGCAGTGGCGCTGGTGGCGTTGGCGAATCTGGCGGCGGTGCGGTTGTACGGCGAAATCGAGTTCTGGTTCGCGATGATCAAAGTcac  <  1:84754/275‑1 (MQ=255)
                                                                                                                                                                                                                                                            aacGAACCGGTAACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCttt                                                                                                                                                                                                                      <  7:341539/58‑1 (MQ=38)
                                                                                                                                                                                                                                                            aacGAACCGGTAACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCttt                                                                                                                                                                                                                      >  8:341539/3‑60 (MQ=38)
                                                                                                                                                                                                                                                             tCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGtt                                                                                                                                                                                         >  5:292930/1‑88 (MQ=255)
                                                                                                                                                                                                                                                             tCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATctct                                                                                                                                                                <  5:690398/113‑1 (MQ=255)
                                                                                                                                                                                                                                                             tCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCa                                                                                                                                                         <  6:175655/120‑1 (MQ=255)
                                                                                                                                                                                                                                                             tCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCa                                                                                                                                                         >  5:175655/1‑120 (MQ=255)
                                                                                                                                                                                                                                                              cGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTta                                                                                                                                                                                                                                      >  5:598074/1‑42 (MQ=255)
                                                                                                                                                                                                                                                              cGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTta                                                                                                                                                                                                                                      <  8:713576/42‑1 (MQ=255)
                                                                                                                                                                                                                                                              cGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTta                                                                                                                                                                                                                                      >  7:713576/1‑42 (MQ=255)
                                                                                                                                                                                                                                                               gAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTATg                                                                                                                                       <  8:252912/136‑1 (MQ=255)
                                                                                                                                                                                                                                                               gAACCGGTAACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTgg                                                                                                                                                                        >  8:401218/1‑103 (MQ=255)
                                                                                                                                                                                                                                                               gAACCGGTAACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTgg                                                                                                                                                                        <  7:401218/103‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ttaccggttCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCt                                                                                                                                                                                                      >  5:474778/1‑66 (MQ=255)
                                                                                                                                                                                                                                                                      ttaccggttCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCt                                                                                                                                                                                                      <  6:474778/66‑1 (MQ=255)
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                    
ACGATAACCGGTAACACCGGAAAGCATGCAAACACAACACGAGGATTTATGGCAGATAACAAACCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCTGAAATGGGCCGGGCCATCCGTATTGTTGGCCTATATCATCGCCGGGCTGTTCGTCTTTTTCATCATGCGTTCAATGGGCGAAATGTTGTTCCTCGAACCGGTTACCGGTTCGTTCGCCGTTTATGCGCATCGTTATATGAGCCCGTTCTTTGGCTATCTCACCGCCTGGTCTTACTGGTTTATGTGGATGGCGGTGGGGATCTCTGAAATCACCGCCATTGGCGTTTATGTCCAGTTCTGGTTCCCGGAGATGGCGCAGTGGATACCCGCATTGATCGCAGTGGCGCTGGTGGCGTTGGCGAATCTGGCGGCGGTGCGGTTGTACGGCGAAATCGAGTTCTGGTTCGCGATGATCAAAGTCAC  >  NC_000913/3980839‑3981362

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