Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,282,330 C→T 5.1% E27K (GAG→AAG)  PP_2883 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,282,3300CT 5.1% 272.7 / 2.3 91E27K (GAG→AAG) PP_2883hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (48/38);  new base T (3/2);  total (51/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CTGCCAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGAAACTG  >  NC_002947/3282184‑3282477
                                                                                                                                                  |                                                                                                                                                    
cTGCCAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAg                                                                                                                                                   <  1:329941/150‑1 (MQ=255)
   ccAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAggtg                                                                                                                                                <  1:1790250/150‑1 (MQ=255)
    cAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAggtgg                                                                                                                                               <  1:1916420/150‑1 (MQ=255)
     aGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAggtgg‑t                                                                                                                                             >  1:1413142/1‑150 (MQ=255)
         ttCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCA‑CCACCTTGAGGTGGCTGCAc                                                                                                                                         <  1:1249037/150‑1 (MQ=255)
               gccggccTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAg                                                                                                                                   >  2:2046316/1‑150 (MQ=255)
                ccggccTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCGCCTCAAGGTGG‑TGCACGCCAAgg                                                                                                                                  <  2:1192757/150‑1 (MQ=255)
                  ggccTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGc                                                                                                                                <  2:1120252/150‑1 (MQ=255)
                  ggccTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGc                                                                                                                                <  1:772776/150‑1 (MQ=255)
                   gccTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGcc                                                                                                                               >  1:1666089/1‑150 (MQ=255)
                     cTGTTCGCCGTCGGGGGCTTGGCAACGATTGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCcac                                                                                                                             <  1:1077210/150‑1 (MQ=255)
                     cTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAggtg                                                                                                                                                >  1:69841/1‑132 (MQ=255)
                     cTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAggtg                                                                                                                                                <  2:69841/132‑1 (MQ=255)
                     cTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCcac                                                                                                                             >  2:2093487/1‑150 (MQ=255)
                             cGTCGGGGGCGTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCa                                                                                                                     <  2:2203859/150‑1 (MQ=255)
                              gTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAg                                                                                                                    <  2:1490736/150‑1 (MQ=255)
                              gTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAg                                                                                                                    <  1:1924808/150‑1 (MQ=255)
                                     gCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTgg                                                                                                             >  2:1177210/1‑150 (MQ=255)
                                       ttGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGAACAACAAGGGCCACATTGCGCAGGCACTGGTa                                                                                                           <  1:1881952/150‑1 (MQ=255)
                                            aacgaacgCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACAt                                                                                                      >  2:491323/1‑150 (MQ=255)
                                            aacgaacgCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACAt                                                                                                      <  1:914294/150‑1 (MQ=255)
                                               gaacgCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATccc                                                                                                   >  1:2323096/1‑150 (MQ=255)
                                                 acgCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCNCgg                                                                                                 <  2:1496884/150‑1 (MQ=255)
                                                           ccAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCCGGGATTATGCTCAAGGTGG‑TGCACGCCAAGGGCCACCTTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTg                                                                                       >  1:2253630/1‑150 (MQ=255)
                                                                  ttCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCt                                                                                >  2:195948/1‑150 (MQ=255)
                                                                  ttCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCt                                                                                >  1:2176457/1‑150 (MQ=255)
                                                                   tCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCtg                                                                               <  1:195948/150‑1 (MQ=255)
                                                                     ggTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCtgtg                                                                             >  2:286861/1‑150 (MQ=255)
                                                                       tATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTggg                                                                           >  2:45928/1‑150 (MQ=255)
                                                                         tCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTc                                                                         >  1:574721/1‑150 (MQ=255)
                                                                           cAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTa                                                                                                           >  2:1982870/1‑114 (MQ=255)
                                                                           cAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTa                                                                                                           <  1:1982870/114‑1 (MQ=255)
                                                                           cAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCaa                                                                       >  1:2312970/1‑150 (MQ=255)
                                                                            aGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAg                                                                      >  2:613865/1‑150 (MQ=255)
                                                                             gggCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCaata                                                                     <  2:2312970/150‑3 (MQ=255)
                                                                               gCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATc                                                                   >  1:670912/1‑150 (MQ=255)
                                                                                ccAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCCCTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCa                                                                  >  2:807881/1‑150 (MQ=255)
                                                                                ccAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCa                                                                  <  2:2323096/150‑1 (MQ=255)
                                                                                 cAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAggg                                                                                                                                 >  2:1307658/1‑86 (MQ=255)
                                                                                 cAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAggg                                                                                                                                 <  1:1307658/86‑1 (MQ=255)
                                                                                      tCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGTAGTCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAAtt                                                            >  1:61373/1‑150 (MQ=255)
                                                                                       cAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTc                                                           <  1:1160751/150‑1 (MQ=255)
                                                                                       cAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTc                                                           <  2:1666089/150‑1 (MQ=255)
                                                                                         gTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGc                                                         <  1:2046316/150‑1 (MQ=255)
                                                                                            gCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCgg                                                      >  1:85148/1‑150 (MQ=255)
                                                                                             ccAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGc                                                     >  2:1289361/1‑150 (MQ=255)
                                                                                              cAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCg                                                    >  1:40305/1‑150 (MQ=255)
                                                                                               aGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCgt                                                   <  2:574721/150‑1 (MQ=255)
                                                                                               aGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCgt                                                   >  2:1856482/1‑150 (MQ=255)
                                                                                                gCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCgtg                                                  >  1:124799/1‑150 (MQ=255)
                                                                                                 cccTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTgg                                                                                                             >  1:675519/1‑90 (MQ=255)
                                                                                                 cccTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTgg                                                                                                             <  2:675519/90‑1 (MQ=255)
                                                                                                    tGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTgg                                              >  1:1110580/1‑150 (MQ=255)
                                                                                                         gCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCtg                                         <  2:124799/150‑1 (MQ=255)
                                                                                                           cGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCg                                                                                                  <  2:390768/91‑1 (MQ=255)
                                                                                                           cGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCg                                                                                                  >  1:390768/1‑91 (MQ=255)
                                                                                                           cGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCtgat                                       <  1:1565811/150‑1 (MQ=255)
                                                                                                           cGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCtgat                                       >  1:1165997/1‑150 (MQ=255)
                                                                                                              cATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACTCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTgatgat                                    >  2:1350849/1‑150 (MQ=255)
                                                                                                                tCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGAtct                                  >  1:1690122/1‑150 (MQ=255)
                                                                                                                  gATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCtga                                        <  2:2297494/142‑1 (MQ=255)
                                                                                                                  gATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCtga                                        <  2:2297151/142‑1 (MQ=255)
                                                                                                                  gATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCtga                                        >  1:2297151/1‑142 (MQ=255)
                                                                                                                  gATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCtga                                        >  1:2297494/1‑142 (MQ=255)
                                                                                                                  gATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCg                                <  1:579724/150‑1 (MQ=255)
                                                                                                                  gATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCg                                >  2:569065/1‑150 (MQ=255)
                                                                                                                   tgcacgccgcgATGGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCgg                               >  1:2108950/6‑150 (MQ=255)
                                                                                                                    tCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGc                              <  2:2253630/150‑1 (MQ=255)
                                                                                                                      aGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCt                            >  2:2183502/1‑150 (MQ=255)
                                                                                                                       ggCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTc                           >  1:196930/1‑150 (MQ=255)
                                                                                                                        gccgccATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCCAGGGCCACATTGCGGAGGCACTGGTAGACATCCCGGGCAATCGGGGAACTGCTGTGGGGCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCa                          >  2:865350/1‑150 (MQ=255)
                                                                                                                         ccgccATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTAAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTGAt                         >  1:1111355/1‑150 (MQ=255)
                                                                                                                             gcccttggcgTGCA‑CCACCTTGAGGTGGCTGCACGCCGAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCgg                     >  1:2394942/7‑150 (MQ=38)
                                                                                                                               tggcggcgTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGc                   <  2:2108950/148‑1 (MQ=255)
                                                                                                                               ccttggcgTGCA‑CCACCTTGAGGTGGCTGCACGCCGAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAg                                                                 >  2:957652/5‑104 (MQ=25)
                                                                                                                               ccttggcgTGCA‑CCACCTTGAGGTGGCTGCACGCCGAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAg                                                                 <  1:957652/100‑1 (MQ=25)
                                                                                                                               ccttggcgTGCA‑CCACCTTGAGGTGGCTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGc                   >  2:2261466/5‑150 (MQ=255)
                                                                                                                                cgcgGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCGCATCCGGGcc                  <  2:1165997/150‑1 (MQ=255)
                                                                                                                                  cggcgTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGa                <  2:670912/150‑1 (MQ=255)
                                                                                                                                   ggcgTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGAt               >  1:1407804/1‑150 (MQ=255)
                                                                                                                                        gCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGAt               <  2:1927101/145‑1 (MQ=255)
                                                                                                                                        gCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGAt               >  1:1927101/1‑145 (MQ=255)
                                                                                                                                         cAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCtt         <  1:613865/150‑1 (MQ=255)
                                                                                                                                              aCCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGc                   >  2:294760/1‑135 (MQ=255)
                                                                                                                                              aCCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGc                   <  1:294760/135‑1 (MQ=255)
                                                                                                                                              aCCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGAAAc    >  2:2220886/1‑150 (MQ=255)
                                                                                                                                              aCCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGAAAc    >  1:1846143/1‑150 (MQ=255)
                                                                                                                                                cTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTc                                                                         <  1:1110620/79‑1 (MQ=255)
                                                                                                                                                cTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTc                                                                         >  2:1110620/1‑79 (MQ=255)
                                                                                                                                                cTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGAAACtg  >  1:1025778/1‑150 (MQ=255)
                                                                                                                                                cTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGAAACtg  >  1:479908/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTGCCAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGAAACTG  >  NC_002947/3282184‑3282477

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