Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,170,762 T→A 5.1% Q1068L (CAG→CTG) ‡ mfd ← transcription‑repair coupling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,170,7620TA 5.1% 523.9 / 10.9 150Q1068L (CAG→CTG) ‡mfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base T (80/62);  new base A (4/4);  total (84/66)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 152 total aligned reads displayed.

CGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCAC  >  NC_000913/1170614‑1170907
                                                                                                                                                    |                                                                                                                                                 
cgTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCt                                                                                                                                                   <  4:874274/149‑1 (MQ=255)
cgTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCt                                                                                                                                                   <  1:665478/149‑1 (MQ=255)
 gTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTg                                                                                                                                                  >  2:201555/1‑149 (MQ=255)
     gACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGgcgc                                                                                                                                              >  1:673442/1‑149 (MQ=255)
         aTCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGAc                                                                                                                                          >  2:346697/1‑149 (MQ=255)
           cAAGGCGGTAATTCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCCAGCAAGCGCAGaaac                                                                                                                                        <  2:100882/149‑4 (MQ=255)
           cAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAATCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTATCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACgg                                                                                                                                        <  3:738203/149‑1 (MQ=255)
           cAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACgg                                                                                                                                        <  4:37727/149‑1 (MQ=255)
           cAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGACTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCCAGCAAGCGCAGaaac                                                                                                                                        <  4:750703/149‑4 (MQ=255)
             aGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCCAGCAAGCGCAGaaactg                                                                                                                                      <  1:14078/149‑6 (MQ=255)
                     aTGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCca                                                                                                                              >  3:683680/1‑149 (MQ=255)
                          gcGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAgg                                                                                                                         <  1:105683/149‑1 (MQ=255)
                          gcGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAgg                                                                                                                         <  1:499987/149‑1 (MQ=255)
                            ggCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGt                                                                                                                       <  4:720150/149‑1 (MQ=255)
                               tGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTAcg                                                                                                                    >  1:104162/1‑149 (MQ=255)
                                        agcaAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATc                                                                                                           >  4:1045862/1‑149 (MQ=255)
                                        agcaAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATc                                                                                                           <  3:356729/149‑1 (MQ=255)
                                        agcaAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATc                                                                                                           >  2:450026/1‑149 (MQ=255)
                                           aaaCCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTgg                                                                                                        <  3:162132/149‑1 (MQ=255)
                                              ccAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGgcag                                                                                                     <  4:333296/149‑1 (MQ=255)
                                                aaTCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGgcagca                                                                                                   >  3:861760/1‑149 (MQ=255)
                                                aaTCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGgcagca                                                                                                   >  3:374576/1‑149 (MQ=255)
                                                 aTCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGcagcag                                                                                                  >  3:655958/1‑149 (MQ=255)
                                                    aaCCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGcc                                                                                               >  2:162133/1‑149 (MQ=255)
                                                       cAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGaa                                                                                            <  1:201555/149‑1 (MQ=255)
                                                       cAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGaa                                                                                            <  2:415352/149‑1 (MQ=255)
                                                       cAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGaa                                                                                            <  1:792075/149‑1 (MQ=255)
                                                         ggCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGc                                                                                          >  2:261968/1‑149 (MQ=255)
                                                            cGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcgat                                                                                       <  4:952810/149‑1 (MQ=255)
                                                             ggATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACgg                                                                                                                                        <  2:578380/99‑1 (MQ=255)
                                                             ggATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACgg                                                                                                                                        >  1:578380/1‑99 (MQ=255)
                                                                  aaCGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATaa                                                                                 <  1:871427/149‑1 (MQ=255)
                                                                   aCGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAg                                                                                >  1:178752/1‑149 (MQ=255)
                                                                      tGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTc                                                                             >  2:343790/1‑149 (MQ=255)
                                                                            ttctCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACttt                                                                       >  3:299102/1‑149 (MQ=255)
                                                                                cGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcgatcga                                                                                    >  1:930903/1‑132 (MQ=255)
                                                                                cGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcgatcga                                                                                    <  2:930903/132‑1 (MQ=255)
                                                                                cGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGAtc                                                                   <  4:655958/149‑1 (MQ=255)
                                                                                  gCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGAtctc                                                                 <  4:374576/149‑1 (MQ=255)
                                                                                   cAAATTCGATCACCCCGCCTTTCTCATTACCTTCCACCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATctct                                                                <  3:1067303/149‑1 (MQ=255)
                                                                                    aaaTTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCtt                                                               >  2:203254/1‑149 (MQ=255)
                                                                                     aaTTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAACATCCCGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTc                                                              >  3:782895/1‑149 (MQ=255)
                                                                                       ttCGATCACCCTGCCTTTCTCATTACCTTCCAGCTTCGTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCa                                                            <  3:703818/149‑1 (MQ=255)
                                                                                       ttCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTc                                                              <  1:415195/147‑1 (MQ=255)
                                                                                       ttCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTc                                                              >  2:415195/1‑147 (MQ=255)
                                                                                          gATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGtt                                                         >  3:387732/1‑149 (MQ=255)
                                                                                           aTCACCCCGCCTTTCTGATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCcagca                                                                                                                           <  4:622672/82‑1 (MQ=255)
                                                                                           aTCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCcagca                                                                                                                           >  3:622672/1‑82 (MQ=255)
                                                                                           aTCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTc                                                        >  2:475503/1‑149 (MQ=255)
                                                                                           aTCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCGGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATAGATAAGCTCGACTTTGATCTCTTCCAGTTc                                                        >  1:714387/1‑149 (MQ=255)
                                                                                              aCCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGtt                                                     >  3:887512/1‑149 (MQ=255)
                                                                                              aCCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGtt                                                     >  1:427791/1‑149 (MQ=255)
                                                                                               ccccGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGttt                                                    >  2:861576/1‑149 (MQ=255)
                                                                                               ccccGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGttt                                                    <  1:949334/149‑1 (MQ=255)
                                                                                                cccGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGc                                                                               <  4:618950/121‑1 (MQ=255)
                                                                                                cccGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGc                                                                               >  3:618950/1‑121 (MQ=255)
                                                                                                cccGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCAAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGtttt                                                   >  2:706639/1‑149 (MQ=255)
                                                                                                    ccTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGt                                               >  3:247401/1‑149 (MQ=255)
                                                                                                     cTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGtt                                              <  2:164508/149‑1 (MQ=255)
                                                                                                      tttCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCGCTTCCAGTTCGTTTTCCGttt                                             >  4:668345/1‑149 (MQ=255)
                                                                                                       ttCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGtttt                                            >  1:326479/1‑149 (MQ=255)
                                                                                                           cATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGgcg                                        >  1:109008/1‑149 (MQ=255)
                                                                                                            aTTACCTTCCAGCTTCTTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCCGTTCGTTTTCCGGTTTGgcgc                                       >  2:669173/1‑149 (MQ=255)
                                                                                                             ttACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCt                                      <  3:937184/149‑1 (MQ=255)
                                                                                                             ttACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCt                                      >  1:112542/1‑149 (MQ=255)
                                                                                                               aCCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTgg                                    <  4:587373/149‑1 (MQ=255)
                                                                                                                 cTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCCCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCa                                  >  1:528226/1‑149 (MQ=255)
                                                                                                                  tttcAGCTTCCTGATCCCTAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAGGCTCGACTTGGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCaa                                 >  4:85076/4‑149 (MQ=255)
                                                                                                                  ttCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCaa                                 <  3:574646/149‑1 (MQ=255)
                                                                                                                  ttCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCaa                                 >  2:791946/1‑149 (MQ=255)
                                                                                                                  ttCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCaa                                 >  1:580447/1‑149 (MQ=255)
                                                                                                                     cAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATAc                              >  4:1046234/1‑149 (MQ=255)
                                                                                                                      aGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACg                             >  4:1070757/1‑149 (MQ=255)
                                                                                                                      aGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACg                             >  4:117255/1‑149 (MQ=255)
                                                                                                                        cTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGtt                           >  3:198320/1‑149 (MQ=255)
                                                                                                                         ttCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGttt                          <  2:631851/149‑1 (MQ=255)
                                                                                                                         ttCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGttt                          <  2:1007791/149‑1 (MQ=255)
                                                                                                                         ttCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGttt                          <  4:247401/149‑1 (MQ=255)
                                                                                                                         ttCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGttt                          <  4:153015/149‑1 (MQ=255)
                                                                                                                         ttCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGttt                          <  1:450026/149‑1 (MQ=255)
                                                                                                                           ccTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGt                        >  1:252474/1‑149 (MQ=255)
                                                                                                                           ccTGATCCCCAGTTTCTGCACTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGt                        >  4:499822/1‑149 (MQ=255)
                                                                                                                             tGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGTGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAg                      >  3:358850/1‑149 (MQ=255)
                                                                                                                               aTCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcgatc                                                                                      <  3:479939/83‑1 (MQ=255)
                                                                                                                               aTCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcgatc                                                                                      >  4:479939/1‑83 (MQ=255)
                                                                                                                               aTCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGtat                          <  1:166925/143‑3 (MQ=255)
                                                                                                                               aTCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGaa                    <  2:231825/149‑1 (MQ=255)
                                                                                                                                  cccAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcgat                                                                                       >  1:138797/1‑79 (MQ=255)
                                                                                                                                   ccAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGAt                                                                    <  3:147408/97‑1 (MQ=255)
                                                                                                                                   ccAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTg                                           <  2:532699/122‑1 (MQ=255)
                                                                                                                                   ccAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTg                                           >  1:532699/1‑122 (MQ=255)
                                                                                                                                   ccAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGAc                <  4:200452/149‑1 (MQ=255)
                                                                                                                                   ccAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGAc                >  1:181779/1‑149 (MQ=255)
                                                                                                                                      ccgtCTGCGCCAGCAAGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGa             >  4:768439/4‑149 (MQ=255)
                                                                                                                                         tCTGCGCCAGCAAGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGAcgc          <  3:768439/149‑1 (MQ=255)
                                                                                                                                         tCTGCGCCAGCAAGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGAcgc          >  4:876899/1‑149 (MQ=255)
                                                                                                                                         tCTGCGCCAGCAAGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGAAAGAcgc          >  2:939534/1‑149 (MQ=255)
                                                                                                                                            gcgcTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCgtg       >  2:1055855/1‑149 (MQ=255)
                                                                                                                                               cTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTc    <  4:526444/149‑1 (MQ=255)
                                                                                                                                                 tGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCAc  <  2:1363/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCAC  >  NC_000913/1170614‑1170907

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