Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,722,492 C→A T847K (ACG→AAG)  pka → protein lysine acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,722,4920CA98.9% 298.6 / NA 88T847K (ACG→AAG) pkaprotein lysine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (51/36);  total (51/37)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACC  >  NC_000913/2722347‑2722636
                                                                                                                                                 |                                                                                                                                                
gCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGAt                                                                                                                                               >  1:253008/1‑149 (MQ=255)
gCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGAt                                                                                                                                               >  6:338230/1‑149 (MQ=255)
 cGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGAt                                                                                                                                               <  7:193686/148‑1 (MQ=255)
 cGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGAt                                                                                                                                               <  5:117763/148‑1 (MQ=255)
 cGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATg                                                                                                                                              >  5:74563/1‑149 (MQ=255)
  gATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGACGACGCTTAATGGATAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGc                                                                                                                                             >  7:132992/1‑149 (MQ=255)
       ccGCTCCTGATAACATCGATGACGTATTTGCTGTACTGGTTCGCTCGGATGTCTAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCaaa                                                                                                                                         <  8:132992/148‑1 (MQ=255)
          aTCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAAt                                                                                                                                     <  3:316903/149‑1 (MQ=255)
               gATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTTAATGGTATTAAGATGCCaa                                                                                                                                          <  6:295382/139‑1 (MQ=255)
               gATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCaa                                                                                                                                          >  5:295382/1‑139 (MQ=255)
                aTAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTg                                                                                                                                 <  2:213610/147‑1 (MQ=255)
                  aaCATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCAt                                                                                                                             >  8:30352/1‑149 (MQ=255)
                  aaCATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCAt                                                                                                                             >  1:353435/1‑149 (MQ=255)
                  aaCATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCAt                                                                                                                             >  4:122838/1‑149 (MQ=255)
                   aCATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCAtg                                                                                                                            <  5:201359/149‑1 (MQ=255)
                    cATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCAtgg                                                                                                                           >  1:43634/1‑149 (MQ=255)
                        gATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCAtggtgg                                                                                                                        <  4:110158/148‑1 (MQ=255)
                         aTGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGgcg                                                                                                                      >  5:317743/1‑149 (MQ=255)
                              gAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGc                                                                                                                 >  3:221432/1‑149 (MQ=255)
                                     cTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCaa                                                                                                           <  3:116039/148‑1 (MQ=255)
                                        tACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTc                                                                                                       >  2:310128/1‑149 (MQ=255)
                                        tACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTc                                                                                                       >  7:104722/1‑149 (MQ=255)
                                           tGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCggg                                                                                                    >  5:234727/1‑149 (MQ=255)
                                             gTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATTGTTTTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGtt                                                                                                  >  5:64407/1‑149 (MQ=255)
                                             gTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGtt                                                                                                  >  3:41000/1‑149 (MQ=255)
                                             gTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGt                                                                                                   <  1:218116/148‑1 (MQ=255)
                                                   tCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCCTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGt                                                                                            >  8:185069/1‑149 (MQ=255)
                                                     ggATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGtt                                                                                           >  4:227200/1‑148 (MQ=255)
                                                     ggATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGtt                                                                                           <  4:41000/148‑1 (MQ=255)
                                                     ggATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTg                                                                                          <  6:74563/149‑1 (MQ=255)
                                                     ggACCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGtt                                                                                           <  8:211623/148‑1 (MQ=255)
                                                            aaaGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATcc                                                                                    <  4:129094/148‑1 (MQ=255)
                                                            aaaGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATc                                                                                     >  3:79879/1‑147 (MQ=255)
                                                             aaGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAg                                                                                  >  8:225932/1‑149 (MQ=255)
                                                             aaGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAg                                                                                  <  2:353435/149‑1 (MQ=255)
                                                               gggTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACCACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCt                                                                                >  1:283287/1‑149 (MQ=255)
                                                               gggTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCt                                                                                <  4:221432/149‑1 (MQ=255)
                                                                 gTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTc                                                                               >  5:156310/1‑148 (MQ=255)
                                                                  ttAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGa                                                                             <  2:17018/149‑1 (MQ=255)
                                                                      gCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAagag                                                                         >  5:5848/1‑149 (MQ=255)
                                                                      gCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAagag                                                                         >  7:315320/1‑149 (MQ=255)
                                                                       cTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCTCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAgg                                                                        <  7:30352/149‑1 (MQ=255)
                                                                        ttAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCa                                                                                   <  4:249411/137‑1 (MQ=255)
                                                                        ttAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCa                                                                                   >  3:249411/1‑137 (MQ=255)
                                                                        ttAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAggg                                                                       >  3:136381/1‑149 (MQ=255)
                                                                        ttAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAggg                                                                       >  2:163347/1‑149 (MQ=255)
                                                                          aGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAggg                                                                       <  3:122838/147‑1 (MQ=255)
                                                                          aGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAggg                                                                       <  1:310128/147‑1 (MQ=255)
                                                                               gACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGtt                                                                >  7:261930/1‑149 (MQ=255)
                                                                                  gCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTggg                                                             <  5:338230/149‑1 (MQ=255)
                                                                                    ttAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAgg                                                                        <  1:292531/136‑1 (MQ=255)
                                                                                    ttAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAgg                                                                        >  2:292531/1‑136 (MQ=255)
                                                                                        tGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACg                                                       >  2:14211/1‑149 (MQ=255)
                                                                                        tGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACg                                                       >  2:341986/1‑149 (MQ=255)
                                                                                            aaaGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGAGGCGCTAGCCCGCAAGCTCGGGTTTAACGATGATATCCAGCACGAAGAGGGGATCGTTGGGCTTACGcaa                                                    >  4:109911/1‑146 (MQ=255)
                                                                                             aaGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTaaa                                                  >  8:29499/1‑149 (MQ=255)
                                                                                               gTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATc                                                <  1:264219/149‑1 (MQ=255)
                                                                                                   aTTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGc                                            >  7:342877/1‑149 (MQ=255)
                                                                                                    ttACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCATCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGcc                                           >  3:45811/1‑149 (MQ=255)
                                                                                                           taCGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGGGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAgc                                       >  3:320757/1‑146 (MQ=255)
                                                                                                           taCGCGAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGGGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAgc                                       <  4:320757/146‑1 (MQ=255)
                                                                                                             cgcgAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAACCAATCGTGGCATGGGGGCGCTAGCCCGCAAGCGCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAg                                  <  2:253008/149‑1 (MQ=255)
                                                                                                             cgcgAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGa                                   <  6:156310/148‑1 (MQ=255)
                                                                                                             cgcgAGATCACGGACTACAACGGCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGa                                   <  6:64407/148‑1 (MQ=255)
                                                                                                              gcgAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAgc                                       <  3:45217/143‑1 (MQ=255)
                                                                                                              gcgAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAgc                                       >  4:45217/1‑143 (MQ=255)
                                                                                                               cgAGATCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAgg                                 >  2:264886/1‑148 (MQ=255)
                                                                                                                gagaTCACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGaa                               >  4:129680/1‑149 (MQ=255)
                                                                                                                     cACGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCAt                           >  4:25633/1‑148 (MQ=255)
                                                                                                                      aCGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCAt                           >  3:331865/1‑147 (MQ=255)
                                                                                                                      aCGGACTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTAGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGa                         >  6:353252/1‑149 (MQ=255)
                                                                                                                          aCTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGt                       >  4:140689/1‑147 (MQ=255)
                                                                                                                           cTACAACGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAg                    <  8:315320/149‑1 (MQ=255)
                                                                                                                               aaCGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTAc                >  4:157750/1‑149 (MQ=255)
                                                                                                                               aaCGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGAGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTAc                >  1:221995/1‑149 (MQ=255)
                                                                                                                                 cGTCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTg              <  6:5848/149‑1 (MQ=255)
                                                                                                                                   tCTGAATGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTgg             >  8:335179/1‑148 (MQ=255)
                                                                                                                                        aTGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATg        <  7:29499/148‑1 (MQ=255)
                                                                                                                                        aTGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGt       >  2:315755/1‑149 (MQ=255)
                                                                                                                                        aTGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGt       <  8:342877/149‑1 (MQ=255)
                                                                                                                                         tGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTCACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGtt      <  6:317743/149‑1 (MQ=255)
                                                                                                                                         tGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGtt      <  3:25633/149‑1 (MQ=255)
                                                                                                                                         tGGTATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGtt      >  2:77115/1‑149 (MQ=255)
                                                                                                                                            tATTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGGTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGAc   <  2:283287/149‑1 (MQ=255)
                                                                                                                                             aTTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAgg                   <  7:134377/132‑1 (MQ=255)
                                                                                                                                             aTTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAgg                   >  8:134377/1‑132 (MQ=255)
                                                                                                                                             aTTAAGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCGAAATCTTGCCCAGCGCGAGGAATCATGAGGAAGGTACTGGAAATGTTGAcc  >  8:141634/1‑149 (MQ=255)
                                                                                                                                             aTTAAGATGCCAAACAATCGTGGCAGGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGAACGTTGGGCTTACGCTAAATCTTGCCCAGCGCGTGGAATCATGAGTAAGGTACTGGAAATGTTGAc   >  1:30960/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCGATTTCCGATCCTGATAACATCGATGCCGAATTTGCTGTACTGGTTCGCTCGGATCTCAAAGGGTTAGGCTTAGGTCGACGCTTAATGGAAAAGTTGATTACCTATACGCGAGATCACGGACTACAACGTCTGAATGGTATTACGATGCCAAACAATCGTGGCATGGTGGCGCTAGCCCGCAAGCTCGGGTTTAACGTTGATATCCAGCTCGAAGAGGGGATCGTTGGGCTTACGCTAAATCTTGCCCAGCGCGAGGAATCATGAGTAAGGTACTGGAAATGTTGACC  >  NC_000913/2722347‑2722636

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