Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,030,800 G→C 8.4% E286D (GAG→GAC) ‡ PP_0887 → sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,030,8000GC 8.4% 211.1 / 7.6 69E286D (GAG→GAC) ‡PP_0887sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (32/31);  new base C (2/4);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TACCCAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCT  >  NC_002947/1030653‑1030940
                                                                                                                                                    |                                                                                                                                            
tACCCAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCTCTGGCGCCCATGAGCGTGCTGATCAACGAAATGCGCCAGGCACACACCGACCCGCTGTTGCAGGAGGACCTGCAACTCCTCCAGGACCAGGTAAACCTGTGCAAagag                                                                                                                                              >  3:27313/1‑148 (MQ=255)
   ccAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACtt                                                                                                                                          >  8:212873/1‑149 (MQ=255)
   ccAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCACTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACtt                                                                                                                                          >  5:231561/1‑149 (MQ=255)
                 ccgcaCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTTGCAAAGTCTCTTTGCACAGCtttac                                                                                                                            <  5:63207/149‑5 (MQ=255)
                  cgccCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTTGCAAAGTCTCTTTGCACAGCtttacc                                                                                                                           <  5:358229/149‑6 (MQ=255)
                   gccCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTgcg                                                                                                                          >  2:95232/1‑149 (MQ=255)
                    cccATGAACTGGGTACGCCACTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTgcgc                                                                                                                         <  6:231561/149‑1 (MQ=255)
                     ccATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGTTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTgcgcg                                                                                                                        >  8:217597/1‑149 (MQ=255)
                      cATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTgcgcg                                                                                                                        >  7:22728/1‑148 (MQ=255)
                       aTGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCa                                                                                                                      <  1:95232/149‑1 (MQ=255)
                          aaCTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGcc                                                                                                                   >  6:231088/1‑149 (MQ=255)
                                gTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGc                                                                                                              >  8:267295/1‑148 (MQ=255)
                                           gCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCt                                                                                                  >  8:163940/1‑149 (MQ=255)
                                           gCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCt                                                                                                  <  5:231088/149‑1 (MQ=255)
                                            cGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTg                                                                                                 <  3:332940/149‑1 (MQ=255)
                                             gACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTgg                                                                                                <  7:212873/149‑1 (MQ=255)
                                             gACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTgg                                                                                                >  7:106160/1‑149 (MQ=255)
                                             gACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTgg                                                                                                >  2:284885/1‑149 (MQ=255)
                                              aCCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGc                                                                                               >  4:300498/1‑149 (MQ=255)
                                                cATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCACCAGCTGGTGCGCGCAGCCGAAGCCCACCGGCGCCTGGc                                                                                               >  8:87182/1‑147 (MQ=255)
                                                  tGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATc                                                                                           >  8:91723/1‑149 (MQ=255)
                                                     gCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTg                                                                                        <  4:56692/149‑1 (MQ=255)
                                                        tGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAg                                                                                     <  7:267295/149‑1 (MQ=255)
                                                              tCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGCAGCAGGAt                                                                               <  8:22728/149‑1 (MQ=255)
                                                                       tGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGcc                                                                      >  6:332296/1‑149 (MQ=255)
                                                                             aGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCctggctg                                                                >  4:340348/1‑149 (MQ=255)
                                                                             aGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCctggctg                                                                <  7:217597/149‑1 (MQ=255)
                                                                                aCCACACCGCCCCGCTGTTGCTGGAGGACGTGCAAATCCTCCAGGACCTGGTATCGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCATCCGCCGCCTTGCGATCGTGGAGGAGGATGGGACCGCCTGGCTGGAc                                                             <  4:27313/149‑1 (MQ=255)
                                                                                 ccACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGAc                                                             >  6:165269/1‑148 (MQ=255)
                                                                                  cacaCCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGAc                                                             >  6:33500/1‑147 (MQ=255)
                                                                                    cacCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAgg                                                         <  7:91723/149‑1 (MQ=255)
                                                                                       cGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGgcgc                                                      <  1:284885/149‑1 (MQ=255)
                                                                                         aCCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGcc                                                                      >  6:280413/1‑131 (MQ=255)
                                                                                         aCCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGcc                                                                      <  5:280413/131‑1 (MQ=255)
                                                                                         aCCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCt                                                     >  5:249194/1‑148 (MQ=255)
                                                                                             gCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAAcc                                                >  4:292192/1‑149 (MQ=255)
                                                                                                 ttGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGcc                                                                                                                   >  1:169998/1‑78 (MQ=255)
                                                                                                 ttGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGcc                                                                                                                   <  2:169998/78‑1 (MQ=255)
                                                                                                 ttGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCt                                             <  4:265719/148‑1 (MQ=255)
                                                                                                   gCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGc                                          >  7:176249/1‑149 (MQ=255)
                                                                                                    cAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCa                                         >  3:120926/1‑149 (MQ=255)
                                                                                                         ggaCCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGa                                    <  5:33500/149‑1 (MQ=255)
                                                                                                         ggaCCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGa                                    >  5:165851/1‑149 (MQ=255)
                                                                                                          gaCCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCctggctg                                                                >  1:321840/1‑120 (MQ=255)
                                                                                                          gaCCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCctggctg                                                                <  2:321840/120‑1 (MQ=255)
                                                                                                            ccTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATgc                                 >  3:192963/1‑149 (MQ=255)
                                                                                                               gCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGcc                              >  7:66404/1‑149 (MQ=255)
                                                                                                                cAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGt                                                                                         >  1:128790/1‑89 (MQ=255)
                                                                                                                cAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGt                                                                                         <  2:128790/89‑1 (MQ=255)
                                                                                                                cAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATgcgc                               <  5:165269/147‑1 (MQ=255)
                                                                                                                  aaTCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCgg                           <  7:190797/149‑1 (MQ=255)
                                                                                                                    tcctccAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGaa                         <  6:165851/149‑1 (MQ=255)
                                                                                                                         cAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCa                     <  6:249194/148‑1 (MQ=255)
                                                                                                                            gACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTa                 <  8:66404/149‑1 (MQ=255)
                                                                                                                            gACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTa                 >  4:129939/1‑149 (MQ=255)
                                                                                                                               cAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCg              <  8:106160/149‑1 (MQ=255)
                                                                                                                                aGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGa                                                           <  4:123713/103‑1 (MQ=255)
                                                                                                                                aGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGa                                                           >  3:123713/1‑103 (MQ=255)
                                                                                                                                aGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCt                                             >  8:64974/1‑117 (MQ=255)
                                                                                                                                aGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCt                                             <  7:64974/117‑1 (MQ=255)
                                                                                                                                   tAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTc          <  7:163940/149‑1 (MQ=255)
                                                                                                                                     caGCTGTTGCAAAGTCTCTTTGCACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGAt                                   <  4:104829/121‑1 (MQ=255)
                                                                                                                                     caGCTGTTGCAAAGTCTCTTTGCACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGAGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGAt                                   >  3:104829/2‑122 (MQ=255)
                                                                                                                                     aaGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCa        <  4:192963/149‑1 (MQ=255)
                                                                                                                                       gctgttgCAAAGTCTCTTTGCACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCt            <  6:167885/138‑1 (MQ=255)
                                                                                                                                       gctgttgCAAAGTCTCTTTGCACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCt            >  5:167885/6‑143 (MQ=255)
                                                                                                                                       gCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAgc      <  3:300498/149‑1 (MQ=255)
                                                                                                                                        cTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAgcg     >  5:39962/1‑149 (MQ=255)
                                                                                                                                            tgCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCt  <  2:164280/149‑1 (MQ=255)
                                                                                                                                            tgCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCt  <  6:39962/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
TACCCAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCT  >  NC_002947/1030653‑1030940

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: