Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,437,911 T→A 5.0% 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 5.0% 404.8 / 5.8 131intergenic (+17/‑113)trkA/mscLNAD‑binding component of TrK potassium transporter/mechanosensitive channel protein, high conductance
Reads supporting (aligned to +/- strand):  ref base T (52/72);  new base A (2/5);  total (54/77)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 133 total aligned reads displayed.

GGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATG  >  NC_000913/3437763‑3438059
                                                                                                                                                    |                                                                                                                                                    
ggAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAAt                                                                                                                                                      <  1:806942/149‑1 (MQ=255)
  aaaCGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGTCGTCTAATGa                                                                                                                                                    >  2:1266124/1‑149 (MQ=255)
     cgacgTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGc                                                                                                                                                 <  2:1055642/149‑1 (MQ=255)
        cgTTATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCAAGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCtt                                                                                                                                              <  1:71664/149‑1 (MQ=255)
         gTGATGATTGCCAATGACAATCTGCTCATTGAGCAAGGCTATCAAGTAATTATGTTACTCACAGATAAAAAGTTTATTACCGAAGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCtta                                                                                                                                             <  1:315462/149‑1 (MQ=255)
         gTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCtta                                                                                                                                             >  1:12983/1‑149 (MQ=255)
         gTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCtta                                                                                                                                             <  1:1510089/149‑1 (MQ=255)
         gTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCATTAGACGCCtta                                                                                                                                             <  1:818692/149‑1 (MQ=255)
         gTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCATTAGACGCCtta                                                                                                                                             <  1:813086/149‑1 (MQ=255)
         gTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCATTAGACGCCtta                                                                                                                                             <  1:146463/149‑1 (MQ=255)
           gatgatTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGATTATTACCGACGTCGAAAGACTCTTCCAGCCCAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCttatt                                                                                                                                           <  1:785044/149‑1 (MQ=255)
            atgatTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCttatta                                                                                                                                          >  1:528642/1‑149 (MQ=255)
            agcacagCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCttatta                                                                                                                                          <  2:797194/143‑1 (MQ=255)
               atTGCCAATGACAATCTGAGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATATGGCGTCTAATGACGCCTTATTAttt                                                                                                                                       <  2:877274/149‑1 (MQ=255)
                  gacAATGACAATCTGCGCATTTAGCAAGGCTATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTccc                                                                                                                                    <  2:1012374/147‑1 (MQ=255)
                  gCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATAAATAAGGCGTCTAATGACGCCTTATTATTTccc                                                                                                                                    <  1:332625/149‑1 (MQ=255)
                   ccAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCt                                                                                                                                   >  2:1464149/1‑149 (MQ=255)
                   ccAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCt                                                                                                                                   <  2:1411138/149‑1 (MQ=255)
                   ccAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCt                                                                                                                                   <  2:1545568/149‑1 (MQ=255)
                   ccAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCt                                                                                                                                   <  2:300806/149‑1 (MQ=255)
                   ccAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCt                                                                                                                                   >  2:1142018/1‑149 (MQ=255)
                   ccAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCt                                                                                                                                   >  1:1495627/1‑149 (MQ=255)
                    cAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCtt                                                                                                                                  <  1:948921/149‑1 (MQ=255)
                     aaTGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCttt                                                                                                                                 <  2:1404480/149‑1 (MQ=255)
                     aaTGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCttt                                                                                                                                 >  2:1340752/1‑149 (MQ=255)
                     aaTGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCttt                                                                                                                                 <  1:427494/149‑1 (MQ=255)
                      aTGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCATTAGACGCCTTATTAattacaagaa                                                                                                                                <  1:801795/149‑8 (MQ=38)
                               tGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAg                                                                                                                       <  2:595966/149‑1 (MQ=255)
                               tGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAg                                                                                                                       <  2:538230/149‑1 (MQ=255)
                                gcgcATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAgg                                                                                                                      <  2:989164/149‑1 (MQ=255)
                                gcgcATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAgg                                                                                                                      >  2:439984/1‑149 (MQ=255)
                                  tcATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGat                                                                                                                    <  2:200870/148‑1 (MQ=255)
                                  gcATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGat                                                                                                                    <  2:1440377/149‑1 (MQ=255)
                                    aTTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGatta                                                                                                                  >  2:1411337/1‑149 (MQ=255)
                                          cAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTaaa                                                                                                            <  2:875251/149‑1 (MQ=255)
                                             ggCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAattc                                                                                                         >  1:158242/1‑149 (MQ=255)
                                                gATCACGTAATTACGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAattcatt                                                                                                      >  2:936378/1‑149 (MQ=255)
                                                 aTCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAattcattc                                                                                                     <  1:952154/149‑1 (MQ=255)
                                                   cACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCt                                                                                                   <  2:516985/149‑1 (MQ=255)
                                                    aCGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTTCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTg                                                                                                  >  2:1445046/1‑149 (MQ=255)
                                                    aCGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTg                                                                                                  <  2:588856/149‑1 (MQ=255)
                                                    aCGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTg                                                                                                  >  2:1166196/1‑149 (MQ=255)
                                                      gTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGc                                                                                                >  1:871235/1‑149 (MQ=255)
                                                         aTTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAgg                                                                                             >  2:811636/1‑149 (MQ=255)
                                                             tGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGaaaa                                                                                         <  1:318154/149‑1 (MQ=255)
                                                              gTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAAt                                                                                        >  2:422919/1‑149 (MQ=255)
                                                                tCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATgg                                                                                      <  1:900085/149‑1 (MQ=255)
                                                                  cTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCt                                                                                    <  1:1014341/149‑1 (MQ=255)
                                                                  cTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCt                                                                                    <  2:322700/149‑1 (MQ=255)
                                                                  cTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCt                                                                                    <  1:949321/149‑1 (MQ=255)
                                                                  cTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCt                                                                                    <  1:1454237/149‑1 (MQ=255)
                                                                   tCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCtt                                                                                   >  1:716208/1‑149 (MQ=255)
                                                                        gATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACAt                                                                              <  1:389810/149‑1 (MQ=255)
                                                                         aTAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACAtt                                                                             >  2:924529/1‑149 (MQ=255)
                                                                           aaaaaGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTg                                                                           <  1:659568/149‑1 (MQ=255)
                                                                           aaaaaGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTg                                                                           <  2:547336/149‑1 (MQ=255)
                                                                             aaaGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGtt                                                                         <  2:1478459/149‑1 (MQ=255)
                                                                              aaGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTa                                                                        <  1:1142018/149‑1 (MQ=255)
                                                                                gTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTAATTACTGGCAGGAAAATTGCTTAACATTTGTTAGa                                                                      >  2:389755/1‑149 (MQ=255)
                                                                                     ttaCCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTAt                                                                 <  2:1510754/149‑1 (MQ=255)
                                                                                     ttaCCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTAt                                                                 <  2:990027/149‑1 (MQ=255)
                                                                                     ttaCCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTAt                                                                 <  2:1120985/149‑1 (MQ=255)
                                                                                      taCCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATg                                                                <  1:363394/149‑1 (MQ=255)
                                                                                      taCCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATg                                                                <  2:871235/149‑1 (MQ=255)
                                                                                         cGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGtt                                                             <  2:945531/149‑1 (MQ=255)
                                                                                           aCGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGt                                                           >  1:465562/1‑149 (MQ=255)
                                                                                           aCGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGt                                                           <  2:158242/149‑1 (MQ=255)
                                                                                            cGTCGAAAGACTCATCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTc                                                          <  2:25377/149‑1 (MQ=255)
                                                                                             gTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCg                                                         <  2:716208/149‑1 (MQ=255)
                                                                                                   aGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCa                                                   <  1:1416847/149‑1 (MQ=255)
                                                                                                       tctTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAgg                                               >  2:336108/1‑149 (MQ=255)
                                                                                                         ttCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGga                                             <  1:426685/149‑1 (MQ=255)
                                                                                                         ttCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGga                                             <  1:936378/149‑1 (MQ=255)
                                                                                                         ttCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTAAGGTTGTCGGCTTCATAGGga                                             >  2:987629/1‑149 (MQ=255)
                                                                                                              gCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATa                                        <  1:1411337/149‑1 (MQ=255)
                                                                                                               ccAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATaa                                       >  2:788186/1‑149 (MQ=255)
                                                                                                                cAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAAc                                      >  1:1261921/1‑149 (MQ=255)
                                                                                                                 aatcccTTTATTCTTGTAATTAATAAGGAGTCTAATGACTCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACa                                     <  2:169364/146‑1 (MQ=255)
                                                                                                                 aaGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACa                                     <  2:1205569/149‑1 (MQ=255)
                                                                                                                     ccTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAg                                 <  2:528642/149‑1 (MQ=255)
                                                                                                                       tttCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTATATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGATTCATAGGGAGAATAACATGAGCa                               >  1:111799/1‑149 (MQ=255)
                                                                                                                       tttCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCa                               >  2:430541/1‑149 (MQ=255)
                                                                                                                            tcttGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCattatt                          >  2:70211/1‑149 (MQ=255)
                                                                                                                               tGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTaaa                       <  2:1511989/149‑1 (MQ=255)
                                                                                                                               tGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTaaa                       >  2:358659/1‑149 (MQ=255)
                                                                                                                               tGTAATTAATAAGGCGTATAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTaaa                       <  1:1266124/149‑1 (MQ=255)
                                                                                                                                ggaaataataaGGCGTCATTAGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAAc                                      >  2:209221/6‑133 (MQ=39)
                                                                                                                                ggaaataataaGGCGTCATTAGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAAc                                      <  1:209221/128‑1 (MQ=39)
                                                                                                                                  aattaatAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGaa                    <  2:1398328/149‑1 (MQ=255)
                                                                                                                                  aattaatAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGaa                    >  2:720676/1‑149 (MQ=255)
                                                                                                                                     taataaGGCGTCATTAGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGATTAACATGAGCATTATTAAAGAAttt                 >  1:214961/1‑149 (MQ=255)
                                                                                                                                       ataaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTcg               >  2:1367510/1‑149 (MQ=255)
                                                                                                                                       ataaGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTcg               >  1:537550/1‑149 (MQ=255)
                                                                                                                                          aGGCATCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAAAATGAGCATTATTAAAGAATTTCGCGa            >  2:211504/1‑149 (MQ=255)
                                                                                                                                           ggCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAAAATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGaa           >  2:1595382/1‑149 (MQ=255)
                                                                                                                                            gCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAAt          >  1:774521/1‑149 (MQ=255)
                                                                                                                                               tCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTg       <  1:358659/149‑1 (MQ=255)
                                                                                                                                                   aTGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGAt   <  1:689981/149‑1 (MQ=255)
                                                                                                                                                   aTGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGAt   >  1:1171710/1‑149 (MQ=255)
                                                                                                                                                    tGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATg  >  2:1072917/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGAAACGACGTGATGATTGCCAATGACAATCTGCGCATTGAGCAAGGCGATCACGTAATTATGTTCCTCACAGATAAAAAGTTTATTACCGACGTCGAAAGACTCTTCCAGCCAAGCCCTTTCTTCTTGTAATTAATAAGGCGTCTAATGACGCCTTATTATTTCCCTTTGATTATCAAGGATTAATTAAATTCATTCCTGGCAGGAAAATGGCTTAACATTTGTTAGACTTATGGTTGTCGGCTTCATAGGGAGAATAACATGAGCATTATTAAAGAATTTCGCGAATTTGCGATG  >  NC_000913/3437763‑3438059

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