Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,170,759 T→A 5.2% Q1069L (CAA→CTA) ‡ mfd ← transcription‑repair coupling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,170,7590TA 5.2% 531.9 / 10.7 152Q1069L (CAA→CTA) ‡mfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base T (79/65);  new base A (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.82e-01

Only 100 of 156 total aligned reads displayed.

GCGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCAC  >  NC_000913/1170611‑1170907
                                                                                                                                                    |                                                                                                                                                    
gcgcGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCtt                                                                                                                                                      <  3:1041993/149‑1 (MQ=255)
   cgTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCt                                                                                                                                                   <  4:874274/149‑1 (MQ=255)
    gTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTg                                                                                                                                                  >  3:853553/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)
              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                                                                                                           <  3:356729/149‑1 (MQ=255)
                                             caAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTg                                                                                                         >  1:631851/1‑149 (MQ=255)
                                                   aaTCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGgcagca                                                                                                   >  3:861760/1‑149 (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:162133/1‑149 (MQ=255)
                                                         ccAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGa                                                                                             <  2:340163/149‑1 (MQ=255)
                                                          cAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGaa                                                                                            <  1:792075/149‑1 (MQ=255)
                                                          cAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGaa                                                                                            <  2:415352/149‑1 (MQ=255)
                                                            ggCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGc                                                                                          <  1:346697/149‑1 (MQ=255)
                                                            ggCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGc                                                                                          >  4:414756/1‑149 (MQ=255)
                                                            ggCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGc                                                                                          >  2:261968/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                                                                                                                                        <  2:578380/99‑1 (MQ=255)
                                                                      aCGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAg                                                                                >  1:178752/1‑149 (MQ=255)
                                                                        gTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCt                                                                              >  4:604180/1‑149 (MQ=255)
                                                                               ttctCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACttt                                                                       >  3:299102/1‑149 (MQ=255)
                                                                                   cGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcgatcga                                                                                    <  2:930903/132‑1 (MQ=255)
                                                                                   cGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcgatcga                                                                                    >  1:930903/1‑132 (MQ=255)
                                                                                    ggCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGAtct                                                                  >  2:429268/1‑149 (MQ=255)
                                                                                     gCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGAtctc                                                                 >  2:220210/1‑149 (MQ=255)
                                                                                     gCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGAtctc                                                                 >  1:246399/1‑149 (MQ=255)
                                                                                       aaaTTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCtt                                                               >  2:203254/1‑149 (MQ=255)
                                                                                          ttCGATCACCCTGCCTTTCTCATTACCTTCCAGCTTCGTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCa                                                            <  3:703818/149‑1 (MQ=255)
                                                                                          ttCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTc                                                              <  1:415195/147‑1 (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)
                                                                                                  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)
                                                                                                         tttCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCGCTTCCAGTTCGTTTTCCGttt                                             >  4:668345/1‑149 (MQ=255)
                                                                                                           tctcATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTg                                           <  1:861576/149‑1 (MQ=255)
                                                                                                             tcATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGgc                                         >  3:153015/1‑149 (MQ=255)
                                                                                                              cATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGgcg                                        >  1:109008/1‑149 (MQ=255)
                                                                                                               aTTACCTTCCAGCTTCTTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCCGTTCGTTTTCCGGTTTGgcgc                                       >  2:669173/1‑149 (MQ=255)
                                                                                                                ttACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCt                                      >  1:112542/1‑149 (MQ=255)
                                                                                                                 tACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTg                                     <  2:1004346/149‑1 (MQ=255)
                                                                                                                  aCCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTc                                                                             >  4:651254/1‑108 (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)
                                                                                                                         aGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACg                             >  4:1070757/1‑149 (MQ=255)
                                                                                                                           cTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGtt                           >  3:198320/1‑149 (MQ=255)
                                                                                                                           cTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGtt                           <  1:429268/149‑1 (MQ=255)
                                                                                                                            ttCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGttt                          <  2:631851/149‑1 (MQ=255)
                                                                                                                            ttCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGttt                          <  1:706639/149‑1 (MQ=255)
                                                                                                                            ttCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGttt                          <  4:247401/149‑1 (MQ=255)
                                                                                                                            ttCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGttt                          <  1:450026/149‑1 (MQ=255)
                                                                                                                            ttCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGttt                          <  3:916751/149‑1 (MQ=255)
                                                                                                                            ttCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGttt                          <  2:1007791/149‑1 (MQ=255)
                                                                                                                            ttCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGttt                          <  4:153015/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)
                                                                                                                               cTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGttt                                             >  3:209378/1‑127 (MQ=255)
                                                                                                                               cTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTa                       <  3:1045862/149‑1 (MQ=255)
                                                                                                                               cTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTa                       <  1:521096/149‑1 (MQ=255)
                                                                                                                                tGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGTGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAg                      >  3:358850/1‑149 (MQ=255)
                                                                                                                                  aTCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcgatc                                                                                      <  3:479939/83‑1 (MQ=255)
                                                                                                                                  aTCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGttt                          >  2:166925/1‑143 (MQ=255)
                                                                                                                                  aTCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGtat                          <  1:166925/143‑3 (MQ=255)
                                                                                                                                  aTCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTAGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGaa                    <  1:875110/149‑1 (MQ=255)
                                                                                                                                     cccAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGcgat                                                                                       <  2:138797/79‑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)
                                                                                                                                       cAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACa               <  2:104162/149‑1 (MQ=255)
                                                                                                                                         ccgtCTGCGCCAGCAAGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGa             >  2:919737/4‑149 (MQ=255)
                                                                                                                                         ccgtCTGCGCCAGCAAGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGa             >  4:768439/4‑149 (MQ=255)
                                                                                                                                            tCTGCGCCAGCAAGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGAcgc          <  3:768439/149‑1 (MQ=255)
                                                                                                                                            tCTGCGCCAGCAAGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGAcgc          >  4:876899/1‑149 (MQ=255)
                                                                                                                                               gcgcTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCgtg       >  2:1055855/1‑149 (MQ=255)
                                                                                                                                                   ttGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCa   >  2:380044/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCGCGTCGGACCATCAAGGCGGTAATGCTGCGGCTGTTTTTGCAGCAAACCAATCAACCAGGCCGGATTAACGTGATTCTTCTCGGCAAATTCGATCACCCCGCCTTTCTCATTACCTTCCAGCTTCCTGATCCCCAGTTTCTGCGCTTGCTGGCGCAGACGGGCAATATCCAGCAGGGTACGCGCCGGATCTGGCAGCAGGCCGAAGCGATCGATAAGCTCGACTTTGATCTCTTCCAGTTCGTTTTCCGTTTTGGCGCTGGCAATACGTTTGTAGAACGACAGACGCGTGTTCAC  >  NC_000913/1170611‑1170907

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: