Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,437,911 T→A 7.1% intergenic (+17/‑113) 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,9110TA 7.1% 258.2 / 5.5 93intergenic (+17/‑113)trkA/mscLNAD‑binding component of TrK potassium transporter/mechanosensitive channel protein, high conductance
Reads supporting (aligned to +/- strand):  ref base T (46/40);  new base A (3/4);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGC  >  NC_000913/3437763‑3438060
                                                                                                                                                    |                                                                                                                                                     
ggAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATg                                                                                                                                                      <  1:1049411/150‑1 (MQ=255)
 gAAAAGAAGTGATGAATGCCAATGAAAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACATATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGa                                                                                                                                                     <  1:1184813/150‑1 (MQ=255)
  aaaCGACGTGATGATTGCCAATGACAATATGCTCATTGAGCAAGGCGATCAAGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGAc                                                                                                                                                    <  1:213970/150‑1 (MQ=255)
     cgacgTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGcc                                                                                                                                                 >  1:374181/1‑150 (MQ=255)
        cgTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCATTAGACGCCtta                                                                                                                                              <  1:71378/150‑1 (MQ=255)
        cgTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCATTAGACGCCtta                                                                                                                                              <  1:77786/150‑1 (MQ=255)
          tgatgaTTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCttatt                                                                                                                                            <  1:1126580/150‑1 (MQ=255)
           gatgatTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCttatta                                                                                                                                           <  2:711208/150‑1 (MQ=255)
             tgatTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTattatt                                                                                                                                         >  1:286327/1‑150 (MQ=255)
              tatTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTAttt                                                                                                                                        <  2:1137201/149‑1 (MQ=255)
                ttGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTcc                                                                                                                                      <  2:907257/150‑1 (MQ=255)
                 tGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTccc                                                                                                                                     <  2:343608/150‑1 (MQ=255)
                  gCCAATGACAATCTGCCCATTGAGCAAGGCGATCACTTTATTATGTTCCTCACAGCTAAAAAGTTTATTACCGACGGCGAAAGTCTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCt                                                                                                                                    <  2:94686/150‑1 (MQ=255)
                     aaTGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAAAGACGCCTTATTATTTCCCTTTg                                                                                                                                 <  1:246313/150‑1 (MQ=255)
                         aCAATCTGCGCATTGAGCAAGGCGAAAACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTACCGCCAAGCCATTTCTTCTTGTAATTAATAAGGCTTCTAATGACGCCTTATTTTTTCCCTTTGATTa                                                                                                                             >  1:618347/1‑150 (MQ=255)
                            aTCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCa                                                                                                                          <  1:172373/150‑1 (MQ=255)
                                  gcATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGTCGCCTTATTATTTCCCTTTGATTATCAAGGatt                                                                                                                    <  1:695523/150‑1 (MQ=255)
                                  gaaaagAGCAAGGCGATCACTTAATTTGGTTCCTCAAAGATTAAAAGTTTATTACCGACGTCGAAAGAATCTTCCAGCCAAGCCCTTTCAGCTAGTAAATAATAAGGCGTCTAATGACGCCTTATTAATTCCCTTTGATTATCAAGGatt                                                                                                                    <  2:184497/145‑1 (MQ=255)
                                   cATTGAGCAAGTCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCATTTCTTCTTGTAATTAATAACGCTACTAATGACTCCTTATTATTTCACTTTAATTATCAAgtaata                                                                                                                   >  2:287912/1‑145 (MQ=255)
                                     ttGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGattaat                                                                                                                 <  1:1212072/150‑1 (MQ=255)
                                     ttGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGattaat                                                                                                                 <  1:349718/150‑1 (MQ=255)
                                      tGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGattaatta                                                                                                               <  2:293241/150‑1 (MQ=255)
                                      tGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGattaatta                                                                                                               <  2:21469/150‑1 (MQ=255)
                                       gAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCACGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGattaatta                                                                                                               >  2:711406/1‑150 (MQ=255)
                                         gCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGTCGTCTAATGACGCCTTATTATTTCCCTGTGATTATCAAGGattatttatt                                                                                                             >  2:1307571/1‑148 (MQ=255)
                                             ggCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAattca                                                                                                         <  1:578037/150‑1 (MQ=255)
                                              gCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATAACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATTAAGCATTATTATTTAAATATGATTATAAAGGATTAATTAAattcat                                                                                                        >  2:570160/1‑150 (MQ=255)
                                               cGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACTCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAattcatt                                                                                                       <  1:636417/150‑1 (MQ=255)
                                                 aTCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCCGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTcc                                                                                                     >  2:1138749/1‑150 (MQ=255)
                                                   cACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTg                                                                                                   <  2:585215/150‑1 (MQ=255)
                                                    aCGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTgg                                                                                                  >  2:488291/1‑150 (MQ=255)
                                                       tAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAg                                                                                               <  2:1256393/150‑1 (MQ=255)
                                                           tATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGTCAGGaaa                                                                                           >  1:1174695/1‑150 (MQ=255)
                                                           tATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGaaa                                                                                           >  1:1180106/1‑150 (MQ=255)
                                                           tATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGaaa                                                                                           >  2:944406/1‑150 (MQ=255)
                                                           tATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCATTCTTATTGTAATTAATAAGGCTTCTAATGACGCATTATTATTTCCCTTTGATTAACAAGGATTAATTAAATACATTCCTGTCAGCaaa                                                                                           >  1:283713/1‑150 (MQ=255)
                                                              tttCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAATCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATg                                                                                        <  1:339614/149‑1 (MQ=255)
                                                               ttCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGTATTAATTAAATTCATTCCTGGCAGGAAAATgg                                                                                       >  1:1016607/1‑150 (MQ=255)
                                                                tCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTGTTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGc                                                                                      <  2:1075276/150‑1 (MQ=255)
                                                                      caGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACa                                                                                >  1:489489/1‑150 (MQ=255)
                                                                      aaGATAAAAAGTTTATAACCGAGGTCGAAAGACTTTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCGTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAATATGGCTTAACa                                                                                <  2:43821/149‑1 (MQ=255)
                                                                          tAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTg                                                                            >  2:333962/1‑150 (MQ=255)
                                                                           aaaaaGTTTATTACCGACGTCTAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGt                                                                           <  1:887874/150‑1 (MQ=255)
                                                                           aaaaaGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGt                                                                           >  2:750807/1‑150 (MQ=255)
                                                                            aaaaGTTTATTACCGTCGTCGAAAGACTCTTCCAGCCATGCCCTTTATTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGtt                                                                          <  2:20113/150‑1 (MQ=255)
                                                                             aaaGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTa                                                                         <  1:1267902/150‑1 (MQ=255)
                                                                             aaaGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTa                                                                         <  1:305133/150‑1 (MQ=255)
                                                                                 tttATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACt                                                                     <  2:284394/150‑1 (MQ=255)
                                                                                   tattaCCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTa                                                                   >  1:1083881/1‑150 (MQ=255)
                                                                                      taCCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGAAGCCATATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATATGTTAGACTTATgt                                                                >  1:482025/1‑149 (MQ=255)
                                                                                       aCCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTGATTATTTCCCTTTGATTATCAAGGATTAATTATATTAATTCCTGGCAGAAAAATGGCTTAACATTTGTTAGACTTATGGt                                                               >  1:189790/1‑150 (MQ=255)
                                                                                             gTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCgg                                                         <  2:733780/150‑1 (MQ=255)
                                                                                                gAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCtt                                                      <  2:405890/150‑1 (MQ=255)
                                                                                                gAAAGACTCTTCCAGCCAAGCCCTTTCCTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCtt                                                      >  2:271070/1‑150 (MQ=255)
                                                                                                 aaaGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTc                                                     >  1:288227/1‑150 (MQ=255)
                                                                                                  aaGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCa                                                    <  1:467729/150‑1 (MQ=255)
                                                                                                  aaGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCa                                                    <  2:97276/150‑1 (MQ=255)
                                                                                                      ctctTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTAGTCTTAATatt                                                >  1:844178/1‑148 (MQ=255)
                                                                                                        tttCGAGCCAAGCCCTTTCTTGTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGAGTTATGGTTGTCGGCTTCATAGGga                                              <  2:423964/149‑1 (MQ=255)
                                                                                                         ttCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGgag                                             >  2:154892/1‑150 (MQ=255)
                                                                                                          tCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGCCGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGgaga                                            >  2:159495/1‑150 (MQ=255)
                                                                                                              gCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATaa                                        >  1:495954/1‑150 (MQ=255)
                                                                                                               ccAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATaaa                                       >  2:1119255/1‑149 (MQ=255)
                                                                                                                cAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGCCGCCTTATTATTTCGCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACa                                      >  1:1127344/1‑150 (MQ=255)
                                                                                                                cAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGATTAACa                                      >  1:1220648/1‑150 (MQ=255)
                                                                                                                cAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACa                                      <  1:90511/150‑1 (MQ=255)
                                                                                                                cAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATAATCAATGATTTATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTATACTTATGGTTGTCGGCTTCATAGGGAGAATAACa                                      <  2:72863/150‑1 (MQ=255)
                                                                                                                 aaGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACAt                                     >  2:891872/1‑150 (MQ=255)
                                                                                                                  aGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATg                                    <  2:1291681/150‑1 (MQ=255)
                                                                                                                  aGCCCTTTCTTCTTGTAATTAATAAGGCGGCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATg                                    <  2:618347/150‑1 (MQ=255)
                                                                                                                    cccTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTGTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAg                                  >  1:844863/1‑150 (MQ=255)
                                                                                                                      cTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCa                                >  1:615682/1‑150 (MQ=255)
                                                                                                                          cttcttGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCattat                            >  2:210920/1‑150 (MQ=255)
                                                                                                                           ttcttGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCattatt                           >  2:1082092/1‑150 (MQ=255)
                                                                                                                            tcttGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCAttatta                          <  1:112719/150‑1 (MQ=255)
                                                                                                                               tGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAg                       >  1:1084010/1‑150 (MQ=255)
                                                                                                                                 taattaatAAGGCGTCGAATGACGCCTTATTATTTCCCTTTGATTATCAATGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGaa                     >  1:90745/1‑150 (MQ=255)
                                                                                                                                  aattaatAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAAt                    >  2:689402/1‑150 (MQ=255)
                                                                                                                                      aataaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTcg                <  1:488291/150‑1 (MQ=255)
                                                                                                                                       ataaGGCGTCATTAGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTATATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGGTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTcgc               >  2:610193/1‑150 (MQ=255)
                                                                                                                                       ataaGGCGTCATTAGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTcgc               >  2:137128/1‑150 (MQ=255)
                                                                                                                                       ataaGGCGTCATTAGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTagc               >  1:37782/1‑150 (MQ=255)
                                                                                                                                        taacgagTCTAAAGACGCCTAATTATTTCCCTTTGATTTTCAATGATTAATTAAATTCAATCCTGGCACGAAAATGGCTTAACATTTGTAAGACTTATGGCTGTCGGCTTCATATGGAGAATAACATGAGCATTATTAAAGAATTTcgcg              <  1:115595/144‑1 (MQ=255)
                                                                                                                                        taaGGCGTCTAATGCCGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATACCATGAGCATTATTAAAGAATTTcgcg              >  2:938398/1‑150 (MQ=255)
                                                                                                                                         aaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCTCGa             >  1:941840/1‑150 (MQ=255)
                                                                                                                                         aaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGCGAGAATAACATGAGCATTATTAAAGAATTTCGCGa             >  1:69217/1‑150 (MQ=255)
                                                                                                                                         aaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATAAAATTCATTCCTGGAAGGAAAATCGCTTAACATTTCTTAGACTAATGGTTGTCGTCTTCATAGGGAGAATAAAATGATCATTATTAAAGAATATCTCGa             >  1:582458/1‑150 (MQ=255)
                                                                                                                                           ggCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAAt           >  2:894608/1‑150 (MQ=255)
                                                                                                                                             cGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAAttt         >  2:304412/1‑150 (MQ=255)
                                                                                                                                              gTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTg        >  2:1164180/1‑150 (MQ=255)
                                                                                                                                                cTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTAGCGAATTTGcc      >  2:489329/1‑149 (MQ=255)
                                                                                                                                                  aaTGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGATTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATACCATGAGCATAATTAAAGAATTTCGCGAATTTGCGAt    >  2:633892/1‑150 (MQ=255)
                                                                                                                                                   aTGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATg   <  1:615775/150‑1 (MQ=255)
                                                                                                                                                    tGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATgc  >  2:578012/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATGC  >  NC_000913/3437763‑3438060

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: