Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,254,998 G→A 5.1% L392L (CTG→TTG)  PP_3728 ← two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,254,9980GA 5.1% 236.0 / 2.3 78L392L (CTG→TTG) PP_3728two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base G (33/40);  new base A (2/2);  total (35/43)
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.77e-01

CGAAGATCTTCGCCTGAACGGCGGGCTCGATACCGATACCATTGTCACGGACCGACAGCCGCAAGCGCTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGTGTCGGCATTCTCCATGGCCATGGGTTCGAGCAACGCCAGCACTTCCTGGCACAACTGAC  >  NC_002947/4254851‑4255141
                                                                                                                                                   |                                                                                                                                               
cGAAGATCTTCCCCTGAACGTCGGGCTCGACACCGATACCATTGTCACGGACCGACAGCCGCAAGCGCTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCcagca                                                                                                                                               <  2:615302/150‑1 (MQ=255)
            ccTGAACGGCGGGCTCGATACCGATACCATTGTCACGGACCGACAGCCGCAAGCGCTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCa                                                                                                                                   >  2:1132409/1‑150 (MQ=255)
             cTGAACGGCGGGCTCGATACCGATACCATTGTCACGGACCGACAGCCGCAAGCGCTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCAc                                                                                                                                  >  1:639839/1‑150 (MQ=255)
              tGAACGGCGGGCTCGATACCGATACCATTGTCACGGACCGACAGCCGCAAGCGCTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCAcc                                                                                                                                 >  2:887111/1‑150 (MQ=255)
                       ggCTCGATACCGATACCATTGTCACGGACCGACAGCCGCAAGCGCTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAgg                                                                                                                        >  1:317968/1‑150 (MQ=255)
                            gataccgataccATTGTCACGGACCGACAGCCGCAAGCGCTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTg                                                                                                                   <  2:639839/150‑1 (MQ=255)
                                 cgataccATTGTCACGGACCGACAGCCGCAAGCGCTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGCACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTg                                                                                                              <  1:667604/150‑1 (MQ=255)
                                   ataccATTGTCACGGACCGACAGCCACAAGCGCTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGAtt                                                                                                                                           <  2:599344/119‑1 (MQ=255)
                                   ataccATTGTCACGGACCGACAGCCACAAGCGCTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGAtt                                                                                                                                           >  1:599344/1‑119 (MQ=255)
                                     accATTGTCACGGACCGACAGCCGCAAGCGCTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGcc                                                                                                          >  2:1184018/1‑150 (MQ=255)
                                         ttGTCACGGACCGACAGCCGCAAGCGCTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTAc                                                                                                      <  2:1015924/150‑1 (MQ=255)
                                              aCGGACCGACAGCCGCAAGCGCTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCGCCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTa                                                                                                 >  2:1093834/1‑150 (MQ=255)
                                                ggACCGACAGCCGCAAGCGCTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTAcg                                                                                               >  2:86119/1‑150 (MQ=255)
                                                  aCCGACAGCCGCAAGCGCTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTAcgcg                                                                                             <  2:126165/150‑1 (MQ=255)
                                                  aCCGACAGCCGCAAGCGCTGCGACCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGAGGGCGTTGGCCATCAGATTGCGCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTAcgcc                                                                                             <  2:317810/150‑2 (MQ=255)
                                                      aaaaccGCAAGAGATGCAATCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCa                                                                                         <  2:1236920/146‑1 (MQ=255)
                                                      aCAGCCGCAAGCGCTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTTGCCAGCAGATTGCTCAGCACCTTTCGCGGGTACTTCTGGGCGCCGTACTGGTAcgccgtcg                                                                                         >  2:589022/1‑145 (MQ=255)
                                                      aCAGCCGCAAGCGCTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTg                                                                                                                   >  2:454193/1‑124 (MQ=255)
                                                      aCAGCCGCAAGCGCTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTg                                                                                                                   <  1:454193/124‑1 (MQ=255)
                                                      aCAGCCGCAAGCGCTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCa                                                                                         <  1:435242/150‑1 (MQ=255)
                                                      aCAGCCGCAAGCGCTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAACAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCa                                                                                         >  2:417788/1‑150 (MQ=255)
                                                           cGCAAGCGCTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGc                                                                                    <  1:589022/150‑1 (MQ=255)
                                                             cAAGCGCTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTc                                                                                  >  2:1203457/1‑150 (MQ=255)
                                                                gcgcTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATc                                                                               <  1:1132409/150‑1 (MQ=255)
                                                                gcgcTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAACAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATc                                                                               <  1:417788/150‑1 (MQ=255)
                                                                         gCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAg                                                                      <  1:1203457/150‑1 (MQ=255)
                                                                           aCTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGtt                                                                    <  1:1184018/150‑1 (MQ=255)
                                                                            cTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTg                                                                   >  2:1190170/1‑150 (MQ=255)
                                                                            cTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTg                                                                   >  1:924540/1‑150 (MQ=255)
                                                                            cTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTg                                                                   <  1:86119/150‑1 (MQ=255)
                                                                             tCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGc                                                                  >  2:294795/1‑150 (MQ=255)
                                                                               gTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAg                                                                <  1:762761/150‑1 (MQ=255)
                                                                                ttATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAgg                                                               <  1:855761/150‑1 (MQ=255)
                                                                                ttATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAgg                                                               <  2:1204115/150‑1 (MQ=255)
                                                                                ttATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTCCGCCGTACTGGTACTCGGCCAGTTGCTCATCGACCTGCAGTTGCAgg                                                               <  1:887111/150‑1 (MQ=255)
                                                                                ttATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGCCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAgg                                                               <  1:764741/150‑1 (MQ=255)
                                                                                      gCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATc                                                                               >  2:793312/1‑128 (MQ=255)
                                                                                      gCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATc                                                                               <  1:793312/128‑1 (MQ=255)
                                                                                      gCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGTGTCg                                                         >  2:1086233/1‑150 (MQ=255)
                                                                                      gCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGTGTCg                                                         >  1:1086300/1‑150 (MQ=255)
                                                                                        aaGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAACAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCt                                                                                   >  2:407682/1‑122 (MQ=255)
                                                                                        aaGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAACAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCt                                                                                   <  1:407682/122‑1 (MQ=255)
                                                                                             ccTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTAGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGGGTCGTCATtct                                                  >  1:661192/1‑150 (MQ=255)
                                                                                              cTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCAc                                                                                                                                  >  1:703146/1‑69 (MQ=255)
                                                                                              cTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCAc                                                                                                                                  <  2:703146/69‑1 (MQ=255)
                                                                                              cTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGTGTCGGCATtctc                                                 >  1:136907/1‑150 (MQ=255)
                                                                                                 tACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCa                                                                 <  2:36349/131‑1 (MQ=255)
                                                                                                 tACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCa                                                                 >  1:36349/1‑131 (MQ=255)
                                                                                                 tACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGTGTCGGCATTCTccat                                              >  1:1234737/1‑150 (MQ=255)
                                                                                                 tACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGTGTCGGCATTCTccat                                              <  1:537988/150‑1 (MQ=255)
                                                                                                 tACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGTGTCGGCATTCTccat                                              <  2:317968/150‑1 (MQ=255)
                                                                                                  aCCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGTGTCGGCATTCTccatg                                             <  1:1086233/150‑1 (MQ=255)
                                                                                                  aCCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGTGTCGGCATCCTccatg                                             <  2:136907/150‑1 (MQ=255)
                                                                                                   ccGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGTGTCGGCATTCTccatgg                                            <  1:1190170/150‑1 (MQ=255)
                                                                                                    cGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTAcgcg                                                                                             >  1:102634/1‑100 (MQ=255)
                                                                                                    cGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTAcgcg                                                                                             <  2:102634/100‑1 (MQ=255)
                                                                                                         gCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCa                                                                 >  1:883223/1‑123 (MQ=255)
                                                                                                         gCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCa                                                                 <  2:883223/123‑1 (MQ=255)
                                                                                                          cTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGgtg                                                            <  1:1194485/127‑1 (MQ=255)
                                                                                                          cTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGgtg                                                            >  2:1194485/1‑127 (MQ=255)
                                                                                                               ccTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGTGTCGGCATTCTCCATGGCCATGGGTTCGa                                <  2:850272/150‑1 (MQ=255)
                                                                                                                      gTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGTGTCGGCATTCTCCATGGCCATGGGTTCGAGCAACGc                         >  1:558327/1‑150 (MQ=255)
                                                                                                                            aGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGTGTCGGCATTCTcc                                                >  1:662801/1‑121 (MQ=255)
                                                                                                                            aGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGTGTCGGCATTCTcc                                                <  2:662801/121‑1 (MQ=255)
                                                                                                                            aGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGTGTCGGCATTCTCCATGGCCATGGGTTCGAGCAACGCCAGCAc                   >  2:248467/1‑150 (MQ=255)
                                                                                                                               gAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGTGTCGGCATTCTCCATGGCCATGGGTTCGAGCAACGCCAGCACTTc                <  1:248467/150‑1 (MQ=255)
                                                                                                                                aaCTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGt                                                                                       <  2:1280660/78‑1 (MQ=255)
                                                                                                                                aaCTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGt                                                                                       >  1:1280660/1‑78 (MQ=255)
                                                                                                                                 aCTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGTGTCGGCATTCTCCATGGCCATGGGTTCGAGCAACGCCAGCACTTCCt              <  1:294795/150‑1 (MQ=255)
                                                                                                                                 aCTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGTGTCGGCATTCTCCATGGCCATGGGTTCGAGCAACGCCAGCACTTCCt              <  2:1086300/150‑1 (MQ=255)
                                                                                                                                     gATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGTGTCGGCATTCTCCATGGCCATGGGTTCGAGCAACGCCAGCACTTCCTGGca          <  1:435796/150‑1 (MQ=255)
                                                                                                                                       tGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAgg                                                               <  2:1250531/95‑1 (MQ=255)
                                                                                                                                       tGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAgg                                                               >  1:1250531/1‑95 (MQ=255)
                                                                                                                                       tGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGTGTCGGCATTCTCCATGGCCATGGGtt                                   <  1:1164131/123‑1 (MQ=255)
                                                                                                                                       tGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGTGTCGGCATTCTCCATGGCCATGGGtt                                   >  2:1164131/1‑123 (MQ=255)
                                                                                                                                       tGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGTGTCGGCATTCTCCATGGCCATGGGTTCGAGCAACGCCAGCACTTCCTGGcaca        >  2:1171583/1‑150 (MQ=255)
                                                                                                                                        ggCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAg                                                                <  1:1134622/93‑1 (MQ=255)
                                                                                                                                        ggCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAg                                                                >  2:1134622/1‑93 (MQ=255)
                                                                                                                                             tGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGTGTCGGCATTCTCCATGGCCATGGGTTCGAGCAACGCCAGCACTTCCTGGCACAACTGAc  <  1:509867/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CGAAGATCTTCGCCTGAACGGCGGGCTCGATACCGATACCATTGTCACGGACCGACAGCCGCAAGCGCTGCGAGCACTCGTTATTGGCAAGCACCTGTACCGTCAGCTGTACCTGGCCGTGCTCAGTGAACTTGATGGCGTTGGCCAGCAGATTGCTCAGCACCTGTCGCAGGTACTGCTCTGCGCCGTACTGGTACGCGGCCAGTTGCTCATCGACCTGCAGTTGCAGGGTGTCGGCATTCTCCATGGCCATGGGTTCGAGCAACGCCAGCACTTCCTGGCACAACTGAC  >  NC_002947/4254851‑4255141

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: