Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,042,958 T→C 5.6% Q418R (CAG→CGG)  etk ← protein‑tyrosine kinase Etk

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,042,9580TC 5.6% 217.8 / 2.9 71Q418R (CAG→CGG) etkprotein‑tyrosine kinase Etk
Reads supporting (aligned to +/- strand):  ref base T (35/32);  new base C (2/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

CATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATACCGCACGGCCCGCTTCTACGTCACGACTTAAACGCAACACTT  >  NC_000913/1042809‑1043099
                                                                                                                                                     |                                                                                                                                             
cATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCt                                                                                                                                               <  1:292705/150‑1 (MQ=255)
cATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCt                                                                                                                                               >  1:1288840/1‑150 (MQ=255)
 aTAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTg                                                                                                                                              >  1:108191/1‑150 (MQ=255)
    aaCGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCgg                                                                                                                                           >  2:491472/1‑150 (MQ=255)
    aaCGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCgg                                                                                                                                           >  2:20301/1‑150 (MQ=255)
     aCGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGc                                                                                                                                          >  2:255486/1‑150 (MQ=255)
       gCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTg                                                                                                                                        <  2:1183487/150‑1 (MQ=255)
         tGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAg                                                                                                                                      <  1:827509/150‑1 (MQ=255)
          gATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGt                                                                                                                                     >  2:681856/1‑150 (MQ=255)
          gATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGt                                                                                                                                     >  2:724773/1‑150 (MQ=255)
          gATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGt                                                                                                                                     <  2:500049/150‑1 (MQ=255)
           aTGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTg                                                                                                                                    <  2:1357655/150‑1 (MQ=255)
            tGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGa                                                                                                                                   >  1:730527/1‑150 (MQ=255)
             gCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGAc                                                                                                                                  >  1:792338/1‑150 (MQ=255)
                 tGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGcc                                                                                                                              <  2:984701/150‑1 (MQ=255)
                       tCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCg                                                                                                                        >  1:150730/1‑150 (MQ=255)
                            ttgttCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAAt                                                                                                                   <  2:849042/150‑1 (MQ=255)
                              gttCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATcc                                                                                                                 <  1:20301/150‑1 (MQ=255)
                                        ttCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTAcc                                                                                                       >  1:1045553/1‑150 (MQ=255)
                                          cTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGCTACCaa                                                                                                     >  2:1409862/1‑150 (MQ=255)
                                            acacCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATg                                                                                                   >  2:509751/1‑150 (MQ=255)
                                            acacCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATg                                                                                                   <  2:545789/150‑1 (MQ=255)
                                                cacGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCAc                                                                                               >  2:621304/1‑150 (MQ=255)
                                                  cgacgCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTg                                                                                             <  1:1012155/150‑1 (MQ=255)
                                                        aaCATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTc                                                                                       <  1:1355026/150‑1 (MQ=255)
                                                               cACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATAc                                                                                >  1:1036854/1‑150 (MQ=255)
                                                               cACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATAc                                                                                >  2:1186727/1‑150 (MQ=255)
                                                                     ccAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACt                                                                          >  2:391758/1‑150 (MQ=255)
                                                                       aGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATg                                                                                                   >  1:550148/1‑123 (MQ=255)
                                                                       aGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATg                                                                                                   <  2:550148/123‑1 (MQ=255)
                                                                       aGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTcc                                                                        <  1:243993/150‑1 (MQ=255)
                                                                       aGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTcc                                                                        >  2:1144910/1‑150 (MQ=255)
                                                                        gCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCct                                                                       >  1:413249/1‑150 (MQ=255)
                                                                          aCGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCctgc                                                                     >  1:1295553/1‑150 (MQ=255)
                                                                            ggCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCctgctg                                                                   <  1:255486/150‑1 (MQ=255)
                                                                             gCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTgg                                                                  >  2:651645/1‑150 (MQ=255)
                                                                              cACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGc                                                                 >  2:445674/1‑150 (MQ=255)
                                                                                   acaGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTa                                                            <  2:792338/150‑1 (MQ=255)
                                                                                   acaGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTa                                                            <  1:1334230/150‑1 (MQ=255)
                                                                                     aGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAg                                                          <  1:491472/150‑1 (MQ=255)
                                                                                      gAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACt                                                                               >  1:631133/1‑128 (MQ=255)
                                                                                      gAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACt                                                                               <  2:631133/128‑1 (MQ=255)
                                                                                      gAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGt                                                         <  1:509751/150‑1 (MQ=255)
                                                                                      gAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGt                                                         >  1:1366945/1‑150 (MQ=255)
                                                                                        aaTAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACg                                                                                                            >  1:475172/1‑97 (MQ=255)
                                                                                        aaTAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACg                                                                                                            <  2:475172/97‑1 (MQ=255)
                                                                                          tAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGAACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAAtt                                                     >  2:49921/1‑150 (MQ=255)
                                                                                              cAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAg                                                 <  2:1288840/150‑1 (MQ=255)
                                                                                              cAGGCCAAGAATAAAACCAAGCACCACGTCCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAg                                                 >  2:900201/1‑150 (MQ=255)
                                                                                                     aGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATAcc                                          >  2:292565/1‑150 (MQ=255)
                                                                                                     aGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATAcc                                          <  1:621304/150‑1 (MQ=255)
                                                                                                        aTAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATACCGCa                                       <  1:136058/150‑1 (MQ=255)
                                                                                                                    accacGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATACCGCACGGCCCGCTTCt                           >  1:1115852/1‑150 (MQ=255)
                                                                                                                     ccacGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATACCGCACGGCCCGCTTCTa                          >  2:83036/1‑150 (MQ=255)
                                                                                                                     ccacGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATACCGCACGGCCCGCTTCTa                          >  1:1109210/1‑150 (MQ=255)
                                                                                                                     ccacGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATACCGCACGGCCCGCTTCTa                          <  1:651645/150‑1 (MQ=255)
                                                                                                                     ccacGTTCAACGCTTTTTTCGGTTTCACTGGCCGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATACCGCACGGCCCGCTTCTa                          >  1:314630/1‑150 (MQ=255)
                                                                                                                     ccacGTCCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATACCGCACGGCCCGCTTCTa                          <  1:900201/150‑1 (MQ=255)
                                                                                                                       acGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATACCGCACGGCCCGCTTCTACg                        <  1:391758/150‑1 (MQ=255)
                                                                                                                         gTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATACCGCACGGCCCGCTTCTACGTc                      >  2:1393730/1‑150 (MQ=255)
                                                                                                                              aCGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATACCGCACGGCCCGCTTCTACGTCACGAc                 <  2:1235780/150‑1 (MQ=255)
                                                                                                                              aCGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATACCGCACGGCCCGCTTCTACGTCACGAc                 >  1:611885/1‑150 (MQ=255)
                                                                                                                                gCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATACCGCACGGCCCGCTTCTACGTCACGACtt               <  2:730527/150‑1 (MQ=255)
                                                                                                                                 cTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATACCGCACGGCCCGCTTCTACGTCACGACTTa              >  1:898075/1‑150 (MQ=255)
                                                                                                                                    ttttCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATACCGCACGGCCCGCTTCTACGTCACGACTTAAAc           <  1:724773/150‑1 (MQ=255)
                                                                                                                                     tttCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATACCGCACGGCCCGCTTCTACGTCACGACTTAAACg          <  1:445674/150‑1 (MQ=255)
                                                                                                                                     tttCGGTTTCACTGGCCGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATACCGCACGGCCCGCTTCTACGTCACGACTTAAACg          <  2:314630/150‑1 (MQ=255)
                                                                                                                                         ggTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATACCGCACGGCCCGCTTCTACGTCACGACTTAAACGCAac      <  2:108191/150‑1 (MQ=255)
                                                                                                                                          gTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTCACAATTGCACTGGATTTCGAAATACTCAACTCCTTCTGGCGTTTAAGTAATTGCAGATTTTCCGCTCTGCGCGCTTCTACGTCACTACTTAAACGCAaca     <  2:1248694/150‑1 (MQ=255)
                                                                                                                                          gTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATACCGCACGGCCCGCTTCTACGTCACGACTTAAACGCAaca     <  1:681856/150‑1 (MQ=255)
                                                                                                                                            ttCACTGGCCGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATACCGCACGGCCCGCTTCTACGTCACGACTTAACCGCaa       <  2:957228/146‑1 (MQ=255)
                                                                                                                                            ttCACTGGCCGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATACCGCACGGCCCGCTTCTACGTCACGACTTAACCGCaa       >  1:957228/1‑146 (MQ=255)
                                                                                                                                             tCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATATCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATACAGCACGGCCCGCTTCTACGCCACGACTTAAACGCAACACtt  >  2:956499/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
CATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCTTCTACACCACGACGCAACATCGCACGCGCCAGCACGGCACCCACAGAAATAAACAGGCCAAGAATAAAACCAAGCACCACGTTCAACGCTTTTTTCGGTTTCACTGGCTGCGGCTGAGTGACTGCCGGGTCGATAATCCGCACGTTACCAATGGCACTGGATTTCGAAATACTCAACTCCTGCTGGCGGTTAAGTAATTGCAGATATACCGCACGGCCCGCTTCTACGTCACGACTTAAACGCAACACTT  >  NC_000913/1042809‑1043099

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