Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,230,081 G→C 5.7% P61R (CCG→CGG) ‡ PP_4607 ← transmembrane sensor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,230,0810GC 5.7% 291.7 / 4.8 92P61R (CCG→CGG) ‡PP_4607transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base G (43/43);  new base C (4/2);  total (47/45)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCACCAGCGACGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGCAGGGCAGCCTGGAA  >  NC_002947/5229932‑5230229
                                                                                                                                                     |                                                                                                                                                    
aaCCACCAGCGACGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGccgccg                                                                                                                                                      <  1:655363/150‑1 (MQ=255)
aaCCACCAGCGACGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCACTGCGTCCGCCAGCCGCGccgccg                                                                                                                                                      >  2:1209487/1‑150 (MQ=255)
   caccaGCGACGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCgggg                                                                                                                                                   >  2:451725/1‑150 (MQ=255)
      caGCGACGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGggcgg                                                                                                                                                >  1:1936666/1‑150 (MQ=255)
        gCGACGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGgcggcg                                                                                                                                              >  2:1950737/1‑150 (MQ=255)
        gCGACGCCCACGGCCAGCACCACGCAGGCTGCAAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGgcggcg                                                                                                                                              <  2:1197450/150‑1 (MQ=255)
          gACGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCACCACCCGCCGCCGGGGCGGCGCt                                                                                                                                            >  1:1945024/1‑150 (MQ=255)
           aCGCCCACGGCGAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCCCGGCGgcgcgg                                                                                                                                           <  1:65178/150‑3 (MQ=255)
           aCGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCtt                                                                                                                                           >  2:1857062/1‑150 (MQ=255)
                cACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGGCCGGGCGGCGCTGTTCTg                                                                                                                                      >  2:1130975/1‑150 (MQ=255)
                cACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTg                                                                                                                                      >  2:1055471/1‑150 (MQ=255)
                cACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTg                                                                                                                                      <  2:761025/150‑1 (MQ=255)
                         cACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCAc                                                                                                                             >  2:1963881/1‑150 (MQ=255)
                          aCCAGGCAGGCTGACAGCGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCCAAGCAcc                                                                                                                            >  2:79/1‑150 (MQ=255)
                             aggcaggcTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTgg                                                                                                                         <  2:748549/150‑1 (MQ=255)
                               gcaggcTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGc                                                                                                                       >  1:821374/1‑150 (MQ=255)
                                 aggcTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCt                                                                                                                     >  2:495116/1‑150 (MQ=255)
                                    cTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGt                                                                                                                  >  2:387795/1‑150 (MQ=255)
                                     tGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTc                                                                                                                 <  1:1404381/150‑1 (MQ=255)
                                        cAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTa                                                                                                              >  1:2192616/1‑150 (MQ=255)
                                        cAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTa                                                                                                              <  1:1192603/150‑1 (MQ=255)
                                         aaCGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAg                                                                                                             <  2:1544485/150‑1 (MQ=255)
                                         aaCGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAg                                                                                                             <  2:621418/150‑1 (MQ=255)
                                          aCGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAgg                                                                                                            <  1:1209487/150‑1 (MQ=255)
                                            gCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGTGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCt                                                                                                          >  2:1103034/1‑150 (MQ=255)
                                            gCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCt                                                                                                                                            >  1:1893875/1‑116 (MQ=255)
                                            gCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCt                                                                                                                                            <  2:1893875/116‑1 (MQ=255)
                                            gCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCt                                                                                                          >  2:223385/1‑150 (MQ=255)
                                            gCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCt                                                                                                          >  1:981512/1‑150 (MQ=255)
                                            gCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCt                                                                                                          <  2:821374/150‑1 (MQ=255)
                                                 gCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGcgcg                                                                                                     >  1:586164/1‑150 (MQ=255)
                                                   cTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGgt                                                                                                   <  1:484701/150‑1 (MQ=255)
                                                    tGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGgt                                                                                                   >  1:302327/1‑149 (MQ=255)
                                                    tGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGgt                                                                                                   <  2:302327/149‑1 (MQ=255)
                                                       cACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCCCGGCGgcgc                                                                                                                                             >  1:2358355/1‑104 (MQ=255)
                                                       cACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCCCGGCGgcgc                                                                                                                                             <  2:2358355/104‑1 (MQ=255)
                                                        aCCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCg                                                                                              >  1:2404410/1‑150 (MQ=255)
                                                         cccGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCgg                                                                                             <  2:1441995/150‑1 (MQ=255)
                                                          ccgccgCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCggg                                                                                            <  1:451725/150‑1 (MQ=255)
                                                           cgccgcCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGATGCGAGgtgttcttgc                                                                                           <  2:1913325/150‑7 (MQ=255)
                                                                ccagccaCCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCgggg                                                                                                                                                   <  1:2195061/89‑1 (MQ=255)
                                                                ccagccaCCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCgggg                                                                                                                                                   >  2:2195061/1‑89 (MQ=255)
                                                                 cagccaCCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTc                                                                                                                 >  1:1997152/1‑122 (MQ=255)
                                                                 cagccaCCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTc                                                                                                                 <  2:1997152/122‑1 (MQ=255)
                                                                 cagccaCCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAg                                                                                                             >  2:1313482/1‑126 (MQ=255)
                                                                 cagccaCCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAg                                                                                                             <  1:1313482/126‑1 (MQ=255)
                                                                  agccaCCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCa                                                                                    <  2:114367/150‑1 (MQ=255)
                                                                  agccaCCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCa                                                                                    >  2:2224569/1‑150 (MQ=255)
                                                                      accacGCGCCGGGGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCaa                                                                                >  1:577262/1‑150 (MQ=255)
                                                                       ccacGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAAc                                                                               <  1:796807/150‑1 (MQ=255)
                                                                                 gtggtCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCGGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGa                                                                     >  2:1542684/1‑150 (MQ=255)
                                                                                  tggtCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGaa                                                                    <  1:1950737/150‑1 (MQ=255)
                                                                                   ggtCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAAt                                                                   <  1:495116/150‑1 (MQ=255)
                                                                                   ggtCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTCCGTCCGCCAGCCTCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAAt                                                                   <  1:1985251/150‑1 (MQ=255)
                                                                                   ggtCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCGCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCTGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAAt                                                                   <  2:2335436/150‑1 (MQ=255)
                                                                                      cTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATgcc                                                                <  1:2230201/150‑1 (MQ=255)
                                                                                       tGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGgt                                                                                                   <  2:657631/114‑1 (MQ=255)
                                                                                       tGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGgt                                                                                                   >  1:657631/1‑114 (MQ=255)
                                                                                        ggCCATGGCATCGAGGTAGCGCTGCAGTGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATgccgc                                                              <  1:1103034/150‑1 (MQ=255)
                                                                                           cATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATgccatggc                                                           >  1:422945/1‑145 (MQ=255)
                                                                                                  tCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAg                                                                                   >  2:632002/1‑119 (MQ=255)
                                                                                                  tCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAg                                                                                   <  1:632002/119‑1 (MQ=255)
                                                                                                       gTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGccc                                               >  1:2335849/1‑150 (MQ=255)
                                                                                                       gTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGccc                                               >  1:1271323/1‑150 (MQ=255)
                                                                                                          gcgcTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCggg                                            <  2:1756300/150‑1 (MQ=255)
                                                                                                            gcTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGtt                                          >  1:314994/1‑150 (MQ=255)
                                                                                                             cTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAGCCGGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTg                                         >  1:1706881/1‑150 (MQ=255)
                                                                                                                cAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTc                                      <  2:1945024/150‑1 (MQ=255)
                                                                                                                 aGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGTCCTGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCg                                     <  1:1879137/150‑1 (MQ=255)
                                                                                                                  gggCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCTTAGGCTACGCGGTGTGAGGGGTCGTCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGt                                    >  1:877130/1‑150 (MQ=255)
                                                                                                                  gggCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGt                                    >  1:681379/1‑150 (MQ=255)
                                                                                                                   ggCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGtt                                   <  1:1556122/150‑1 (MQ=255)
                                                                                                                      gTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGAt                                >  2:38034/1‑150 (MQ=255)
                                                                                                                        ccTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTg                                                                                                                    >  2:1874897/1‑64 (MQ=255)
                                                                                                                        ccTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTg                                                                                                                    <  1:1874897/64‑1 (MQ=255)
                                                                                                                         gtGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCg                             >  2:1167269/2‑150 (MQ=255)
                                                                                                                          tGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTTGGCTTGCCCGGGTTGCTCGTTGAtaagg                            >  2:2230112/1‑146 (MQ=255)
                                                                                                                           gCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGt                           <  2:1743119/150‑1 (MQ=255)
                                                                                                                            ctctTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCgg                                                                                             <  1:1480442/83‑1 (MQ=255)
                                                                                                                            ctctTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCgg                                                                                             >  2:1480442/1‑83 (MQ=255)
                                                                                                                            ctctTCGTCCGCCAGCCGCGCCGCAGGGTAATCTCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGtt                          >  1:663500/1‑150 (MQ=255)
                                                                                                                              ctTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCACCCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGTTCCGGTTGa                        <  2:1271323/150‑1 (MQ=255)
                                                                                                                               ttCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGAATAACTTCCTCGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGaa                       <  2:663500/150‑1 (MQ=255)
                                                                                                                                  gTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAAc                                                                               >  2:2067120/1‑91 (MQ=255)
                                                                                                                                  gTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAAc                                                                               <  1:2067120/91‑1 (MQ=255)
                                                                                                                                  gTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCa                    <  1:38034/150‑1 (MQ=255)
                                                                                                                                   tCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCAt                   >  2:558783/1‑150 (MQ=255)
                                                                                                                                      gCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGACAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGc                >  2:1075187/1‑150 (MQ=255)
                                                                                                                                          gccgcgccgcCGGGGCGGCGCTTTTCTGCCAAAGCACCt                                                                                                                           >  2:823444/1‑39 (MQ=255)
                                                                                                                                          gccgcgccgcCGGGGCGGCGCTTTTCTGCCAAAGCACCt                                                                                                                           <  1:823444/39‑1 (MQ=255)
                                                                                                                                          gccgcgccgcCGGGGCGGCGCTTTTCTGCCAAAGCACCTTGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGCAggg            >  2:2125355/1‑150 (MQ=255)
                                                                                                                                          aaagcgccgcCCCGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGACCCATTGCAggg            >  2:56342/4‑150 (MQ=255)
                                                                                                                                            cgcgccgcCGGGGCGGCGCTTTTCTGCCACAGCCCCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGCAGGGCa          <  1:1542684/150‑1 (MQ=255)
                                                                                                                                            cgcgccgcCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGCAGGGCa          <  2:1521255/150‑1 (MQ=255)
                                                                                                                                            cgcgccgcCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGCAGGGCa          <  2:1024690/150‑1 (MQ=255)
                                                                                                                                             gcgccgcCCCGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGACCGGGTTGCTCGTTGATCCGGTTGAACCATTGCAGGGCAg         >  2:227086/1‑150 (MQ=255)
                                                                                                                                              cgccgcCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGCAGGGCAGc        <  1:1055471/150‑1 (MQ=255)
                                                                                                                                                   ccgGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGt                                           <  2:1021737/110‑1 (MQ=255)
                                                                                                                                                   ccgGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGt                                           >  1:1021737/1‑110 (MQ=255)
                                                                                                                                                    cgGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGCAGGGCAGCCTGGaa  >  1:1686982/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AACCACCAGCGACGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGCAGGGCAGCCTGGAA  >  NC_002947/5229932‑5230229

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: