Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,437,908 T→A 7.0% intergenic (+14/‑116) trkA → / → mscL NAD‑binding component of TrK potassium transporter/mechanosensitive channel protein, high conductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,437,9080TA 7.0% 240.2 / 6.4 71intergenic (+14/‑116)trkA/mscLNAD‑binding component of TrK potassium transporter/mechanosensitive channel protein, high conductance
Reads supporting (aligned to +/- strand):  ref base T (31/35);  new base A (2/3);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

TTGCGCCGCGGCGTAGCAGAAGCTATTGAAGCCGTTGCTCACGGTGATGAAAGCACCTCACGCGTTGTCGGCAGAGTCATTGACGAAATCAAGCTACCGCCAGGAACGATTATTGGAGCGGTGGTACGTGGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGCGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTG  >  NC_000913/3437634‑3438158
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                          
ttGCGCCGCGGCGTAGCAGAAGCTATTGAAGCCGTTGCTCACGGTGATGAAAGCACCTCACGCGTTGTCGGCAGAGTCATTGACGAAATCAAGCTACCGCCAGGAACGATTATTGGAGCGGTGGTACGTGGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCATCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCATTTCTTCTTGTAATTAATAAGGAGTCt                                                                                                                                                                                                                                                            >  1:99372/1‑275 (MQ=255)
          gcgTAGCAGAAGTTATTGAGTCCGTTGCTCACGGGGAATAAAGCACCTCACGCGTTGTCGGCAGAGTCATTGACGAAATCAAGCGACCGCCAGGAACGATTATTGGAGCGGTGGTACGTGGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCt                                                                                                                                                                                                                                                  <  4:469047/275‑1 (MQ=255)
          gcgTAGCAGAAGCTATTGAAGCCGGTGCTCACGGTGATGAAAGGACCTCATGGGTAGTCGGCAGAGTCATTGACGAAATCAAGCTACCGACAGGAACGATTATTGGAGCGGTGGTACGTGGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTTAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCGTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCt                                                                                                                                                                                                                                                  <  4:160251/275‑1 (MQ=255)
                                               tgaAAGCACCTCACGCGTTGTCGGCAGAGTCATTGACGAAATCAAGCTACCGCCAGGAACGATTATTGGAGCGGTGGTACGTGGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAACGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGAAGCCTTATTATTTCCCt                                                                                                                                                                                                                                      >  4:239330/1‑250 (MQ=255)
                                                   aGCACCTCACGCGTTGTCGGCAGAGTCATTGACGAAATCAAGCTACCGCCAGGAACGATTATTGGAGCGGTGGTACGTGGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGa                                                                                                                                                                                                                                  <  2:239334/250‑1 (MQ=255)
                                                      aCCTCACGCGTTGTCGGCAGAGTCATTGACGAAATCAAGTTACCGCCAGGAACGATTATTGGAGCGGTGGTACGTGGAAATGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATACAAAGTTTATTACCGACGTCTAAAGACTCATCTAGCCAACCCCTTTCTTCATGTAATTAATAAGGTCACTAATGACGCCTTATTATTTCCCTTTGCTTATCCACGATTCATTAAattcattc                                                                                                                                                                                                        >  4:369124/1‑273 (MQ=255)
                                                                  gTCGGCAGAGTCATTGACGAAATCAAGCTACCGCCAGGAACGATTATTGGAGCGGTGGTACGTGGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCATGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATg                                                                                                                                                                                          >  2:45987/1‑275 (MQ=255)
                                                                                                 ggctaGGAAGGCTTATTTGTGCGGTGGTAAGTGGAACCGACGTGCTGCTTGCCAATGACAATGTGCGCATTTAGCAATGCGATCACGTAATTATGTTCCTGACAGATAAACAGTTTATTACCGCCGTCGAAAGACTCTTGCAGCCGAACCCTGTCTTCTTGTTAGTAATAATGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAACTTCATTCCTGGCAGGAAAAAGGCTTAACATTTGTTAGACTTTTTGTGGTCg                                                                                                                                                            <  4:304145/270‑1 (MQ=255)
                                                                                                     aGGAACGATTATTGGAGCGGTGGTACGTGGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCATGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTc                                                                                                                                                       <  2:369129/275‑1 (MQ=255)
                                                                                                             ttattGGAGCGGTGGTACGTGGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGgag                                                                                                                                               <  4:217095/275‑1 (MQ=255)
                                                                                                                  ggAGCGGTGGTACGTGGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAgg                                                                                                                                                                                                >  4:209990/1‑221 (MQ=255)
                                                                                                                      cGGTGGTACGTGGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGaaaa                                                                                                                                                                                            <  2:209993/221‑1 (MQ=255)
                                                                                                                        ggtgtACGGTGAAAAGACGTGATGATTGCCAATTGCAAGCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTGCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTAATTACCTTTGATTAGCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAg                                                                                                                                    <  4:213763/272‑1 (MQ=255)
                                                                                                                          tgtACGGTGAAACGACGTGATGATTGCCAACGACAATATGCGCATTGGGCAAGGCGATAACGTAATTATGTTCCTCAAAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTGTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTTCCTTTGATTATCAAGGATTAATTAAATTCATTCCTTGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGTAGAATAACATGAGCa                                                                                                                                  <  4:45987/274‑1 (MQ=255)
                                                                                                                                     acgacgTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCTGCTTCATAGGGAGAATAACATGAGAATAATTAAAGaa                                                                                                                       >  1:100257/1‑275 (MQ=255)
                                                                                                                                        acgTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAATAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCt                                                                                                                                                                                                                                                  <  8:269083/149‑1 (MQ=255)
                                                                                                                                         cgTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCtt                                                                                                                                                                                                                                                 <  6:394962/149‑1 (MQ=255)
                                                                                                                                             atgatTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCCACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCATTAGACGCCttatta                                                                                                                                                                                                                                             <  7:185649/149‑1 (MQ=255)
                                                                                                                                                atTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTAttt                                                                                                                                                                                                                                          <  7:130125/149‑1 (MQ=255)
                                                                                                                                                        tGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACCTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCtttta                                                                                                                                                                                                                                  <  8:342562/149‑3 (MQ=255)
                                                                                                                                                                  cgcATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAgg                                                                                                                                                  >  2:210096/1‑219 (MQ=255)
                                                                                                                                                                   gcATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCATTAGACGCCtta                                                                                                                                                                                                                                                >  8:425294/1‑124 (MQ=255)
                                                                                                                                                                   gcATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCATTAGACGCCtta                                                                                                                                                                                                                                                <  7:425294/124‑1 (MQ=255)
                                                                                                                                                                      ttGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGgaga                                                                                                                                              <  4:210093/219‑1 (MQ=255)
                                                                                                                                                                        gAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTGATGACGCCTTATTATTTCCCTTTGATTATCAAGGattaatt                                                                                                                                                                                                                  >  5:421298/1‑149 (MQ=255)
                                                                                                                                                                                 gATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAattcatt                                                                                                                                                                                                         >  6:243409/1‑149 (MQ=255)
                                                                                                                                                                                    cACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCt                                                                                                                                                                                                      <  8:307325/149‑1 (MQ=255)
                                                                                                                                                                                                   cTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGattaatta                                                                                                                                                                                                                 >  5:147019/1‑123 (MQ=255)
                                                                                                                                                                                                   cTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGattaatta                                                                                                                                                                                                                 <  6:147019/123‑1 (MQ=255)
                                                                                                                                                                                                   cTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACAtt                                                                                                                                                                                >  4:159439/1‑156 (MQ=255)
                                                                                                                                                                                                   cTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGa                                                                                                                                                                         >  2:205885/1‑163 (MQ=255)
                                                                                                                                                                                                   cTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGCGAGGGAACGTGGTGGATTTGGCGGTTGGTGTCATTATCGGTGc                                                         >  2:286112/1‑275 (MQ=255)
                                                                                                                                                                                                      acaGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTaa                                                                                                                                                                                    >  6:315735/1‑149 (MQ=255)
                                                                                                                                                                                                       caGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGtt                                                                                                                                                                            <  2:159441/156‑1 (MQ=255)
                                                                                                                                                                                                       caGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTa                                                                                                                                                                     <  4:205882/163‑1 (MQ=255)
                                                                                                                                                                                                              aaaGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAAt                                                                                                                                                                                           >  6:339621/1‑134 (MQ=255)
                                                                                                                                                                                                              aaaGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAAt                                                                                                                                                                                           <  5:339621/134‑1 (MQ=255)
                                                                                                                                                                                                               aaGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTGATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTa                                                                                                                                                                           <  6:421298/149‑1 (MQ=255)
                                                                                                                                                                                                               aaGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTa                                                                                                                                                                           <  7:308443/149‑1 (MQ=255)
                                                                                                                                                                                                                          cGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTa                                                                                                                                                                           <  7:239275/138‑1 (MQ=255)
                                                                                                                                                                                                                          cGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTa                                                                                                                                                                           >  8:239275/1‑138 (MQ=255)
                                                                                                                                                                                                                            aCGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCtt                                                                                                                                                        >  2:503536/1‑155 (MQ=255)
                                                                                                                                                                                                                               tCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGc                                                                                                                                                                                                                                                    >  8:440040/1‑60 (MQ=255)
                                                                                                                                                                                                                               tCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGc                                                                                                                                                                                                                                                    <  7:440040/60‑1 (MQ=255)
                                                                                                                                                                                                                                cGNAAGANTCTTCCAGCCAAGCCCTTTCTTCNTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATa                                                                                                                                                    <  4:503531/155‑1 (MQ=255)
                                                                                                                                                                                                                                   aaGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTc                                                                                                                                                       <  5:243409/149‑1 (MQ=255)
                                                                                                                                                                                                                                        tctTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAgg                                                                                                                                                  >  5:204483/1‑149 (MQ=255)
                                                                                                                                                                                                                                                  gatcttttttttttTGTAATTAATACGTCNTTTAATTCCGCTTTATTATTTTCATNTTATTATGAAGGATTAATTCAATTCCTTCCTTGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTTTCGGCTTCATNGGGCTAGTGACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGCGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGc          <  3:99372/269‑1 (MQ=255)
                                                                                                                                                                                                                                                       cTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGCGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGcctcct     <  2:323870/275‑1 (MQ=255)
                                                                                                                                                                                                                                                          tcttctTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCatt                                                                                                                                <  5:315735/149‑1 (MQ=255)
                                                                                                                                                                                                                                                          tcttctTGTAATTAATAAGGCGTCTAATGACCCCTTATTATTTCTCTTTGATTATCAAGGATTAATTAAATTCATTCCGGGCAGCAAAATGGCTTAACATCTGTTAGACTAATGGTTTTCGGCTGCATAGGGAGCATACCATGACCATTATTAAAGAATTTCGCGAATTTGCGATGCTCAGGAACGTGGTGGCTTTGGCGGTGGGTTTCATTCCCGGTCCGTCATTCTGGAAGTTTGCATTTTCCCTGGTTGCCGTTATCACCTTGCCATctctt  >  4:319877/1‑274 (MQ=255)
                                                                                                                                                                                                                                                              cttGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCTCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCAttatta                                                                                                                            >  8:276684/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              cttGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACAt                                                                                                                                                                                 >  2:408088/1‑96 (MQ=255)
                                                                                                                                                                                                                                                               ttGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAat                                                                                                                                                                                                              >  5:101675/1‑66 (MQ=255)
                                                                                                                                                                                                                                                               ttGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAat                                                                                                                                                                                                              <  6:101675/66‑1 (MQ=255)
                                                                                                                                                                                                                                                                tGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTGAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTaaa                                                                                                                          <  7:276684/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                  taattaatAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGt                                                                                                                                                                             <  4:408083/96‑1 (MQ=255)
                                                                                                                                                                                                                                                                  taattaatAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGa                                                                                                                        >  6:254594/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    attaatAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCgg                                                                                                                                                           >  1:54167/1‑112 (MQ=255)
                                                                                                                                                                                                                                                                      taataaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCAt                                                                                                                                                     >  2:516713/1‑116 (MQ=255)
                                                                                                                                                                                                                                                                       aataaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTc                                                                                                                   <  6:204483/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ataaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATg                                                                                                                                                                   >  6:324393/1‑100 (MQ=255)
                                                                                                                                                                                                                                                                        ataaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATg                                                                                                                                                                   <  5:324393/100‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ataaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTc                                                                                                                                                       <  3:54167/112‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ataaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAAGATGAGCATTATTAGAGAATTTcg                                                                                                                  >  7:298094/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         taaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTcgc                                                                                                                 >  7:337810/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         taaGGCGTCATTAGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTc                                                                                                                   >  7:9986/1‑147 (MQ=255)
                                                                                                                                                                                                                                                                         taaGGCGTCATTAGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTc                                                                                                                   <  8:9986/147‑1 (MQ=255)
                                                                                                                                                                                                                                                                          aaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGATTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGCATTTcgcg                                                                                                                >  5:11898/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                          aaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAggg                                                                                                                                                 <  4:516708/116‑1 (MQ=255)
                                                                                                                                                                                                                                                                               gTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCGTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAAttt                                                                                                           >  6:375165/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                                tCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCa                                                                                                                                                      >  2:458373/1‑105 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                          
TTGCGCCGCGGCGTAGCAGAAGCTATTGAAGCCGTTGCTCACGGTGATGAAAGCACCTCACGCGTTGTCGGCAGAGTCATTGACGAAATCAAGCTACCGCCAGGAACGATTATTGGAGCGGTGGTACGTGGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGCGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTG  >  NC_000913/3437634‑3438158

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