Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,720,209 T→C 9.5% L86P (CTG→CCG) ‡ pka → protein lysine acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,720,2090TC 9.5% 491.0 / 32.6 157L86P (CTG→CCG) ‡pkaprotein lysine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (69/73);  new base C (6/9);  total (75/82)
Fisher's exact test for biased strand distribution p-value = 5.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01

Only 100 of 177 total aligned reads displayed.

CGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTACTGGCTCCCTGGCAA  >  NC_000913/2720062‑2720354
                                                                                                                                                   |                                                                                                                                                 
cGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTg                                                                                                                                                  <  3:934060/149‑1 (MQ=255)
   gAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGaa                                                                                                                                               <  4:690903/149‑1 (MQ=255)
   gAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGaa                                                                                                                                               <  3:1234787/149‑1 (MQ=255)
        tttAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCCAGAAGAGca                                                                                                                                          >  3:313089/1‑148 (MQ=255)
         ttAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCCAGAAGAGcaa                                                                                                                                         <  3:517525/149‑3 (MQ=255)
         ttAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCCAGAAGAGcaa                                                                                                                                         <  1:1616447/149‑3 (MQ=255)
          tAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCg                                                                                                                                        >  2:152560/1‑149 (MQ=255)
           aaCGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCgg                                                                                                                                       >  1:1427456/1‑149 (MQ=255)
              ggACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCga                                                                                                                                    >  3:842929/1‑149 (MQ=255)
              ggACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCga                                                                                                                                    <  1:152560/149‑1 (MQ=255)
                    gTACTCCCGGTGACGCCAGCCTGGAAAGCGGAGTTGGGTGTGTTGGCCTGGCCGGGTATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAgg                                                                                                                              <  1:538952/149‑1 (MQ=255)
                     tACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGc                                                                                                                             <  2:276493/149‑1 (MQ=255)
                        tCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTg                                                                                                                                                  <  2:11884/125‑1 (MQ=255)
                        tCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTg                                                                                                                                                  >  1:11884/1‑125 (MQ=255)
                          ccGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTaa                                                                                                                        >  1:832553/1‑149 (MQ=255)
                               gACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCt                                                                                                                   <  3:1195201/149‑1 (MQ=255)
                                       ccTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCattatt                                                                                                           <  2:1427456/149‑1 (MQ=255)
                                         tGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAaga                                                                                                                                             >  1:27919/1‑113 (MQ=255)
                                         tGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAAc                                                                                                                     >  1:85752/1‑137 (MQ=255)
                                         tGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAAc                                                                                                                     <  2:85752/137‑1 (MQ=255)
                                         tGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCt                                                                                                         >  3:365297/1‑149 (MQ=255)
                                          ggAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCGAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCtt                                                                                                        <  1:1780728/149‑1 (MQ=255)
                                            aaaGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTc                                                                                                      <  2:832553/149‑1 (MQ=255)
                                            aaaGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTc                                                                                                      >  2:1116373/1‑149 (MQ=255)
                                             aaGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTcc                                                                                                     >  2:535844/1‑149 (MQ=255)
                                                     ttGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCtt                                                                                                        >  3:276570/1‑138 (MQ=255)
                                                     ttGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCtt                                                                                                        <  4:276570/138‑1 (MQ=255)
                                                         gtgtgtTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGaa                                                                                                                                               >  1:147011/1‑95 (MQ=255)
                                                            tgtTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCaa                                                                                      <  4:842929/149‑1 (MQ=255)
                                                              ttGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAaca                                                                                    <  1:1270391/149‑1 (MQ=255)
                                                                 gcctggccGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGa                                                                                 >  4:1054731/1‑149 (MQ=255)
                                                                  cctggccGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGGATTATTCTTTCCGCCCCGGCATCGCAACACGaa                                                                                <  1:535844/149‑1 (MQ=255)
                                                                    tggccGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCttt                                                                                                       <  1:109072/124‑1 (MQ=255)
                                                                    tggccGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCttt                                                                                                       >  2:109072/1‑124 (MQ=255)
                                                                        cGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGAtctc                                                                          <  3:577293/149‑1 (MQ=255)
                                                                         ggATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAACCTTGCTCTTCTGGAAGAGCACGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACATGAAGATCTcc                                                                         <  3:1167203/149‑1 (MQ=255)
                                                                            tatTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCcgcg                                                                      >  2:1014810/1‑149 (MQ=255)
                                                                                 ccAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTc                                                                                                      <  3:921439/112‑1 (MQ=255)
                                                                                       tGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGcctgcgccctg                                                           <  4:1477881/149‑1 (MQ=255)
                                                                                          ccTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTc                                                                                                      <  3:1516775/103‑1 (MQ=255)
                                                                                          ccTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTc                                                                                                      >  4:1516775/1‑103 (MQ=255)
                                                                                          ccTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCg                                                                                         <  4:1403619/116‑1 (MQ=255)
                                                                                          ccTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCg                                                                                         >  3:1403619/1‑116 (MQ=255)
                                                                                          ccTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGcctgcgccctgcgc                                                        >  2:609610/1‑149 (MQ=255)
                                                                                           cTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCctgcgccctgcgcc                                                       >  2:439484/1‑149 (MQ=255)
                                                                                                 acCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACa                                                 >  4:1040269/1‑149 (MQ=255)
                                                                                                 acCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACa                                                 >  3:990034/1‑149 (MQ=255)
                                                                                                 acCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACa                                                 >  4:240734/1‑149 (MQ=255)
                                                                                                 acCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCCAGAAGAGcaa                                                                                                                                         >  2:1378704/1‑59 (MQ=255)
                                                                                                 acCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCCAGAAGAGcaa                                                                                                                                         <  1:1378704/61‑3 (MQ=38)
                                                                                                  cccGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCttt                                                                                                       >  2:355889/1‑94 (MQ=255)
                                                                                                  cccGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCttt                                                                                                       <  1:355889/94‑1 (MQ=255)
                                                                                                  cccGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACAt                                                >  4:510967/1‑149 (MQ=255)
                                                                                                   ccGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCg                                                                                         >  2:1437587/1‑107 (MQ=255)
                                                                                                   ccGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATg                                               >  4:1235768/1‑149 (MQ=255)
                                                                                                   ccGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATg                                               >  4:664103/1‑149 (MQ=255)
                                                                                                   ccGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATg                                               >  2:1270073/1‑149 (MQ=255)
                                                                                                    cGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATgc                                              <  1:439484/149‑1 (MQ=255)
                                                                                                         ttGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTg                                         >  1:1690926/1‑149 (MQ=255)
                                                                                                         ttGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTg                                         >  1:739031/1‑149 (MQ=255)
                                                                                                          tGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGc                                        <  3:240734/149‑1 (MQ=255)
                                                                                                             ggTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTg                                     >  1:267862/1‑149 (MQ=255)
                                                                                                             ggTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTg                                     >  3:70164/1‑149 (MQ=255)
                                                                                                              gTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTgg                                    >  2:1809930/1‑149 (MQ=255)
                                                                                                               ttttATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCcgc                                                                       <  4:1112038/113‑1 (MQ=255)
                                                                                                               ttttATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATaa                                                   <  4:1252396/133‑1 (MQ=255)
                                                                                                                tttATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGAc                                  <  1:1077360/149‑1 (MQ=255)
                                                                                                                  tATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCa                                >  3:1472512/1‑149 (MQ=255)
                                                                                                                  tATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCa                                >  2:285712/1‑149 (MQ=255)
                                                                                                                    tGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCaaa                              <  4:990034/149‑1 (MQ=255)
                                                                                                                    tGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCaaa                              <  3:1040269/149‑1 (MQ=255)
                                                                                                                      tACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACa                            >  3:775881/1‑149 (MQ=255)
                                                                                                                        ccAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGt                          <  1:1816067/149‑1 (MQ=255)
                                                                                                                        ccAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGt                          <  3:664103/149‑1 (MQ=255)
                                                                                                                        ccAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGt                          <  2:739031/149‑1 (MQ=255)
                                                                                                                         cAATGCCAGCCGTAATCTTTCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAATACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGcctgcgccctg                                                           >  4:1306423/1‑115 (MQ=255)
                                                                                                                         cAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGcctgcgccctg                                                           <  3:1306423/115‑1 (MQ=255)
                                                                                                                          aaTGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCt                        <  1:1014810/149‑1 (MQ=255)
                                                                                                                           aTGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGa                                   >  3:402773/1‑137 (MQ=255)
                                                                                                                           aTGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGa                                   <  4:402773/137‑1 (MQ=255)
                                                                                                                            tGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTgg                      <  3:106358/149‑1 (MQ=255)
                                                                                                                            tGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTgg                      <  1:1809930/149‑1 (MQ=255)
                                                                                                                               cAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGaa                                                                                >  4:370471/1‑88 (MQ=255)
                                                                                                                                 gCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTa                 >  4:167148/1‑149 (MQ=255)
                                                                                                                                 gCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTa                 <  4:365297/149‑1 (MQ=255)
                                                                                                                                 gCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTa                 <  1:1270073/149‑1 (MQ=255)
                                                                                                                                  ccGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAAc                             <  1:574237/136‑1 (MQ=255)
                                                                                                                                  ccGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAAc                             >  2:574237/1‑136 (MQ=255)
                                                                                                                                    gTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGc                                                                                                                             <  4:926793/38‑1 (MQ=255)
                                                                                                                                         cgagCTCTTCCAGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATgc                                              <  1:317301/109‑1 (MQ=255)
                                                                                                                                         cgagCTCTTCCAGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATgc                                              >  2:317301/4‑112 (MQ=255)
                                                                                                                                          gagCTCTTCCAGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGcct                                                                   >  2:191853/3‑90 (MQ=255)
                                                                                                                                          gagCTCTTCCAGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGcct                                                                   <  1:191853/88‑1 (MQ=255)
                                                                                                                                              tcttctGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTgg                      >  2:1664580/1‑131 (MQ=255)
                                                                                                                                              tcttctGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTgg                      >  1:886789/1‑131 (MQ=255)
                                                                                                                                              tcttctGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTACTGGCTCCCTGGc    >  3:1722763/1‑149 (MQ=255)
                                                                                                                                                ttctGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCaaa                              >  2:956711/1‑121 (MQ=255)
                                                                                                                                                ttctGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGtt                   >  2:24878/1‑132 (MQ=255)
                                                                                                                                                  ctGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTACTGGCTCCCt       <  3:1445885/142‑1 (MQ=255)
                                                                                                                                                  ctGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTACTGGCTCCCt       >  4:1445885/1‑142 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTACTGGCTCCCTGGCAA  >  NC_000913/2720062‑2720354

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