Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,437,910 A→T 7.0% intergenic (+16/‑114) 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,9100AT 7.0% 249.4 / 7.1 71intergenic (+16/‑114)trkA/mscLNAD‑binding component of TrK potassium transporter/mechanosensitive channel protein, high conductance
Reads supporting (aligned to +/- strand):  ref base A (31/35);  new base T (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 = 7.50e-01

GCGTAGCAGAAGCTATTGAAGCCGTTGCTCACGGTGATGAAAGCACCTCACGCGTTGTCGGCAGAGTCATTGACGAAATCAAGCTACCGCCAGGAACGATTATTGGAGCGGTGGTACGTGGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGCGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTG  >  NC_000913/3437644‑3438158
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        
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                                                                                                                                                                                                                 <  6:147019/123‑1 (MQ=255)
                                                                                                                                                                                         cTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGattaatta                                                                                                                                                                                                                 >  5:147019/1‑123 (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                                                                                                                                                                                                                                                    <  7:440040/60‑1 (MQ=255)
                                                                                                                                                                                                                     tCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGc                                                                                                                                                                                                                                                    >  8:440040/1‑60 (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                                                                                                                                                                                                              <  6:101675/66‑1 (MQ=255)
                                                                                                                                                                                                                                                     ttGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAat                                                                                                                                                                                                              >  5:101675/1‑66 (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)
                                                                                                                                                                                                                                                                          aTGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAgg                                                                                                                                                  <  4:458368/105‑1 (MQ=255)
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        
GCGTAGCAGAAGCTATTGAAGCCGTTGCTCACGGTGATGAAAGCACCTCACGCGTTGTCGGCAGAGTCATTGACGAAATCAAGCTACCGCCAGGAACGATTATTGGAGCGGTGGTACGTGGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGCGCGGGAACGTGGTGGATTTGGCGGTGGGTGTCATTATCGGTGCGGCATTCGGGAAGATTGTCTCTTCACTGGTTGCCGATATCATCATGCCTCCTCTG  >  NC_000913/3437644‑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: