Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,437,910 A→T 5.3% 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 5.3% 489.1 / 9.8 161intergenic (+16/‑114)trkA/mscLNAD‑binding component of TrK potassium transporter/mechanosensitive channel protein, high conductance
Reads supporting (aligned to +/- strand):  ref base A (79/73);  new base T (5/4);  total (84/77)
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.26e-01

Only 100 of 165 total aligned reads displayed.

GTGGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATG  >  NC_000913/3437761‑3438059
                                                                                                                                                     |                                                                                                                                                     
gTGGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTaa                                                                                                                                                       >  2:793053/1‑150 (MQ=255)
 tGGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAAt                                                                                                                                                      >  1:1479576/1‑150 (MQ=255)
  ggAAACGCCGTGAAGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATCTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTAGTAATTAATAAGGAGTCTAATg                                                                                                                                                     >  2:658694/1‑150 (MQ=255)
   gAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATgc                                                                                                                                                    >  2:505637/1‑149 (MQ=255)
    aaaCGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGAc                                                                                                                                                   <  1:889826/150‑1 (MQ=255)
         acgTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCtt                                                                                                                                              <  1:1006154/150‑1 (MQ=255)
           gttattatTGCCAATGACAATCTGCTCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCttat                                                                                                                                            <  1:292347/144‑1 (MQ=255)
           gTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGCCGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCttat                                                                                                                                            <  1:1082796/150‑1 (MQ=255)
           gTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCttat                                                                                                                                            >  2:873162/1‑150 (MQ=255)
           gTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCATTAGACGCCttat                                                                                                                                            <  1:247231/150‑1 (MQ=255)
            tgatgaTTGCCAATGACAATCTGCGGATTGAGCAAGGCGTTCACGTAATTATGTTCCTCTCAGATTAAATGTTTATTACCGACGTCGAAAGACTCTTCCAGGAAATCCATTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCttatt                                                                                                                                           <  1:1347829/150‑1 (MQ=255)
            tgatgaTTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCttatt                                                                                                                                           <  2:374473/150‑1 (MQ=255)
             gatgatTGCCAATGACAATCTGCGCATTGAGCAAGGCGATTACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCATTAGACGCCttatta                                                                                                                                          <  2:1633705/150‑1 (MQ=255)
                   tGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTccc                                                                                                                                    <  1:941062/150‑1 (MQ=255)
                       aaTGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTg                                                                                                                                <  2:132648/150‑1 (MQ=255)
                        aTGACAATCTGAGCATTGAGCAAGGCGATCACGTAATTATGTTCCTGACAGATAAAAAGGTTATTAGCGACGGCGAATGACTCTTCGAGCCAAGCCCTTTCTTCTCGTAATTAATAAGGCTTCTAATGACGCCTTATTATTTCCCTTTGa                                                                                                                               <  1:849008/150‑1 (MQ=255)
                           aCAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTa                                                                                                                            <  2:111066/150‑1 (MQ=255)
                              aTCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTATAATGACGCCGTATTATTTCACTTTGATTATCa                                                                                                                         >  1:773317/1‑150 (MQ=255)
                                cTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAg                                                                                                                       <  1:954626/150‑1 (MQ=255)
                                cTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAg                                                                                                                       <  2:208350/150‑1 (MQ=255)
                                cTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAg                                                                                                                       <  1:783408/150‑1 (MQ=255)
                                    gcATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGatt                                                                                                                   <  2:1516976/150‑1 (MQ=255)
                                     cATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGatta                                                                                                                  >  1:837580/1‑150 (MQ=255)
                                       ttGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGattaat                                                                                                                <  1:1576600/150‑1 (MQ=255)
                                         gAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGattaatta                                                                                                              >  2:42807/1‑150 (MQ=255)
                                         gAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGattaatta                                                                                                              <  1:1012794/150‑1 (MQ=255)
                                            cAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAat                                                                                                           <  1:373321/150‑1 (MQ=255)
                                            cAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACAGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAat                                                                                                           <  2:1279189/150‑1 (MQ=255)
                                             aaGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAatt                                                                                                          >  2:565013/1‑150 (MQ=255)
                                              aGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAattc                                                                                                         <  1:1343027/150‑1 (MQ=255)
                                                  gATAACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCATCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGTCGTCTAAAGACGCCTTATTATTTCCCTTTGATTATCAATGATTAATTAAattcatta                                                                                                     >  2:801094/1‑149 (MQ=255)
                                                     cACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAg                                                                                                                       >  2:1561584/1‑129 (MQ=255)
                                                                gTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATg                                                                                       <  1:1239880/150‑1 (MQ=255)
                                                                gTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATg                                                                                       <  2:606029/150‑1 (MQ=255)
                                                                  tCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCACTTTGATTATAAAGGATTAATTAAATTCATTACTGGCAGGAAAATGGc                                                                                     >  2:364275/1‑150 (MQ=255)
                                                                   ccTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCt                                                                                    >  2:255992/1‑150 (MQ=255)
                                                                    cTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCtt                                                                                   >  2:685967/1‑150 (MQ=255)
                                                                      cacaGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTaa                                                                                 >  2:31647/1‑150 (MQ=255)
                                                                      cacaGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCAAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTaa                                                                                 <  1:41298/150‑1 (MQ=255)
                                                                        caGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACa                                                                               >  1:749692/1‑150 (MQ=255)
                                                                        caGATAAAAAGTTTAATACCGACGTCGAAAGACTCTTCCATCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACa                                                                               <  2:219216/150‑1 (MQ=255)
                                                                          gATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACAtt                                                                             >  2:1003489/1‑150 (MQ=255)
                                                                            tAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTg                                                                           >  1:880481/1‑150 (MQ=255)
                                                                              aaaaGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGtt                                                                         <  2:855506/150‑1 (MQ=255)
                                                                               aaaGTTTATTTCCGACGTCGAAAGCCTCTTCCAGCCAAGCCCGTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTACCTTTGATTATCAAGGAGTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAGCATTTGttt                                                                        >  2:1231683/1‑149 (MQ=255)
                                                                               aaaGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTa                                                                        >  1:354604/1‑150 (MQ=255)
                                                                               aaaGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTa                                                                        <  1:549338/150‑1 (MQ=255)
                                                                               aaaGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTa                                                                        <  1:873771/150‑1 (MQ=255)
                                                                                  gTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTAATTCCTGGCAGGAAAATGGCTTAACATATGTTAGAc                                                                     >  1:1263886/1‑150 (MQ=255)
                                                                                    ttattaCCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACtt                                                                   >  2:1726681/1‑150 (MQ=255)
                                                                                     tattaCCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTa                                                                  >  1:935059/1‑150 (MQ=255)
                                                                                      attaCCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTAt                                                                 >  2:889689/1‑150 (MQ=255)
                                                                                         aacgtcGTCGAAATACTCTTCCTGACAAGACCTTTCTTCTTGGAATTAATAAGGCGTCTATTGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTTAATTCATTCCTGGCAGGAAATTGGCTTAACATTTGTTAGACTTATGGt                                                              <  1:1634649/145‑1 (MQ=255)
                                                                                          ccGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGtt                                                             >  2:1201861/1‑150 (MQ=255)
                                                                                               gTCGAAAGCCTCTTCCAGCCAAGCCCGTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAAGTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACGTATTGTTGTCgg                                                        >  2:939378/1‑150 (MQ=255)
                                                                                                tCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGTCTTAACATTTGTTAGACTTATGGTTGTCGGc                                                       >  1:812871/1‑150 (MQ=255)
                                                                                                 cGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCt                                                      >  2:1451506/1‑150 (MQ=255)
                                                                                                   aaaGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGAGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTc                                                    <  2:428781/150‑1 (MQ=255)
                                                                                                     aGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCAt                                                  <  2:1544740/150‑1 (MQ=255)
                                                                                                      gACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATa                                                 <  1:558629/150‑1 (MQ=255)
                                                                                                        ctctTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAgg                                               >  1:1617238/1‑150 (MQ=255)
                                                                                                        ctctTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAgg                                               <  1:1112997/150‑1 (MQ=255)
                                                                                                          ctTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGga                                             >  1:1519164/1‑150 (MQ=255)
                                                                                                           ttCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGgag                                            <  2:1173519/150‑1 (MQ=255)
                                                                                                              cAGCCAAGCCGTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAGATTCATTCCTGGCAGGAAAATGGATTAACATTTGTTAGACATATGGTTTTCGGCTTCATAGGGAGAAt                                         >  1:342503/1‑150 (MQ=255)
                                                                                                              cAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAAt                                         >  1:683673/1‑150 (MQ=255)
                                                                                                               aGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTAATTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATa                                        <  1:278786/150‑1 (MQ=255)
                                                                                                               aGCAAATCCCTTTCTTATTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATa                                        <  1:339703/150‑1 (MQ=255)
                                                                                                                  cAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACa                                     >  2:1283685/1‑150 (MQ=255)
                                                                                                                  cAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACa                                     >  2:1699437/1‑150 (MQ=255)
                                                                                                                     gCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGa                                  >  2:381628/1‑150 (MQ=255)
                                                                                                                      cccTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGACTTATTATTTCCCTTTGATTATCTAGGATTAATTAAATTAATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGat                                 >  1:89129/1‑149 (MQ=255)
                                                                                                                       ccTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGc                                >  2:10625/1‑150 (MQ=255)
                                                                                                                         tttCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCat                              <  2:1489606/150‑1 (MQ=255)
                                                                                                                            cttcttGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCattat                           >  1:549339/1‑150 (MQ=255)
                                                                                                                            cttcttGTAATTAATAAGGCGTCTAATGAAGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCATGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGTAGAATAACATGAGCattat                           >  1:1369605/1‑150 (MQ=255)
                                                                                                                             ttcttGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCATGAATATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCattatt                          <  1:1092544/150‑1 (MQ=255)
                                                                                                                              tcttGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATa                                        >  1:1643353/1‑135 (MQ=255)
                                                                                                                              tcttGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATa                                        <  2:1643353/135‑1 (MQ=255)
                                                                                                                              tcttGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCAttatta                         >  1:1236842/1‑150 (MQ=255)
                                                                                                                                 tGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAg                      <  1:987118/150‑1 (MQ=255)
                                                                                                                                  gTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGa                     >  1:322803/1‑150 (MQ=255)
                                                                                                                                     attaatAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAAtt                  >  2:633441/1‑150 (MQ=255)
                                                                                                                                     aataataaGGCGTCATTAGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAAtt                  >  2:1168447/3‑150 (MQ=255)
                                                                                                                                      ataataaGGCGTCATTAGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAAttt                 >  2:963939/2‑150 (MQ=255)
                                                                                                                                       taataaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTc                <  2:836674/150‑1 (MQ=255)
                                                                                                                                         ataaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTcgg              >  2:670849/1‑149 (MQ=255)
                                                                                                                                         ataaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTcgc              >  1:485716/1‑150 (MQ=255)
                                                                                                                                          taaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTcgcg             <  2:269885/150‑1 (MQ=255)
                                                                                                                                          taaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTcgcg             >  2:937136/1‑150 (MQ=255)
                                                                                                                                          taaGGAGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAATGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTcgcg             <  1:389216/150‑1 (MQ=255)
                                                                                                                                           aaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGa            >  2:1593572/1‑150 (MQ=255)
                                                                                                                                                gTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTg       >  2:1205538/1‑150 (MQ=255)
                                                                                                                                                gTCTAATGACGCCTTATTATTTACCTTTTATTATCAAGTATTTATTAAATTCATTCCTGGCAGGAAAATGGCTTAACAATTGTTAGACTTATGGTTGTCGGCCTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTg       <  1:1564240/150‑1 (MQ=255)
                                                                                                                                                 tCTAATTACGACTTATTATTTCCCTTTGATTATCTAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCTGCTTCATAGGGTGAATAACATGAGCTTTATTAAAGAATTTCGCGAATTTGc      <  1:801094/150‑1 (MQ=255)
                                                                                                                                                  cTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCg     <  2:195144/150‑1 (MQ=255)
                                                                                                                                                   tAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGa    <  1:649644/150‑1 (MQ=255)
                                                                                                                                                   tAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATCGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGa    <  1:316602/150‑1 (MQ=255)
                                                                                                                                                     aTGACTCCTTATTATTTGACTTTGATTATCAAGTATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATg  <  1:931107/150‑1 (MQ=255)
                                                                                                                                                     aTGACGCCTTATTATTTCCCTTTGATTATCAAGTATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATg  <  1:1556836/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATG  >  NC_000913/3437761‑3438059

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: