Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,170,758 T→C 5.2% Q1069Q (CAA→CAG) ‡ mfd ← transcription‑repair coupling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,170,7580TC 5.2% 509.7 / 9.9 152Q1069Q (CAA→CAG) ‡mfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base T (79/65);  new base C (4/4);  total (83/69)
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.93e-01

Only 100 of 155 total aligned reads displayed.

GCGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCA  >  NC_000913/1170611‑1170906
                                                                                                                                                   |                                                                                                                                                    
gcgcGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCtt                                                                                                                                                     <  3:1041993/149‑1 (MQ=255)
   cgTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCt                                                                                                                                                  <  4:874274/149‑1 (MQ=255)
    gTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTg                                                                                                                                                 >  3:853553/1‑149 (MQ=255)
    gTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTg                                                                                                                                                 >  2:201555/1‑149 (MQ=255)
        gACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGgcgc                                                                                                                                             >  1:673442/1‑149 (MQ=255)
            aTCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGAc                                                                                                                                         >  2:346697/1‑149 (MQ=255)
             tCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACg                                                                                                                                        >  3:685062/1‑149 (MQ=255)
              cAAGGCGGTAATTCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCCAGCAAGCGCAGaaac                                                                                                                                       <  2:100882/149‑4 (MQ=255)
              cAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACgg                                                                                                                                       <  4:37727/149‑1 (MQ=255)
              cAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGACTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCCAGCAAGCGCAGaaac                                                                                                                                       <  4:750703/149‑4 (MQ=255)
                aGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCCAGCAAGCGCAGaaactg                                                                                                                                     <  1:14078/149‑6 (MQ=255)
                             gcGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAgg                                                                                                                        <  1:105683/149‑1 (MQ=255)
                             gcGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAgg                                                                                                                        <  1:542671/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                                                                                                          >  2:450026/1‑149 (MQ=255)
                                           agcaAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATc                                                                                                          <  3:356729/149‑1 (MQ=255)
                                            gcaAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCt                                                                                                         <  3:373560/149‑1 (MQ=255)
                                             caAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTg                                                                                                        >  1:631851/1‑149 (MQ=255)
                                              aaaCCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTgg                                                                                                       <  3:162132/149‑1 (MQ=255)
                                                aCCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGgca                                                                                                     >  2:521096/1‑149 (MQ=255)
                                                 ccAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGgcag                                                                                                    <  4:333296/149‑1 (MQ=255)
                                                   aaTCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGgcagca                                                                                                  >  4:370921/1‑149 (MQ=255)
                                                   aaTCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGgcagca                                                                                                  >  3:374576/1‑149 (MQ=255)
                                                   aaTCAACCAGGCAGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGgcagca                                                                                                  >  4:106520/1‑149 (MQ=255)
                                                    aTCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGcagcag                                                                                                 >  3:655958/1‑149 (MQ=255)
                                                       aaCCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGcc                                                                                              >  2:792075/1‑149 (MQ=255)
                                                       aaCCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGcc                                                                                              >  2:162133/1‑149 (MQ=255)
                                                          cAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGaa                                                                                           <  1:792075/149‑1 (MQ=255)
                                                          cAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGaa                                                                                           <  1:201555/149‑1 (MQ=255)
                                                            ggCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGc                                                                                         >  2:261968/1‑149 (MQ=255)
                                                            ggCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGc                                                                                         <  1:346697/149‑1 (MQ=255)
                                                            ggCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGc                                                                                         >  4:414756/1‑149 (MQ=255)
                                                             gCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcg                                                                                        >  3:751115/1‑149 (MQ=255)
                                                               cGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcgat                                                                                      <  3:985871/149‑1 (MQ=255)
                                                               cGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcgat                                                                                      <  4:952810/149‑1 (MQ=255)
                                                                ggATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACgg                                                                                                                                       >  1:578380/1‑99 (MQ=255)
                                                                     aaCGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATaa                                                                                <  1:871427/149‑1 (MQ=255)
                                                                      aCGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAg                                                                               >  1:178752/1‑149 (MQ=255)
                                                                      aCGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAg                                                                               >  3:44775/1‑149 (MQ=255)
                                                                        gTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCt                                                                             >  4:604180/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)
                                                                                   cGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGAtc                                                                  <  4:655958/149‑1 (MQ=255)
                                                                                    ggCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGAtct                                                                 >  2:429268/1‑149 (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)
                                                                                              aTCACCCCGCCTTTCTGATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCcagca                                                                                                                          <  4:622672/82‑1 (MQ=255)
                                                                                              aTCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCGGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATAGATAAGCTCGACTTTGATCTCTTCCAGTTc                                                       >  1:714387/1‑149 (MQ=255)
                                                                                                 aCCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGtt                                                    >  3:887512/1‑149 (MQ=255)
                                                                                                  ccccGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGttt                                                   <  1:949334/149‑1 (MQ=255)
                                                                                                  ccccGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGttt                                                   >  2:861576/1‑149 (MQ=255)
                                                                                                   cccGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGc                                                                              >  3:618950/1‑121 (MQ=255)
                                                                                                   cccGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGtttt                                                  >  4:522180/1‑149 (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)
                                                                                                           tctcATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTg                                          <  1:861576/149‑1 (MQ=255)
                                                                                                              cATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGgcg                                       >  1:109008/1‑149 (MQ=255)
                                                                                                                ttACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCt                                     >  1:112542/1‑149 (MQ=255)
                                                                                                                  aCCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTc                                                                            <  3:651254/108‑1 (MQ=255)
                                                                                                                  aCCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTgg                                   <  4:587373/149‑1 (MQ=255)
                                                                                                                    cTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCCCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCa                                 >  1:528226/1‑149 (MQ=255)
                                                                                                                     tttcAGCTTCCTGATCCCTAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAGGCTCGACTTGGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCaa                                >  4:85076/4‑149 (MQ=255)
                                                                                                                     ttCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCaa                                >  2:690381/1‑149 (MQ=255)
                                                                                                                         aGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACg                            >  4:1070757/1‑149 (MQ=255)
                                                                                                                           cTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGtt                          <  1:429268/149‑1 (MQ=255)
                                                                                                                           cTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGtt                          >  3:198320/1‑149 (MQ=255)
                                                                                                                            ttCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGttt                         <  1:706639/149‑1 (MQ=255)
                                                                                                                            ttCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGttt                         <  1:450026/149‑1 (MQ=255)
                                                                                                                            ttCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGttt                         <  3:916751/149‑1 (MQ=255)
                                                                                                                            ttCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGttt                         <  2:631851/149‑1 (MQ=255)
                                                                                                                            ttCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGttt                         <  2:1007791/149‑1 (MQ=255)
                                                                                                                              ccTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGt                       >  1:252474/1‑149 (MQ=255)
                                                                                                                              ccTGATCCCCAGTTTCTGCACTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGt                       >  4:499822/1‑149 (MQ=255)
                                                                                                                               cTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGttt                                            <  4:209378/127‑1 (MQ=255)
                                                                                                                               cTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTa                      <  1:521096/149‑1 (MQ=255)
                                                                                                                                  aTCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcgatc                                                                                     <  3:479939/83‑1 (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)
                                                                                                                                      ccAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGAt                                                                   >  4:147408/1‑97 (MQ=255)
                                                                                                                                      ccAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTg                                          <  2:532699/122‑1 (MQ=255)
                                                                                                                                      ccAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGAc               >  1:181779/1‑149 (MQ=255)
                                                                                                                                      ccAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGAc               <  4:200452/149‑1 (MQ=255)
                                                                                                                                      ccAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGAc               <  2:580447/149‑1 (MQ=255)
                                                                                                                                         ccgtCTGCGCCAGCAAGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGa            >  2:919737/4‑149 (MQ=255)
                                                                                                                                         ccgtCTGCGCCAGCAAGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGa            >  4:768439/4‑149 (MQ=255)
                                                                                                                                            tCTGCGCCAGCAAGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGAcgc         >  4:876899/1‑149 (MQ=255)
                                                                                                                                            tCTGCGCCAGCAAGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGAAAGAcgc         >  2:939534/1‑149 (MQ=255)
                                                                                                                                              tGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCgt       <  3:604180/149‑1 (MQ=255)
                                                                                                                                               gcgcTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCgtg      >  2:1055855/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCA  >  NC_000913/1170611‑1170906

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