Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,561,287 G→C 6.2% I102M (ATC→ATG PP_4048 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,561,2870GC 6.2% 361.2 / 12.4 113I102M (ATC→ATGPP_4048hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (56/49);  new base C (4/3);  total (60/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

Only 100 of 120 total aligned reads displayed.

AAACGCAACCACTGTCTGCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGCACA  >  NC_002947/4561140‑4561434
                                                                                                                                                   |                                                                                                                                                   
aaaCGCAACCACTGTCTGCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCgatgat                                                                                                                                                   <  2:1061076/150‑1 (MQ=255)
  aCGCAACCACTGTCTGCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCg                                                                                                                                                 <  1:2208195/150‑1 (MQ=255)
   cGCAACCACTGTCTGCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCgg                                                                                                                                                <  1:1486801/150‑1 (MQ=255)
        ccACTGTCTGCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATCATCGgcgctg                                                                                                                                           <  2:2021482/150‑3 (MQ=255)
         cACTGTCTGCAGGATGACCTGACTTTCACGTATTAGTTTCAAGCACAACCCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGAt                                                                                                                                          >  1:803672/1‑150 (MQ=255)
                  cAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAgg                                                                                                                                      <  2:2326239/145‑1 (MQ=255)
                  cAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAgg                                                                                                                                      >  1:2326239/1‑145 (MQ=255)
                    ggATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGtcatca                                                                                                                               <  2:743712/150‑1 (MQ=255)
                          ccTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATcaggcagg                                                                                                                         >  2:2022937/1‑150 (MQ=255)
                           cTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGtca                                                                                                                                  <  2:579985/140‑1 (MQ=255)
                           cTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGtca                                                                                                                                  >  1:579985/1‑140 (MQ=255)
                           cTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCaggcaggc                                                                                                                        >  1:2213028/1‑150 (MQ=255)
                            tGACTTTCACGCATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGcc                                                                                                                       >  2:1449059/1‑150 (MQ=255)
                                tttCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGccacca                                                                                                                   >  1:68590/1‑150 (MQ=255)
                                tttCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGccacca                                                                                                                   >  2:1924449/1‑150 (MQ=255)
                                tttCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGccacca                                                                                                                   <  1:1460572/150‑1 (MQ=255)
                                tttCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGccacca                                                                                                                   >  2:1170184/1‑150 (MQ=255)
                                tttCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGccacca                                                                                                                   >  1:1487483/1‑150 (MQ=255)
                                      gTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATcaggcag                                                                                                                          >  2:2219887/1‑137 (MQ=255)
                                      gTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATcaggcag                                                                                                                          <  1:2219887/137‑1 (MQ=255)
                                       tATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAacac                                                                                                            <  2:1487483/150‑1 (MQ=255)
                                         tCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATcaggcag                                                                                                                          >  1:1217526/1‑134 (MQ=255)
                                         tCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATcaggcag                                                                                                                          <  2:1217526/134‑1 (MQ=255)
                                         tCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCaggcaggc                                                                                                                        <  1:187377/136‑1 (MQ=255)
                                           aaCAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCAc                                                                                                        >  2:1114093/1‑150 (MQ=255)
                                            aCAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGtcatc                                                                                                                                >  1:109211/1‑125 (MQ=255)
                                            aCAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGtcatc                                                                                                                                <  2:109211/125‑1 (MQ=255)
                                            aCAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATcaggcag                                                                                                                          <  2:1311134/131‑1 (MQ=255)
                                            aCAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATcaggcag                                                                                                                          >  1:1311134/1‑131 (MQ=255)
                                            aCAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCAct                                                                                                       >  1:1770035/1‑150 (MQ=255)
                                                 aagcaCAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGc                                                                                                  <  2:1770035/150‑1 (MQ=255)
                                                 aagcaCAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGGTCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGc                                                                                                  <  1:418751/150‑1 (MQ=255)
                                                     acaAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCAtt                                                                                              >  2:809092/1‑150 (MQ=255)
                                                      caAACCCCGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTc                                                                                             <  1:1114093/150‑1 (MQ=255)
                                                        aaCCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCTCCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCcg                                                                                           >  2:1826993/1‑150 (MQ=255)
                                                         aCCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCaggcaggc                                                                                                                        <  2:1554119/120‑1 (MQ=255)
                                                         aCCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCaggcaggc                                                                                                                        >  1:1554119/1‑120 (MQ=255)
                                                          ccACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCcggc                                                                                         <  2:714206/150‑1 (MQ=255)
                                                          ccACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCcggc                                                                                         <  1:773757/150‑1 (MQ=255)
                                                                   aaCGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTgcggc                                                                                >  2:1655933/1‑150 (MQ=255)
                                                                        tgctgctTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGc                                                                           >  1:1447614/1‑150 (MQ=255)
                                                                            gctTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTgg                                                                       >  2:712274/1‑150 (MQ=255)
                                                                             ctTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCCCAGTTCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGc                                                                      >  1:62664/1‑150 (MQ=255)
                                                                                  tGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGAGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCg                                                                 >  1:1012901/1‑150 (MQ=255)
                                                                                    cccAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATcaggcag                                                                                                                          >  1:1321959/1‑91 (MQ=255)
                                                                                    cccAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATcaggcag                                                                                                                          <  2:1321959/91‑1 (MQ=255)
                                                                                    cccAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCaggcaggc                                                                                                                        <  1:873946/93‑1 (MQ=255)
                                                                                    cccAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCaggcaggc                                                                                                                        >  2:873946/1‑93 (MQ=255)
                                                                                      cAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGcc                                                                  <  2:584851/145‑1 (MQ=255)
                                                                                      cAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGcc                                                                  >  1:584851/1‑145 (MQ=255)
                                                                                      cAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGGAGGTAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTccg                                                             >  1:2222247/1‑150 (MQ=255)
                                                                                       aaCCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAacac                                                                                                            <  1:1091586/102‑1 (MQ=255)
                                                                                            gTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAac                                                                                                              <  2:2103126/95‑1 (MQ=255)
                                                                                            gTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAac                                                                                                              >  1:2103126/1‑95 (MQ=255)
                                                                                            gTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGt                                                       >  1:210644/1‑150 (MQ=255)
                                                                                                gcCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAgcg                                                   >  1:1255702/1‑150 (MQ=255)
                                                                                                   tGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCtt                                                >  2:2195559/1‑150 (MQ=255)
                                                                                                     ccATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATcagg                                                                                                                             >  1:614949/1‑71 (MQ=255)
                                                                                                     ccATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATcagg                                                                                                                             <  2:614949/71‑1 (MQ=255)
                                                                                                      cATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTg                                                                                                   >  2:1859461/1‑96 (MQ=255)
                                                                                                       aTCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGccc                                            >  1:1647458/1‑150 (MQ=255)
                                                                                                        tCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCa                                                                                                         <  1:973759/88‑1 (MQ=255)
                                                                                                           aGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGgc                                        <  1:1900488/150‑1 (MQ=255)
                                                                                                             cccACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGgcgc                                      >  2:1423655/1‑150 (MQ=255)
                                                                                                             cccACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGgcgc                                      >  1:2393134/1‑150 (MQ=255)
                                                                                                               cACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGACCAGGCGCAg                                    <  2:1460502/150‑1 (MQ=255)
                                                                                                                aCGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGccac                                                                                                                     <  1:1063819/68‑1 (MQ=255)
                                                                                                                aCGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGccac                                                                                                                     >  2:1063819/1‑68 (MQ=255)
                                                                                                                aCGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAgg                                   <  2:1373127/150‑1 (MQ=255)
                                                                                                                aCGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAgg                                   >  1:1503775/1‑150 (MQ=255)
                                                                                                                    gAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGcaca                               >  1:2316956/1‑150 (MQ=255)
                                                                                                                        gCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATcag                                                                                                                              >  1:1385307/1‑51 (MQ=255)
                                                                                                                        gCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATcag                                                                                                                              <  2:1385307/51‑1 (MQ=255)
                                                                                                                        gCCCACCGCCCCCACCACAGCGCCGATGATCGGCGAGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGc                           <  2:1012901/150‑1 (MQ=255)
                                                                                                                         cccACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTg                                                                        >  2:2135699/1‑104 (MQ=255)
                                                                                                                         cccACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTg                                                                        <  1:2135699/104‑1 (MQ=255)
                                                                                                                         ccaaCCGCTCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCa                          >  2:836613/4‑150 (MQ=255)
                                                                                                                           cACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGc                        >  1:1205927/1‑150 (MQ=255)
                                                                                                                           cACCGCCCCCACCACAGCGCCGAGTATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGc                        <  2:557951/150‑1 (MQ=255)
                                                                                                                            aCCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGcc                       >  2:1055688/1‑150 (MQ=255)
                                                                                                                              cttccccTCCACAGTTCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAg                     <  2:2316956/147‑1 (MQ=255)
                                                                                                                                   ccctgataGCGCCGATCATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCAc                >  2:710637/8‑150 (MQ=255)
                                                                                                                                   ccaccaCAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCACGCACACCGCAGCCAGGTCAc                >  1:1196904/1‑150 (MQ=255)
                                                                                                                                    caccacAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGc                                                                           <  2:937675/90‑1 (MQ=255)
                                                                                                                                     accacAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACgg              <  1:1423655/150‑1 (MQ=255)
                                                                                                                                          cgcgcCGATCATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTc                                                               <  1:2314385/95‑1 (MQ=255)
                                                                                                                                          cgcgcCGATCATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTc                                                               >  2:2314385/2‑96 (MQ=255)
                                                                                                                                           gcgcCGATCATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGcc                                                                  >  2:557874/1‑92 (MQ=255)
                                                                                                                                           gcgcCGATCATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGcc                                                                  <  1:557874/92‑1 (MQ=255)
                                                                                                                                            cgcCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCCCTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCa                                           <  1:1658494/114‑1 (MQ=255)
                                                                                                                                            cgcCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCa                                           >  2:1658494/1‑114 (MQ=255)
                                                                                                                                             gcCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAg      <  1:1655933/150‑1 (MQ=255)
                                                                                                                                              ccGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGc     <  2:1255702/150‑1 (MQ=255)
                                                                                                                                               cGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGca    <  1:1837283/150‑1 (MQ=255)
                                                                                                                                               cGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTACCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGc     <  1:1259059/149‑1 (MQ=255)
                                                                                                                                                gatgatCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGcac   >  1:2180920/1‑150 (MQ=255)
                                                                                                                                                 atgatCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGcaca  >  2:919670/1‑150 (MQ=255)
                                                                                                                                                 atgatCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGca    <  1:772983/148‑1 (MQ=255)
                                                                                                                                                 atgatCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGca    >  2:772983/1‑148 (MQ=255)
                                                                                                                                                   gatCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCcggc                                                                                         >  2:1262170/1‑61 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAACGCAACCACTGTCTGCAGGATGACCTGACTTTCACGTATCAACAGCAAGCACAAACCACGCAAGAACGCTGCTGCTTCCTGCCCAACCAGTGCGCCTGCCATCGAGCCCACGGGAGTGCCCACCGCCCCCACCACAGCGCCGATGATCGGCGCGATCAGGGTCATCAGGCAGGCCACCAGCAACACCCACTCTTGCCATTCCGGCGGCTGCGGCGTTGCCTGGCTGCCGTCCGCCGCGTAGCGCTTGCCCAGGCGCAGGCACACCGCAGCCAGGTCACGGGCACTGAGCACA  >  NC_002947/4561140‑4561434

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: