Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,650,091 T→A 8.4% intergenic (+30/+104) PP_2320 → / ← oprI hypothetical protein/major outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,650,0910TA 8.4% 237.4 / 10.8 83intergenic (+30/+104)PP_2320/oprIhypothetical protein/major outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (38/38);  new base A (5/2);  total (43/40)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

GGACATGGTGCGTGACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCG  >  NC_002947/2649944‑2650233
                                                                                                                                                   |                                                                                                                                              
gcaccgggTGCGTGACCTGGTGGCGCCAGAAGATGGTAGGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCCCACCGAGGGTGACGAGCATTCCGCCCGAACTGCAGAAAGGTAATTTCGGCGTACCCAGGCCTGCCCttc                                                                                                                                               <  4:311087/143‑2 (MQ=255)
  aCATGGTGCGTGACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGt                                                                                                                                             <  7:316174/149‑1 (MQ=255)
   cATGGTGCGTGACGTGGTGGCTGCAGAAGATGGTATGACGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGTCCGCACCGATGGTGACGAGCATTCCTCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTg                                                                                                                                            <  3:20503/149‑1 (MQ=255)
    aTGGTGCGTGACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGACAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGc                                                                                                                                           <  7:189415/149‑1 (MQ=255)
      ggTGCGTGACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAg                                                                                                                                         >  8:202423/1‑149 (MQ=255)
        tGCGTGACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGc                                                                                                                                       <  4:174162/149‑1 (MQ=255)
          cGTGACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGcgt                                                                                                                                     >  8:329020/1‑147 (MQ=255)
           gTGACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGATGGAAAACCAGGGGGCCGCACCGATTGTGACGAGCATTCCGCCAGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCACTAAGGGCAGGCCtg                                                                                                                                    <  8:310109/149‑2 (MQ=255)
             gACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCtttc                                                                                                                                  >  8:100812/1‑148 (MQ=255)
                 tggtggCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTa                                                                                                                              >  3:233908/1‑149 (MQ=255)
                   gtggCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTAcc                                                                                                                            >  2:98360/1‑149 (MQ=255)
                   gtggCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTAcc                                                                                                                            >  6:175185/1‑149 (MQ=255)
                   gtggCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTAc                                                                                                                             <  7:202423/148‑1 (MQ=255)
                        tGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACctgct                                                                                                                        >  2:139743/1‑148 (MQ=255)
                        tGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACctgct                                                                                                                        >  8:14267/1‑148 (MQ=255)
                          cAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGa                                                                                                                      <  7:143255/148‑1 (MQ=255)
                           agaagaTGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGAcc                                                                                                                    >  5:116255/1‑149 (MQ=255)
                           agaagaTGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGAcc                                                                                                                    <  3:280087/149‑1 (MQ=255)
                                aTGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGtt                                                                                                               <  6:9595/149‑1 (MQ=255)
                                        ccGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGtt                                                                                                       >  3:219296/1‑149 (MQ=255)
                                        ccGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGt                                                                                                        <  2:308131/148‑1 (MQ=255)
                                         cGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGttt                                                                                                      >  3:291004/1‑149 (MQ=255)
                                          ggTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACctg                                                                                                                          >  5:260583/1‑128 (MQ=255)
                                          ggTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACctg                                                                                                                          <  6:260583/128‑1 (MQ=255)
                                           gTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAg                                                                                                                                         >  5:320948/1‑112 (MQ=255)
                                           gTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAg                                                                                                                                         <  6:320948/112‑1 (MQ=255)
                                           gTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACcgc                                                                                                                  <  7:332438/135‑1 (MQ=255)
                                           gTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACcgc                                                                                                                  >  8:332438/1‑135 (MQ=255)
                                           gTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTc                                                                                                     >  3:93724/1‑148 (MQ=255)
                                           gTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTc                                                                                                     >  5:193457/1‑148 (MQ=255)
                                            tAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGc                                                                                                                                           >  2:31702/1‑109 (MQ=255)
                                                aaTTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAgg                                                                                                >  7:6272/1‑148 (MQ=255)
                                                aaTTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGc                                                                                               <  7:14267/149‑1 (MQ=255)
                                                  ttGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCaa                                                                                             >  6:87593/1‑149 (MQ=255)
                                                    gccgtgccgGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGGCGGTTTCAAAGGCAAt                                                                                            >  8:57972/1‑148 (MQ=255)
                                                          ccgGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATaaaaaa                                                                                      <  4:227687/148‑1 (MQ=255)
                                                          ccgGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATaaaaaa                                                                                      <  4:291004/148‑1 (MQ=255)
                                                          ccgGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATaaaaaa                                                                                      <  5:175185/148‑1 (MQ=255)
                                                           cgGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGc                                                                                    <  8:6272/149‑1 (MQ=255)
                                                                gAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGAcc                                                                               <  3:42039/149‑1 (MQ=255)
                                                                  aaaCCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGACGACCca                                                                             >  6:347794/1‑149 (MQ=255)
                                                                   aaCCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCca                                                                             >  1:2379/1‑148 (MQ=255)
                                                                          gggCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGAcc                                                                               >  5:194364/1‑139 (MQ=255)
                                                                          gggCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGAcc                                                                               <  6:194364/139‑1 (MQ=255)
                                                                            gCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTc                                                                   <  4:233908/149‑1 (MQ=255)
                                                                              cGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCATTAAAAAAGCCGACCCACTTGAGGGTCgg                                                                 <  1:200237/149‑1 (MQ=255)
                                                                              cGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCgg                                                                 >  5:1344/1‑149 (MQ=255)
                                                                               gCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGc                                                                <  2:2379/149‑1 (MQ=255)
                                                                                cACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTCCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACTTGTGGGTCGGct                                                               <  7:100812/149‑1 (MQ=255)
                                                                                 aCCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGctt                                                              <  6:193457/149‑1 (MQ=255)
                                                                                   cGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGcttct                                                            <  7:212235/149‑1 (MQ=255)
                                                                                   cGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGcttct                                                            <  4:204761/149‑1 (MQ=255)
                                                                                     ctGGTCACGAGCATTCCGCCCGACCTGCACTAAGGTAATTTCGCCGTTCCCAGGCCGGCCCTTCGGGCAGGCCTTTTCCTACCTGCCGACCGCGTTTGCCGGTTTCAAAGGCTATAAAACACCCGACCCACAAGTGGGGCGGCTTCTTa                                                          <  2:286039/148‑1 (MQ=255)
                                                                                     aTGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGcttctt                                                           >  2:50942/1‑148 (MQ=255)
                                                                                         tGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTTACa                                                       >  8:120665/1‑148 (MQ=255)
                                                                                              aGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGt                                                 <  4:219296/149‑1 (MQ=255)
                                                                                                cATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGt                                                 >  4:46764/1‑147 (MQ=255)
                                                                                                cATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGa                                               <  1:50942/149‑1 (MQ=255)
                                                                                                cATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGa                                               <  4:93724/149‑1 (MQ=255)
                                                                                                      gCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGatta                                         >  3:343690/1‑149 (MQ=255)
                                                                                                      gCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGatta                                         >  1:90653/1‑149 (MQ=255)
                                                                                                                 aGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGc                                 <  6:97381/146‑1 (MQ=255)
                                                                                                                 aGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGc                                 >  5:97381/1‑146 (MQ=255)
                                                                                                                   tAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCggctg                            <  1:139743/149‑1 (MQ=255)
                                                                                                                      ggTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCAGAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCggctggct                         <  7:57972/149‑1 (MQ=255)
                                                                                                                       gTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGTTGGCtt                        >  6:162743/1‑149 (MQ=255)
                                                                                                                        tAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCttt                       <  2:90653/149‑1 (MQ=255)
                                                                                                                        tAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCttt                       <  4:343690/149‑1 (MQ=255)
                                                                                                                            ttCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTc                    >  3:196282/1‑148 (MQ=255)
                                                                                                                              cGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCa                  <  5:347794/148‑1 (MQ=255)
                                                                                                                                gCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTc                                                                   <  2:344175/97‑1 (MQ=255)
                                                                                                                                gCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTc                                                                   >  1:344175/1‑97 (MQ=255)
                                                                                                                                gCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTg                                  >  6:365879/1‑130 (MQ=255)
                                                                                                                                gCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCa               >  5:183863/1‑149 (MQ=255)
                                                                                                                                gCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTCCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTg                                  <  5:365879/130‑1 (MQ=255)
                                                                                                                                    gaaaaGGCCTGCACTAAGGGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACg           >  7:212099/5‑149 (MQ=255)
                                                                                                                                      ccAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACg           >  3:160671/1‑147 (MQ=255)
                                                                                                                                       aaGGCCTGCACTAAGGGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAg        >  5:319653/2‑149 (MQ=255)
                                                                                                                                        aGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAgg       >  4:273425/1‑149 (MQ=255)
                                                                                                                                         ggCCTGCACTAAGGGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAggg      >  7:197550/1‑149 (MQ=255)
                                                                                                                                          gCCTGCACTAAGGGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGgc     >  2:292195/1‑149 (MQ=255)
                                                                                                                                             tGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGgcgcg  >  3:145756/1‑149 (MQ=255)
                                                                                                                                             tGCACTAAGGGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGc                >  7:36653/1‑135 (MQ=255)
                                                                                                                                             tGCACTAAGGGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGc                <  8:36653/135‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
GGACATGGTGCGTGACGTGGTGGCTGCAGAAGATGGTATGCCGGTAGAAATTGCCGTGCCGGTGGAAAACCAGGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCG  >  NC_002947/2649944‑2650233

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: