Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,720,210 G→A 8.6% L86L (CTG→CTA) ‡ pka → protein lysine acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,720,2100GA 8.6% 548.2 / 32.0 174L86L (CTG→CTA) ‡pkaprotein lysine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (79/80);  new base A (6/9);  total (85/89)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

Only 100 of 179 total aligned reads displayed.

CGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTACTGGCTCCCTGGCAA  >  NC_000913/2720062‑2720354
                                                                                                                                                    |                                                                                                                                                
cGGGAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTg                                                                                                                                                  <  3:934060/149‑1 (MQ=255)
   gAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGaa                                                                                                                                               <  3:695659/149‑1 (MQ=255)
   gAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGaa                                                                                                                                               <  4:690903/149‑1 (MQ=255)
   gAGGCTTTAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGaa                                                                                                                                               <  3:1234787/149‑1 (MQ=255)
        tttAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCt                                                                                                                                          >  3:198863/1‑149 (MQ=255)
        tttAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCt                                                                                                                                          >  1:276493/1‑149 (MQ=255)
        tttAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCCAGAAGAGca                                                                                                                                          >  3:313089/1‑148 (MQ=255)
        tttAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCCAGAAGAGca                                                                                                                                          <  3:1172405/149‑2 (MQ=255)
         ttAACGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCCAGAAGAGcaa                                                                                                                                         <  4:313089/149‑3 (MQ=255)
           aaCGGACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCgg                                                                                                                                       >  1:1427456/1‑149 (MQ=255)
              ggACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCga                                                                                                                                    <  1:152560/149‑1 (MQ=255)
              ggACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCga                                                                                                                                    >  3:842929/1‑149 (MQ=255)
              ggACCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCga                                                                                                                                    >  3:93438/1‑149 (MQ=255)
                aCCGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCgaga                                                                                                                                  <  1:1792643/149‑1 (MQ=255)
                  cGGTACTCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGGTAGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACTGTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCTGCGAGaaa                                                                                                                                >  4:1167203/1‑149 (MQ=255)
                    gTACTCCCGGTGACGCCAGCCTGGAAAGCGGAGTTGGGTGTGTTGGCCTGGCCGGGTATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAgg                                                                                                                              <  1:538952/149‑1 (MQ=255)
                        tCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTg                                                                                                                                                  <  2:11884/125‑1 (MQ=255)
                        tCCCGGTGACGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTg                                                                                                                                                  >  1:11884/1‑125 (MQ=255)
                                aCGCCAGCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTg                                                                                                                  <  3:359484/149‑1 (MQ=255)
                                      gCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCt                                                                                                                            <  1:408148/133‑1 (MQ=255)
                                      gCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCt                                                                                                                            >  2:408148/1‑133 (MQ=255)
                                      gCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCattat                                                                                                            <  1:188504/149‑1 (MQ=255)
                                      gCCTGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCattat                                                                                                            >  2:1588831/1‑149 (MQ=255)
                                         tGGAAAGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAaga                                                                                                                                             <  2:27919/113‑1 (MQ=255)
                                             aaGCGGTGTTGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTcc                                                                                                     >  3:996264/1‑149 (MQ=255)
                                                     ttGGGTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCtt                                                                                                        <  4:276570/138‑1 (MQ=255)
                                                       gggTGTGTTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCAt                                                                                           >  3:1477881/1‑149 (MQ=255)
                                                         gtgtgtTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGaa                                                                                                                                               >  1:147011/1‑95 (MQ=255)
                                                         gtgtgtTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGaaa                                                                                                                                <  3:1326191/110‑1 (MQ=255)
                                                         gtgtgtTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCCAGAAGAGcaa                                                                                                                                         >  3:47393/1‑99 (MQ=255)
                                                            tgtTGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCaa                                                                                      <  4:842929/149‑1 (MQ=255)
                                                              ttGGCCTGGCCGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAaca                                                                                    <  1:1270391/149‑1 (MQ=255)
                                                                 gcctggccGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGa                                                                                 >  4:1054731/1‑149 (MQ=255)
                                                                    tggccGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCttt                                                                                                       <  1:109072/124‑1 (MQ=255)
                                                                    tggccGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCCAGAAGAGcaa                                                                                                                                         >  3:264287/1‑88 (MQ=255)
                                                                        cGGATATTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGAtctc                                                                          <  3:577293/149‑1 (MQ=255)
                                                                            tatTGCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCcgcg                                                                      >  2:1014810/1‑149 (MQ=255)
                                                                                gCCAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGcctg                                                                  >  2:1077360/1‑149 (MQ=255)
                                                                                 ccAGCTTGCCCTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTc                                                                                                      <  3:921439/112‑1 (MQ=255)
                                                                                          ccTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTc                                                                                                      >  4:1516775/1‑103 (MQ=255)
                                                                                          ccTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCg                                                                                         <  4:1403619/116‑1 (MQ=255)
                                                                                          ccTTTACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGcct                                                                   >  2:560514/1‑138 (MQ=255)
                                                                                            tttACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGc                                                                                                                                      <  1:83355/69‑1 (MQ=255)
                                                                                            tttACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGc                                                                                                                                      >  2:83355/1‑69 (MQ=255)
                                                                                            tttACACCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCa                                                      >  2:1399623/1‑149 (MQ=255)
                                                                                                 acCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCaa                                                                                      <  2:1706231/112‑1 (MQ=255)
                                                                                                 acCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCaa                                                                                      >  1:1706231/1‑112 (MQ=255)
                                                                                                 acCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGcctgc                                                                 <  1:82398/133‑1 (MQ=255)
                                                                                                 acCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGcctgc                                                                 >  2:82398/1‑133 (MQ=255)
                                                                                                 acCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACa                                                 >  4:240734/1‑149 (MQ=255)
                                                                                                 acCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCCAGAAGAGcaa                                                                                                                                         <  1:1378704/61‑3 (MQ=38)
                                                                                                 acCCGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCCAGAAGAGcaa                                                                                                                                         >  2:1378704/1‑59 (MQ=255)
                                                                                                  cccGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCttt                                                                                                       <  1:355889/94‑1 (MQ=255)
                                                                                                  cccGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACAt                                                >  4:510967/1‑149 (MQ=255)
                                                                                                   ccGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCg                                                                                         >  2:1437587/1‑107 (MQ=255)
                                                                                                   ccGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATg                                               >  2:1270073/1‑149 (MQ=255)
                                                                                                   ccGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATg                                               >  4:664103/1‑149 (MQ=255)
                                                                                                   ccGACCTTGCGGTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATg                                               >  2:67706/1‑149 (MQ=255)
                                                                                                              gTTTTATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTgg                                    >  2:1809930/1‑149 (MQ=255)
                                                                                                               ttttATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCcgc                                                                       >  3:1112038/1‑113 (MQ=255)
                                                                                                               ttttATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCcgc                                                                       <  4:1112038/113‑1 (MQ=255)
                                                                                                               ttttATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATaa                                                   >  3:1252396/1‑133 (MQ=255)
                                                                                                               ttttATGTACCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGa                                   <  2:267862/149‑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)
                                                                                                                       aCCAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCg                                                                                               >  3:1425996/1‑81 (MQ=255)
                                                                                                                        ccAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGt                          <  3:664103/149‑1 (MQ=255)
                                                                                                                        ccAATGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGt                          <  1:1816067/149‑1 (MQ=255)
                                                                                                                         cAATGCCAGCCGTAATCTTTCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAATACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGcctgcgccctg                                                           >  4:1306423/1‑115 (MQ=255)
                                                                                                                          aaTGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCt                        <  1:1014810/149‑1 (MQ=255)
                                                                                                                           ttGCCAGCCGTAATCTTGCTCTTCTGAAAGAGCTCGGCGAGAAAGGCTGCAACCCCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCCCCGGGACTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTg                       <  4:1472512/148‑1 (MQ=255)
                                                                                                                           aTGCCAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGa                                   <  4:402773/137‑1 (MQ=255)
                                                                                                                               cAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGaa                                                                                <  3:370471/88‑1 (MQ=255)
                                                                                                                               cAGCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGaa                                                                                >  4:370471/1‑88 (MQ=255)
                                                                                                                                 gCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTa                 <  1:1270073/149‑1 (MQ=255)
                                                                                                                                 gCCGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTa                 >  4:167148/1‑149 (MQ=255)
                                                                                                                                  ccGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAAc                             >  2:574237/1‑136 (MQ=255)
                                                                                                                                  ccGTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAAc                             <  1:574237/136‑1 (MQ=255)
                                                                                                                                    gTAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGc                                                                                                                             <  4:926793/38‑1 (MQ=255)
                                                                                                                                     tAATCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTACTgg             <  1:609610/149‑1 (MQ=255)
                                                                                                                                        tCTTGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTACTGGCTc          >  3:721004/1‑149 (MQ=255)
                                                                                                                                         cgagCTCTTCCAGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATgc                                              >  2:317301/4‑112 (MQ=255)
                                                                                                                                         cgagCTCTTCCAGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATgc                                              <  1:317301/109‑1 (MQ=255)
                                                                                                                                          gagCTCTTCCAGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGcct                                                                   >  2:191853/3‑90 (MQ=255)
                                                                                                                                           tGCTCTTCTGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTACTGGCTCCCt       <  3:1054731/149‑1 (MQ=255)
                                                                                                                                              tcttctGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTcc                                                                         >  2:11033/1‑80 (MQ=255)
                                                                                                                                              tcttctGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTgg                      >  1:886789/1‑131 (MQ=255)
                                                                                                                                              tcttctGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTgg                      <  2:886789/131‑1 (MQ=255)
                                                                                                                                              tcttctGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTgg                      <  1:1664580/131‑1 (MQ=255)
                                                                                                                                              tcttctGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTACTGGCTCCCTGGc    >  1:460967/1‑149 (MQ=255)
                                                                                                                                                ttctGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCaaa                              >  2:956711/1‑121 (MQ=255)
                                                                                                                                                ttctGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGtt                   <  1:24878/132‑1 (MQ=255)
                                                                                                                                                ttctGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTACTg              >  3:1590819/1‑137 (MQ=255)
                                                                                                                                                ttctGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTACTGGCTCCCTGGCaa  <  2:460967/149‑1 (MQ=255)
                                                                                                                                                ttctGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTACTGGCTCCCTGGCaa  >  4:1248670/1‑149 (MQ=255)
                                                                                                                                                ttctGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTACTGGCTCCCTGGCaa  <  1:1588831/149‑1 (MQ=255)
                                                                                                                                                  ctGGAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTACTGGCTCCCt       >  4:1445885/1‑142 (MQ=255)
                                                                                                                                                    ggAAGAGCTCGGCGAGAAAGGCTGTAAAACCTGCATTATTCTTTCCGCCCCGGCATCGCAACACGAAGATCTCCGCGCCTGCGCCCTGCGCCATAACATGCGCCTGCTTGGACCAAACAGTCTGGGTTTACTgg             <  2:15166/134‑1 (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 < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: