Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,257,489 T→C 7.9% intergenic (+27/+216) ureG → / ← PP_2851 urease accessory protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,257,4890TC 7.9% 966.5 / 46.1 316intergenic (+27/+216)ureG/PP_2851urease accessory protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (148/142);  new base C (17/8);  total (165/151)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

Only 100 of 319 total aligned reads displayed.

CGACACCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGG  >  NC_002947/3257342‑3257638
                                                                                                                                                   |                                                                                                                                                     
cgACACCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTcc                                                                                                                                                     <  1:6765383/150‑1 (MQ=255)
 gACACCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCg                                                                                                                                                    >  1:5373047/1‑150 (MQ=255)
 gACACCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCg                                                                                                                                                    >  1:5787543/1‑150 (MQ=255)
   cacCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGca                                                                                                                                                  <  1:1587409/150‑1 (MQ=255)
   cacCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGca                                                                                                                                                  <  2:4506587/150‑1 (MQ=255)
    acCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGcac                                                                                                                                                 <  2:3198114/150‑1 (MQ=255)
      ccAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGTCTCTGTGCGGGCATGTCCGCACag                                                                                                                                               >  2:6997301/1‑150 (MQ=255)
      ccAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACag                                                                                                                                               >  2:2695672/1‑150 (MQ=255)
      ccAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACag                                                                                                                                               >  1:4076638/1‑150 (MQ=255)
      ccAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACag                                                                                                                                               >  2:5898756/1‑150 (MQ=255)
       cAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACaga                                                                                                                                              <  1:2617120/150‑1 (MQ=255)
       cAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACaga                                                                                                                                              <  1:732858/150‑1 (MQ=255)
       cAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACaga                                                                                                                                              <  2:3930655/150‑1 (MQ=255)
         gcgcATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAgaga                                                                                                                                            >  2:1101603/1‑150 (MQ=255)
                cGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGTCTCTGTGCGGACATGCCCGCACAGAGTcctcgat                                                                                                                                     <  1:1504701/150‑6 (MQ=39)
                cGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGAcctcgat                                                                                                                                     <  2:5325487/150‑6 (MQ=255)
                           gCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCa                                                                                                                          <  1:6077367/150‑1 (MQ=255)
                             ccTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCAcc                                                                                                                        >  2:1890684/1‑150 (MQ=255)
                               tGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGTGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGGGGCTGGCACCgg                                                                                                                      >  1:2076482/1‑150 (MQ=255)
                               tGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTGGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCgg                                                                                                                      <  1:5147147/150‑1 (MQ=255)
                                 gACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGct                                                                                                                    <  2:4076638/150‑1 (MQ=255)
                                 gACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGct                                                                                                                    >  2:4678852/1‑150 (MQ=255)
                                  aCCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGctc                                                                                                                   >  1:4919675/1‑150 (MQ=255)
                                    cTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTg                                                                                                                 >  1:87062/1‑150 (MQ=255)
                                    cTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTg                                                                                                                 >  1:7013091/1‑150 (MQ=255)
                                        aGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCgg                                                                                                             >  1:1827542/1‑150 (MQ=255)
                                            aCCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTg                                                                                                                 <  1:3879431/142‑1 (MQ=255)
                                             ccTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTc                                                                                                        >  2:7235577/1‑150 (MQ=255)
                                             ccTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGACACCGGCAGAGCCGgtgcca                                                                                                        <  1:5321878/150‑4 (MQ=255)
                                              cTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTcg                                                                                                       >  2:5151275/1‑150 (MQ=255)
                                               tGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGtg                                                                                                                                 <  2:4233116/123‑1 (MQ=255)
                                               tGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTcgc                                                                                                      >  2:3624647/1‑150 (MQ=255)
                                               tGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTcgc                                                                                                      >  2:2261757/1‑150 (MQ=255)
                                                gaagaaGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTcgcg                                                                                                     >  2:4871602/1‑150 (MQ=255)
                                                gaagaaGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACAAACGCGGtggct                                                                                                                              >  1:196145/1‑125 (MQ=255)
                                                 aagaagGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCgg                                                                                                    <  1:4224160/150‑1 (MQ=255)
                                                 aagaagGGCGAGGGGTTGCAGGCCGTGATCGATGTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCgg                                                                                                    <  2:4359240/150‑1 (MQ=255)
                                                   gaagGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGt                                                                                                  >  2:3568084/1‑150 (MQ=255)
                                                      gggCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTa                                                                                               >  1:6307077/1‑150 (MQ=255)
                                                           aGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCcgc                                                                                          >  1:1069600/1‑150 (MQ=255)
                                                           aGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCcgc                                                                                          <  2:3388488/150‑1 (MQ=255)
                                                            ggggTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCcgcg                                                                                         >  1:60074/1‑150 (MQ=255)
                                                             gggTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGa                                                                                        >  2:4298362/1‑150 (MQ=255)
                                                               gTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAAt                                                                                      >  1:6653879/1‑150 (MQ=255)
                                                               gTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAAt                                                                                      >  1:3757636/1‑150 (MQ=255)
                                                                ttGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATa                                                                                     >  2:3801968/1‑150 (MQ=255)
                                                                 tGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTATTCGCGGGTAAACCCGCGAacac                                                                                    <  1:4762895/150‑1 (MQ=255)
                                                                     ggCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTa                                                                                <  1:6608535/150‑1 (MQ=255)
                                                                     ggCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTa                                                                                >  2:1392541/1‑150 (MQ=255)
                                                                     ggCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTa                                                                                <  2:2031043/150‑1 (MQ=255)
                                                                       ccGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGTGGCTGACATATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGc                                                                              <  1:6237769/150‑1 (MQ=255)
                                                                           gATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGcc                                                                          >  2:251561/1‑150 (MQ=255)
                                                                            aTCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCg                                                                         >  2:2819667/1‑150 (MQ=255)
                                                                               gATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGGACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGcaa                                                                      >  2:1959321/1‑150 (MQ=255)
                                                                               gATTGTCTCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGcaa                                                                      <  2:6686218/150‑1 (MQ=255)
                                                                               caTTTTATCGACGAGCGCGGGATGCTGGGGGTGCGAGGCTGACACATCGAGGACTCTGTGCGGGCATGTCCCCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGcaa                                                                      <  1:4627825/149‑1 (MQ=255)
                                                                                  tttATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGcaacaa                                                                   <  2:310723/150‑1 (MQ=255)
                                                                                   ttATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGTTCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAAt                                                                  <  1:3514133/150‑1 (MQ=255)
                                                                                      tcgtcgAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGaa                                                               >  2:5718129/1‑150 (MQ=255)
                                                                                         tcgAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAAcgg                                                            >  1:3937862/1‑150 (MQ=255)
                                                                                         tcgAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAAcgg                                                            >  1:3337260/1‑150 (MQ=255)
                                                                                         tcgAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAAcgg                                                            >  2:223879/1‑150 (MQ=255)
                                                                                          cgAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAAcggc                                                           >  2:2983295/1‑150 (MQ=255)
                                                                                           gAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTCCCCGCGAATACGGTAGCGGCCGCAACAATGAAcggcg                                                          <  2:6394981/150‑1 (MQ=255)
                                                                                           gAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTCGCGGCCGCAACAATGAAcggcg                                                          >  2:3268793/1‑150 (MQ=255)
                                                                                              cgcgGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGt                                                       <  1:2437144/150‑1 (MQ=255)
                                                                                                   gATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGaaa                                                  <  1:3577044/150‑1 (MQ=255)
                                                                                                   gATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGACATGCCCGcaca                                                                                                                                                >  1:3047835/1‑56 (MQ=37)
                                                                                                      ggtGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCTTGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTg                                               <  1:5902933/148‑1 (MQ=255)
                                                                                                      gCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTg                                               >  2:6065862/1‑150 (MQ=255)
                                                                                                      gCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTg                                               <  1:4298362/150‑1 (MQ=255)
                                                                                                      gCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTg                                               >  1:2792332/1‑150 (MQ=255)
                                                                                                        tGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGc                                             <  2:1735077/150‑1 (MQ=255)
                                                                                                        tGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGc                                             >  1:5198465/1‑150 (MQ=255)
                                                                                                           gggTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCcag                                          <  2:6801915/150‑1 (MQ=255)
                                                                                                                  aGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGggcc                                   >  1:2995278/1‑150 (MQ=255)
                                                                                                                            aTCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTcg                         >  1:3786914/1‑150 (MQ=255)
                                                                                                                             tcgtgtaCTCGGGGCGGGCTTGTCCGCACAGAGACAAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTcgc                        <  1:5023849/144‑1 (MQ=37)
                                                                                                                             tCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAAAGAAACGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTcgc                        >  1:1444649/1‑150 (MQ=255)
                                                                                                                              gttggtctctGTGCGGACATGCCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTcgcg                       >  2:1167772/7‑150 (MQ=255)
                                                                                                                              cGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTcgcg                       <  1:7120511/150‑1 (MQ=255)
                                                                                                                               gAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCgg                      <  1:4438023/150‑1 (MQ=255)
                                                                                                                               gAGGACTCTGTGCGAGCATGTCCGCACAGAGACCACCGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCgg                      <  1:1256901/150‑1 (MQ=255)
                                                                                                                               gAGGACTCGGTGCGGGCCGGTCCGCACAGAGACCACCGCGGTGGCTGGCACCGGCTCGGCCGGTGTTCGCGGGTTGACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCgg                      <  2:3716392/150‑1 (MQ=37)
                                                                                                                                    ctctGTGCGGACATGCCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTaaa                 >  1:738743/1‑150 (MQ=255)
                                                                                                                                     tctGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAc                >  2:3698409/1‑150 (MQ=255)
                                                                                                                                     tctGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAc                >  2:4453323/1‑150 (MQ=255)
                                                                                                                                     tctGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTGTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAc                <  2:1274396/150‑1 (MQ=255)
                                                                                                                                     tctGTGCGGACATGCCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAc                >  1:4505704/1‑150 (MQ=255)
                                                                                                                                     tctGTGCGGACATGCCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAc                >  2:689064/1‑150 (MQ=255)
                                                                                                                                       tgtgCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAGTACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAccc              >  2:1599760/1‑150 (MQ=255)
                                                                                                                                        gtgCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTCAACCCg             >  2:5616363/1‑150 (MQ=255)
                                                                                                                                         tgCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGc            >  1:2189968/1‑150 (MQ=255)
                                                                                                                                         tgCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGc            >  2:4471296/1‑150 (MQ=255)
                                                                                                                                          gCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGCACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCt           <  2:4751721/150‑1 (MQ=255)
                                                                                                                                          gCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCt           <  1:4848238/150‑1 (MQ=255)
                                                                                                                                            gggCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTcc         >  1:6643172/1‑150 (MQ=255)
                                                                                                                                             ggCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCt        >  1:5445721/1‑150 (MQ=255)
                                                                                                                                                 tGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACgg    <  1:6486135/150‑1 (MQ=255)
                                                                                                                                                   tCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACgggg  <  1:5100558/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGACACCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGG  >  NC_002947/3257342‑3257638

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