Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,526,411 A→C 5.9% K219Q (AAG→CAG) ‡ PP_2215 → acetyl‑CoA acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,526,4110AC 5.9% 216.2 / 3.7 68K219Q (AAG→CAG) ‡PP_2215acetyl‑CoA acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (34/30);  new base C (2/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

CCCAGGCCAATGCCTTCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATGCGCCGCTCCGAG  >  NC_002947/2526265‑2526557
                                                                                                                                                  |                                                                                                                                                  
cccAGGCCAATGCCTTCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGGCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCaa                                                                                                                                                   <  2:282687/148‑1 (MQ=255)
 ccAGGCCAATGCCTTCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAg                                                                                                                                                  >  3:321203/1‑148 (MQ=255)
 ccAGGCCAATGCCTTCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAg                                                                                                                                                  >  3:51300/1‑148 (MQ=255)
   aGGCCAATGCCTTCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGAt                                                                                                                                               >  7:36362/1‑149 (MQ=255)
     gCCAATGCCTTCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGa                                                                                                                                             <  2:376950/149‑1 (MQ=255)
     gCCAATGCCTTCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGa                                                                                                                                             >  2:157736/1‑149 (MQ=255)
      ccAATGCCTTCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATGAAGGTTGag                                                                                                                                            <  2:341198/149‑2 (MQ=255)
       cAATGCCTTCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGAc                                                                                                                                            >  6:301481/1‑148 (MQ=255)
       cAATGCCTTCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGACGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGAc                                                                                                                                            >  5:196338/1‑148 (MQ=255)
             cTTCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGagcagc                                                                                                                                     >  1:13303/1‑149 (MQ=255)
              ttCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGGTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAgcc                                                                                                                                    >  2:224932/1‑149 (MQ=255)
              ttCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAgcc                                                                                                                                    >  4:250106/1‑149 (MQ=255)
                 aGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAgccgcc                                                                                                                                 >  2:204722/1‑149 (MQ=255)
                 aGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAgccgc                                                                                                                                  <  8:2559/148‑1 (MQ=255)
                 aGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAgccgc                                                                                                                                  <  4:216923/148‑1 (MQ=255)
                  gCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAgccgcc                                                                                                                                 >  6:44094/1‑148 (MQ=255)
                    cgcgAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGccc                                                                                                                                >  1:353694/1‑147 (MQ=255)
                      cgAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAgg                                                                                                                            >  5:130910/1‑149 (MQ=255)
                         ggCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGgcgc                                                                                                                         <  4:359428/149‑1 (MQ=255)
                                gACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGa                                                                                                                  <  1:157736/149‑1 (MQ=255)
                                 aCCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGAc                                                                                                                 <  3:82131/149‑1 (MQ=255)
                                  ccAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCCTTGATGACgcg                                                                                                                                         <  5:372014/124‑3 (MQ=255)
                                  ccAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCCTTGATGACgcg                                                                                                                                         >  6:372014/1‑122 (MQ=255)
                                  ccAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGAc                                                                                                                 >  4:68806/1‑148 (MQ=255)
                                     gTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGa                                                                                                             <  2:353694/149‑1 (MQ=255)
                                     gTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGACTGGAc                                                                                                                 >  3:75018/1‑145 (MQ=255)
                                         gCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTg                                                                                                          <  4:51300/148‑1 (MQ=255)
                                           cATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTgcg                                                                                                        >  2:200479/1‑148 (MQ=255)
                                           cATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTgcg                                                                                                        >  7:117264/1‑148 (MQ=255)
                                           cATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTgcg                                                                                                        >  7:182192/1‑148 (MQ=255)
                                               gCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGc                                                                                                    <  8:182192/148‑1 (MQ=255)
                                               gCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCt                                                                                                   <  7:162913/149‑1 (MQ=255)
                                                    gCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCaaa                                                                                               <  6:148309/148‑1 (MQ=255)
                                                     cTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAAcc                                                                                             >  3:329211/1‑149 (MQ=255)
                                                        aCCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGTATGACGAGCAGCCGCCCAAGTCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGc                                                                                          <  3:95874/149‑1 (MQ=255)
                                                          ccGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCa                                                                                                      >  6:37604/1‑135 (MQ=255)
                                                          ccGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCa                                                                                                      <  5:37604/135‑1 (MQ=255)
                                                           cGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGtt                                                                                       <  4:3898/149‑1 (MQ=255)
                                                           cGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGt                                                                                        >  5:155298/1‑148 (MQ=255)
                                                           cGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGt                                                                                        <  6:155298/148‑1 (MQ=255)
                                                              gcgcAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCg                                                                                    <  8:36362/149‑1 (MQ=255)
                                                              gcgcAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCg                                                                                    >  3:60031/1‑149 (MQ=255)
                                                               cgcAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGt                                                                                   <  4:321203/149‑1 (MQ=255)
                                                                 cAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTg                                                                                  >  3:319124/1‑148 (MQ=255)
                                                                 cAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGa                                                                                 <  3:68806/149‑1 (MQ=255)
                                                                       gCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAggcgg                                                                           <  4:329211/149‑1 (MQ=255)
                                                                          aTCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGt                                                                                   <  3:122781/138‑1 (MQ=255)
                                                                          aTCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGt                                                                                   >  4:122781/1‑138 (MQ=255)
                                                                           tCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCAcc                                                                       <  4:319124/149‑1 (MQ=255)
                                                                           tCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCAcc                                                                       >  4:24729/1‑149 (MQ=255)
                                                                             agcagcGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGt                                                                     <  8:117264/149‑1 (MQ=255)
                                                                                             ccgccgAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGc                                                      <  2:13303/148‑1 (MQ=255)
                                                                                             ccgccgAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGc                                                      >  3:209141/1‑148 (MQ=255)
                                                                                                     aTCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGAt                                             <  5:301481/149‑1 (MQ=255)
                                                                                                      tCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGAtt                                            <  5:93346/149‑1 (MQ=255)
                                                                                                         tGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTcc                                         <  6:196338/149‑1 (MQ=255)
                                                                                                              gTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGAcgg                                    >  3:259064/1‑149 (MQ=255)
                                                                                                                  aaGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGAcggcgct                                >  4:222683/1‑149 (MQ=255)
                                                                                                                    gTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGAcggcgctgc                              >  3:119446/1‑149 (MQ=255)
                                                                                                                       aCCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGAcggcgctgcggc                           >  6:384311/1‑149 (MQ=255)
                                                                                                                        ccGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGAcgccgctgcggcg                          >  7:158458/1‑149 (MQ=255)
                                                                                                                         cGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGAcggcgctgcggcgc                         <  6:130910/149‑1 (MQ=255)
                                                                                                                           aGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGAcggcgctgcggcgctg                       <  3:250106/149‑1 (MQ=255)
                                                                                                                               cAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGAcggcgctgcggcgctg                       <  2:29062/145‑1 (MQ=255)
                                                                                                                               cAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGAcggcgctgcggcgctg                       >  1:29062/1‑145 (MQ=255)
                                                                                                                                aaGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGAc                                                                  >  6:12164/1‑101 (MQ=255)
                                                                                                                                     aaaGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATgc             >  3:335493/1‑149 (MQ=255)
                                                                                                                                         ctcgTCATCCTTGATGACGCGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGAtt                                            >  7:108865/3‑114 (MQ=255)
                                                                                                                                         ctcgTCATCCTTGATGACGCGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGAtt                                            <  8:108865/112‑1 (MQ=255)
                                                                                                                                             tcatcaAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATGCGCCGCTcc     <  5:44094/149‑1 (MQ=255)
                                                                                                                                              catcaAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATGCGCCGCTCCg    <  7:285413/149‑1 (MQ=255)
                                                                                                                                               atcaAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATGCGCCGCTcc     >  4:229500/1‑147 (MQ=255)
                                                                                                                                                tcaAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATGCGCCGCTCCGAg  >  3:12034/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CCCAGGCCAATGCCTTCAGCCGCGAGGCCCAGGACCAGTTCGCCATCGCCTCGCTGACCCGAGCGCAGGAAGCCATCAGCAGCGGCCGTTTTGCCGCCGAGATCGTGCCGGTGGAAGTCACCGAGGGCAAGGAAAAGCGCGTCATCAAGGATGACGAGCAGCCGCCCAAGGCGCGTCTGGACAAGATTGCGCAGCTCAAACCGGCGTTTCGTGAAGGCGGCACCGTGACGGCGGCCAACGCCAGTTCGATTTCCGACGGCGCTGCGGCGCTGGTACTGATGCGCCGCTCCGAG  >  NC_002947/2526265‑2526557

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: