Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,823,216 A→G 5.2% intergenic (+123/+14) PP_4242 → / ← pvdL conserved protein of unknown function with Rhs repeat‑associated core domain/non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,823,2160AG 5.2% 939.1 / 23.2 288intergenic (+123/+14)PP_4242/pvdLconserved protein of unknown function with Rhs repeat‑associated core domain/non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (153/120);  new base G (5/10);  total (158/130)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 290 total aligned reads displayed.

CTAAGGGAACCGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCG  >  NC_002947/4823067‑4823363
                                                                                                                                                     |                                                                                                                                                   
  aaGGGAACCGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAgc                                                                                                                                                   <  2:2082002/150‑1 (MQ=255)
  aaGGGAACCGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAgc                                                                                                                                                   >  1:2391072/1‑150 (MQ=255)
   aGGGAACCGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAgcg                                                                                                                                                  >  2:5379719/1‑150 (MQ=255)
    gggAACCGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAgcgc                                                                                                                                                 <  1:4066099/150‑1 (MQ=255)
    gggAACCGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAgcgc                                                                                                                                                 >  1:3279310/1‑150 (MQ=255)
     ggAACCGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGcc                                                                                                                                                <  2:2821569/150‑1 (MQ=255)
       aaCCGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATAAAAGGGGCGCCAAGCGcccc                                                                                                                                              >  1:3931943/1‑150 (MQ=255)
        aCCGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCTTGGCGCCCCt                                                                                                                                             <  2:3760815/150‑1 (MQ=255)
        aCCGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCt                                                                                                                                             >  2:857405/1‑150 (MQ=255)
        aCCGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCt                                                                                                                                             <  2:6022552/150‑1 (MQ=255)
           gTTTGGCTATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCACTGAAAGGGGCGCTTGGCGCCCCtttc                                                                                                                                          <  2:5983628/150‑4 (MQ=255)
                gCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATca                                                                                                                                     <  2:1723107/150‑1 (MQ=255)
                 cGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATcac                                                                                                                                    >  2:4613314/1‑150 (MQ=255)
                       ttaGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGcccc                                                                                                                              >  1:3636700/1‑150 (MQ=255)
                         aGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTc                                                                                                                            <  2:1654138/150‑1 (MQ=255)
                            ccTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAg                                                                                                                         >  1:4554580/1‑150 (MQ=255)
                            ccTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAg                                                                                                                         <  1:6586871/150‑1 (MQ=255)
                             cTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGt                                                                                                                        >  1:730966/1‑150 (MQ=255)
                              tAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGtt                                                                                                                       <  2:6212449/150‑1 (MQ=255)
                              tAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGtt                                                                                                                       >  2:5304870/1‑150 (MQ=255)
                              tAGCTTAGAGGGCCTCATCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCAGCAAATCCGGATGAACGACTCAGCCGCAGTCAGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGtt                                                                                                                       >  2:6807246/1‑150 (MQ=255)
                                     gAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAt                                                                                                                <  2:566998/150‑1 (MQ=255)
                                       gggCCTCTTCACTAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCa                                                                                                              <  1:4781765/150‑1 (MQ=255)
                                       gggCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCa                                                                                                              <  1:5527317/150‑1 (MQ=255)
                                          ccTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGt                                                                                                           >  2:3206685/1‑150 (MQ=255)
                                          ccTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGt                                                                                                           <  1:3838669/150‑1 (MQ=255)
                                          ccTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGt                                                                                                           >  1:2141232/1‑150 (MQ=255)
                                            tctTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCa                                                                                                         >  2:2196708/1‑150 (MQ=255)
                                            tctTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCa                                                                                                         <  2:5597629/150‑1 (MQ=255)
                                            tctTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCa                                                                                                         >  2:5927847/1‑150 (MQ=255)
                                            tctTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCa                                                                                                         >  1:2853504/1‑150 (MQ=255)
                                              ttCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACt                                                                                                       <  2:4888987/150‑1 (MQ=255)
                                              ttCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACt                                                                                                       <  2:1296833/150‑1 (MQ=255)
                                               tCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTg                                                                                                      >  2:5438745/1‑150 (MQ=255)
                                                 aCGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGcc                                                                                                                                                >  1:4183887/1‑106 (MQ=255)
                                                 aCGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGcc                                                                                                                                                <  2:4183887/106‑1 (MQ=255)
                                                  cGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGc                                                                                                   <  1:140989/150‑1 (MQ=255)
                                                   gAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGcc                                                                                                  >  2:950951/1‑150 (MQ=255)
                                                      aaTAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGt                                                                                               <  2:5331076/150‑1 (MQ=255)
                                                         aGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGa                                                                                            >  1:4752200/1‑150 (MQ=255)
                                                          gCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGAc                                                                                           <  2:6485341/150‑1 (MQ=255)
                                                           caaGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGAcc                                                                                          >  1:6306143/3‑150 (MQ=255)
                                                            cAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCt                                                                                         >  2:1075106/1‑150 (MQ=255)
                                                             aGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCtt                                                                                        >  1:1671209/1‑150 (MQ=255)
                                                                aaCGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGtc                                                                                     >  1:5154878/1‑150 (MQ=255)
                                                                aaCGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGtc                                                                                     <  1:2883485/150‑1 (MQ=255)
                                                                   gAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGtcatc                                                                                  <  2:4792259/150‑1 (MQ=255)
                                                                   gAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGtcatc                                                                                  >  1:5406214/1‑150 (MQ=255)
                                                                    aGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGtcatca                                                                                 <  1:4028316/150‑1 (MQ=255)
                                                                    aGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGtcatca                                                                                 >  1:2769589/1‑150 (MQ=255)
                                                                        aCTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCa                                                                             <  2:4047582/150‑1 (MQ=255)
                                                                          tGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAgc                                                                           <  1:6272450/150‑1 (MQ=255)
                                                                          tGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAgc                                                                           <  1:1493080/150‑1 (MQ=255)
                                                                          tGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCAGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAgc                                                                           >  1:6118420/1‑150 (MQ=255)
                                                                          tGCCTGATTGATCTGCAAATCCGGACGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAgc                                                                           <  1:876591/150‑1 (MQ=255)
                                                                               gattgatCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCTTGGCGCCCCtttca                                                                                                                                         <  2:5239615/83‑5 (MQ=38)
                                                                               gattgatCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCTTGGCGCCCCtttca                                                                                                                                         >  1:5239615/1‑79 (MQ=38)
                                                                               gattgatCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACt                                                                      <  2:6339454/150‑1 (MQ=255)
                                                                                attgatCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTg                                                                     >  2:6781192/1‑150 (MQ=255)
                                                                                    atCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCt                                                                 <  1:1234521/150‑1 (MQ=255)
                                                                                    atCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCt                                                                 <  1:2079576/150‑1 (MQ=255)
                                                                                     tCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTg                                                                >  2:745338/1‑150 (MQ=255)
                                                                                     tCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACACGACCAATGAAAGGGGCGCCAAGCGCCCCAGCCCATCACAGCCCCTCCAGTTCCGCCATCAGCTCACTGAGCCGGACGACCTAGTCATCATCCAGCGCACAGCCCTg                                                                >  2:5794760/1‑150 (MQ=255)
                                                                                       tGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCa                                                              >  2:6599913/1‑150 (MQ=255)
                                                                                         cAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCGTCATCCAGCGCACTGCCCTGCAGc                                                            >  1:2306243/1‑150 (MQ=255)
                                                                                         cAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGc                                                            <  2:4374139/150‑1 (MQ=255)
                                                                                           aaTCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGccc                                                          <  1:6398406/150‑1 (MQ=255)
                                                                                           aaTCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGccc                                                          <  2:1874574/150‑1 (MQ=255)
                                                                                            aTCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGcccc                                                         <  2:4865822/150‑1 (MQ=255)
                                                                                             tCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCt                                                        >  2:3614845/1‑150 (MQ=255)
                                                                                             tCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCt                                                        >  1:3456235/1‑150 (MQ=255)
                                                                                                ggATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTa                                                     >  1:4386950/1‑150 (MQ=255)
                                                                                                ggATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTa                                                     >  1:6830940/1‑150 (MQ=255)
                                                                                                 gATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTAc                                                    <  1:2963213/150‑1 (MQ=255)
                                                                                                   tggaCGACTCAGCCGCAGTGTGCGCAGCCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGGTCCGGCATCAGGTCACTGCGCCGGTCGACGTGGTCATCATCCTGCGCTCTGCCCTGCAGCCCCTGTGCAt                                                  >  2:2479233/4‑150 (MQ=255)
                                                                                                     aaCGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATAc                                                >  2:5469208/1‑150 (MQ=255)
                                                                                                      aCGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATAAGGTCACTGAGCCGGTCGACCTTGTCATCATACAGCGCACTGCCCTGCAGCCCCTGTACATACt                                               >  1:4947770/1‑150 (MQ=255)
                                                                                                       cGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTc                                              <  1:4116969/150‑1 (MQ=255)
                                                                                                         aCTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCgg                                            <  1:5656482/150‑1 (MQ=255)
                                                                                                           tCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGcc                                          >  1:1110605/1‑150 (MQ=255)
                                                                                                           tCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGcc                                          >  1:53430/1‑150 (MQ=255)
                                                                                                              gCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAAc                                       <  1:4293142/150‑1 (MQ=255)
                                                                                                                 gCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACAcc                                    <  1:6667105/150‑1 (MQ=255)
                                                                                                                      gtgCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCACCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCAc                               <  2:117445/150‑1 (MQ=255)
                                                                                                                         cgcAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGt                            >  1:2580123/1‑150 (MQ=255)
                                                                                                                          gcAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGtg                           >  2:4276070/1‑150 (MQ=255)
                                                                                                                          gcAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGtg                           >  1:2655153/1‑150 (MQ=255)
                                                                                                                              cccGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGtgctgc                       >  2:2470262/1‑150 (MQ=255)
                                                                                                                                 gACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACt                    >  1:3917090/1‑150 (MQ=255)
                                                                                                                                 gACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACt                    >  2:1929086/1‑150 (MQ=255)
                                                                                                                                    cAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGa                 >  2:3439282/1‑150 (MQ=255)
                                                                                                                                    cAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCAGCCGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACAGACTCGGCCAACGCCTCGACCGTGCTGGGCTCGa                 >  2:2291409/1‑150 (MQ=255)
                                                                                                                                       tGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACa              <  2:3279310/150‑1 (MQ=255)
                                                                                                                                         gaaGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATc            >  1:1591929/2‑150 (MQ=255)
                                                                                                                                          caGGGGCGCTTGGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCg           >  2:4469747/2‑150 (MQ=255)
                                                                                                                                          caGGGGCGCTTGGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCg           >  1:5917460/2‑150 (MQ=255)
                                                                                                                                          aaGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACATTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCg           >  1:2492820/1‑150 (MQ=255)
                                                                                                                                             gggCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGccc        <  2:5856882/150‑1 (MQ=255)
                                                                                                                                                cgcCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCa     <  2:2391072/150‑1 (MQ=255)
                                                                                                                                                   cAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCATCCACCGTGCTGCACTCGAACATCGCCCGAAGCg  <  2:2753972/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CTAAGGGAACCGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCG  >  NC_002947/4823067‑4823363

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