Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,642,656 C→T 7.2% Q185* (CAG→TAG) ‡ PP_4107 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,642,6560CT 7.2% 201.6 / 4.1 69Q185* (CAG→TAG) ‡PP_4107LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base T (2/3);  total (34/35)
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.64e-01

CATTCGACCTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCCGCGCGGCCATCTC  >  NC_002947/4642507‑4642804
                                                                                                                                                      |                                                                                                                                                    
cATTCGACCTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGccc                                                                                                                                                      >  2:1201918/1‑150 (MQ=255)
 aTTCGACCTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCa                                                                                                                                                     >  1:1178574/1‑150 (MQ=255)
       ccTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCGCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGct                                                                                                                                               >  1:549237/1‑150 (MQ=255)
          ggTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCTGGGCGctgta                                                                                                                                            <  2:135465/150‑3 (MQ=255)
              cAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCTGGGCGctgtcgaaa                                                                                                                                        <  1:1014100/150‑7 (MQ=255)
              cAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAg                                                                                                                                        <  2:858861/150‑1 (MQ=255)
                    aaGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGc                                                                                                                                  >  1:18972/1‑150 (MQ=255)
                         ggACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGtt                                                                                                                             <  2:738145/150‑1 (MQ=255)
                             gTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCaa                                                                                                                         >  1:1034499/1‑150 (MQ=255)
                             gTCGCCTTTTTCTACCGCAAGCGAACCTGGCCCCGGTCGACCTGCATCGAGTTGTTCAGAGAAGTAGGACGATGTGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCaa                                                                                                                         <  1:542639/150‑1 (MQ=255)
                                gCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGctga                                                                                                                                             <  2:954355/127‑1 (MQ=255)
                                gCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGctga                                                                                                                                             >  1:954355/1‑127 (MQ=255)
                                    ttttCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCgtgt                                                                                                                  <  2:142986/150‑1 (MQ=255)
                                     tttCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCaa                                                                                                                         >  2:162236/1‑142 (MQ=255)
                                     tttCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCaa                                                                                                                         <  1:162236/142‑1 (MQ=255)
                                                aGGGCACCTGGCCCGGGGCGGCCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAgg                                                                                                      <  1:1261854/150‑1 (MQ=255)
                                                 gggCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGc                                                                                                     >  1:404421/1‑150 (MQ=255)
                                                     aCCTGGCCCGGGGCGATCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTgg                                                                                                 >  1:70242/1‑150 (MQ=255)
                                                           cccGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGaa                                                                                           >  1:592512/1‑150 (MQ=255)
                                                                     ccTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTgg                                                                                 <  2:18972/150‑1 (MQ=255)
                                                                      cTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCt                                                                                   >  1:1054841/1‑147 (MQ=255)
                                                                      cTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCt                                                                                   <  2:1054841/147‑1 (MQ=255)
                                                                      cTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTggg                                                                                >  1:217279/1‑150 (MQ=255)
                                                                      cTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTggg                                                                                >  2:793355/1‑150 (MQ=255)
                                                                      cTGCATCGAGTGGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTggg                                                                                >  1:142918/1‑150 (MQ=255)
                                                                             gAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCa                                                                                                                                                     >  1:1097603/1‑74 (MQ=255)
                                                                             gAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCa                                                                                                                                                     <  2:1097603/74‑1 (MQ=255)
                                                                               gttgttCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGtt                                                                       <  1:793355/150‑1 (MQ=255)
                                                                                  gttCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCa                                                                    <  2:1034499/150‑1 (MQ=255)
                                                                                   ttCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCTGGGCGctgtc                                                                                                                                            >  2:1239576/1‑75 (MQ=255)
                                                                                   ttCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCTGGGCGctgtc                                                                                                                                            <  1:1239576/77‑3 (MQ=255)
                                                                                   ttCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCAt                                                                   >  2:341894/1‑150 (MQ=255)
                                                                                    tCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATa                                                                  <  1:1280478/150‑1 (MQ=255)
                                                                                      aGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGc                                                                                                            >  2:438675/1‑106 (MQ=255)
                                                                                      aGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGc                                                                                                            <  1:438675/106‑1 (MQ=255)
                                                                                      aGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGc                                                                <  2:70242/150‑1 (MQ=255)
                                                                                      aGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGc                                                                >  2:747554/1‑150 (MQ=255)
                                                                                          aagaagTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCCCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGa                                                                                            <  1:288019/118‑1 (MQ=255)
                                                                                          aagaagTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGa                                                                                            >  2:288019/1‑118 (MQ=255)
                                                                                           agaagTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGaa                                                           <  2:404421/150‑1 (MQ=255)
                                                                                            gaagTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGGCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAAc                                                          >  1:1103831/1‑150 (MQ=255)
                                                                                               gtggtgCCGGTGTGCAGCCCGCAGCTGCCGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGc                                                       >  1:482484/1‑150 (MQ=255)
                                                                                                   tgCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATa                                                                  <  2:604386/135‑1 (MQ=255)
                                                                                                   tgCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATa                                                                  >  1:604386/1‑135 (MQ=255)
                                                                                                      cGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACg                                                >  2:1193406/1‑150 (MQ=255)
                                                                                                       ggTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACgg                                               >  2:1203304/1‑150 (MQ=255)
                                                                                                       ggTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACgg                                               <  2:1178574/150‑1 (MQ=255)
                                                                                                       ggTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACgg                                               >  1:673530/1‑150 (MQ=255)
                                                                                                         gtgtGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGAGAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACggg                                              >  2:738496/1‑150 (MQ=255)
                                                                                                               aGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCgcgc                                        <  2:549237/150‑1 (MQ=255)
                                                                                                                 cccGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCCGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCtt                                      >  2:1150769/1‑150 (MQ=255)
                                                                                                                  ccGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCACAACAGCTACCACGGGCCGCGCTTc                                     <  1:747554/150‑1 (MQ=255)
                                                                                                                   cGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCtt                                      <  1:62594/148‑1 (MQ=255)
                                                                                                                   cGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCtt                                      >  2:62594/1‑148 (MQ=255)
                                                                                                                            tgGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCt                           <  2:243550/150‑1 (MQ=255)
                                                                                                                                cagccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGc                                                                <  2:229716/109‑1 (MQ=255)
                                                                                                                                cagccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGc                                                                >  1:229716/1‑109 (MQ=255)
                                                                                                                                cagccaCCGTTGCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGTTTGCATAGCCAGAACAGCTACCACGGGCCGCGCCTCCACACGTTCTACCt                       <  2:142918/150‑1 (MQ=255)
                                                                                                                                  gccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCtgt                     <  1:1003055/150‑1 (MQ=255)
                                                                                                                                  gccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCtgt                     <  2:482484/150‑1 (MQ=255)
                                                                                                                                   ccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTggg                                                                         >  2:114665/1‑97 (MQ=255)
                                                                                                                                   ccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTggg                                                                         <  1:114665/97‑1 (MQ=255)
                                                                                                                                    caCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATGGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTgtgt                   <  1:1201918/150‑1 (MQ=255)
                                                                                                                                       cGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATc                <  1:1193406/150‑1 (MQ=255)
                                                                                                                                         tttCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCcg              >  1:995720/1‑150 (MQ=255)
                                                                                                                                             gACAGCGCCCAGGCGCTGACTGAGCACCGCCTGtt                                                                                                                             <  2:1031161/35‑1 (MQ=255)
                                                                                                                                             gACAGCGCCCAGGCGCTGACTGAGCACCGCCTGtt                                                                                                                             >  1:1031161/1‑35 (MQ=255)
                                                                                                                                              tcAGCGCCTGGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCAGGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAAGAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCCGCGCgg         >  1:683595/2‑150 (MQ=255)
                                                                                                                                              aCAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCtgt                     <  2:57114/138‑1 (MQ=255)
                                                                                                                                              aCAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCtgt                     >  1:57114/1‑138 (MQ=255)
                                                                                                                                                     ccAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCCGCGCGGCCAtctc  <  1:1150769/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CATTCGACCTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCCGCGCGGCCATCTC  >  NC_002947/4642507‑4642804

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: