Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,230,082 G→C 5.3% P61A (CCG→GCG) ‡ PP_4607 ← transmembrane sensor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,230,0820GC 5.3% 299.7 / 6.1 95P61A (CCG→GCG) ‡PP_4607transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base G (46/44);  new base C (3/2);  total (49/46)
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.44e-01

CACCAGCGACGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGCAGGGCAGCCTGGAA  >  NC_002947/5229935‑5230229
                                                                                                                                                   |                                                                                                                                                   
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:1192603/150‑1 (MQ=255)
                                     cAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTa                                                                                                              >  1:2192616/1‑150 (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:821374/150‑1 (MQ=255)
                                         gCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCt                                                                                                          >  1:981512/1‑150 (MQ=255)
                                         gCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCt                                                                                                          >  2:223385/1‑150 (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                                                                                   <  1:632002/119‑1 (MQ=255)
                                                                                               tCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAg                                                                                   >  2:632002/1‑119 (MQ=255)
                                                                                                    gTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGccc                                               >  1:1271323/1‑150 (MQ=255)
                                                                                                    gTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGccc                                               >  1:2335849/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)
                                                                                                                                                   |                                                                                                                                                   
CACCAGCGACGCCCACGGCCAGCACCAGGCAGGCTGCCAACGCCAGGCCTGCCACCCGCCGCCAGCCACCACGCGCCGGTGGTCTGGCCATGGCATCGAGGTAGCGCTGCAGGGCGTCCTGCTCTTCGTCCGCCAGCCGCGCCGCCGGGGCGGCGCTTTTCTGCCAAAGCACCTGGGCCTGGTCGTAGGCTGCGCGGTGTGCGGGGTCGGCCAGCAACCAGCGCTTGAATGCCGCCCCTTCGGCTTGCCCGGGTTGCTCGTTGATCCGGTTGAACCATTGCAGGGCAGCCTGGAA  >  NC_002947/5229935‑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: