Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,650,161 A→T 5.7% intergenic (+100/+34) 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,1610AT 5.7% 208.7 / 2.5 70intergenic (+100/+34)PP_2320/oprIhypothetical protein/major outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (36/29);  new base T (2/2);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAGCCTTACGGTAGGCT  >  NC_002947/2650016‑2650306
                                                                                                                                                 |                                                                                                                                                 
gggggCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTg                                                                                                                                                >  4:86240/1‑149 (MQ=255)
    gCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTc                                                                                                                                            >  8:275955/1‑149 (MQ=255)
     ccGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGt                                                                                                                                                 <  5:174693/143‑1 (MQ=255)
     ccGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGt                                                                                                                                                 >  6:174693/1‑143 (MQ=255)
      cGCACCGATGGTGACTAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGGTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACTTGTGGGTCgg                                                                                                                                          <  6:272451/149‑1 (MQ=255)
      cGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACTTGTGGGTCgg                                                                                                                                          <  5:113704/149‑1 (MQ=255)
         aCCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCTGACCCACAAGTGGGTCGGct                                                                                                                                        >  2:238875/1‑148 (MQ=255)
            gATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGcttctt                                                                                                                                    <  8:15743/149‑1 (MQ=255)
             aTGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGcttc                                                                                                                                      <  5:136659/146‑1 (MQ=255)
             aTGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGcttc                                                                                                                                      >  6:136659/1‑146 (MQ=255)
              tGGTGACGAGCATTCTTCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTaa                                                                                                                                  <  1:275352/149‑1 (MQ=255)
                 tGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACaa                                                                                                                               >  7:128782/1‑149 (MQ=255)
                    cGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATcc                                                                                                                            >  8:240086/1‑149 (MQ=255)
                     gAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGcttttttat                                                                                                                                 <  3:263805/143‑4 (MQ=255)
                     gAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGcttttttat                                                                                                                                 >  4:263805/1‑140 (MQ=255)
                        cATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGa                                                                                                                        <  1:201703/149‑1 (MQ=255)
                          ttCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAg                                                                                                                       >  8:176308/1‑148 (MQ=255)
                                ccGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGATTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGattatt                                                                                                                <  7:240086/149‑1 (MQ=255)
                                  gAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTAc                                                                                                              <  1:168982/149‑1 (MQ=255)
                                  gAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTAc                                                                                                              <  3:86240/149‑1 (MQ=255)
                                    aCTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCACTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCTGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACtt                                                                                                            <  1:238875/149‑1 (MQ=255)
                                     cTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTg                                                                                                           >  1:259574/1‑149 (MQ=255)
                                             aGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCggctgg                                                                                                    <  5:132273/148‑1 (MQ=255)
                                               gTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCtt                                                                                                 >  2:45542/1‑149 (MQ=255)
                                                             cccAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGc                                                                                   <  6:119376/149‑1 (MQ=255)
                                                              ccAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACg                                                                                    >  5:245410/1‑147 (MQ=255)
                                                              ccAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACg                                                                                    >  1:258106/1‑147 (MQ=255)
                                                                aGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAgg                                                                                                                      <  2:151628/111‑1 (MQ=255)
                                                                       cccTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGgcgcg                                                                           <  8:128782/147‑1 (MQ=255)
                                                                              tGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGc                                                                  <  7:176308/149‑1 (MQ=255)
                                                                              tGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGc                                                                  >  4:139944/1‑149 (MQ=255)
                                                                                cAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCACGGCCCGCTCGTTGGCtt                                                                >  5:219568/1‑149 (MQ=255)
                                                                                 aGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTc                                                               <  2:259574/149‑1 (MQ=255)
                                                                                 aGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTc                                                               >  5:264727/1‑149 (MQ=255)
                                                                                 aGGCCTTTTCGTACCTGCTGACCGAGTTAGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGTCGCGCTCGTTGGCTTc                                                               >  6:183477/1‑149 (MQ=255)
                                                                                   gCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGt                                                             >  1:316255/1‑149 (MQ=255)
                                                                                   gCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGt                                                             >  5:191784/1‑149 (MQ=255)
                                                                                              ccTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTctgc                                                  >  2:51327/1‑149 (MQ=255)
                                                                                              ccTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTctgc                                                  >  3:182809/1‑149 (MQ=255)
                                                                                                       cgcgTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTc                                         <  6:133543/149‑1 (MQ=255)
                                                                                                       cgcgTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTc                                         >  5:133543/1‑149 (MQ=255)
                                                                                                               ccGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCgg                                  <  3:139944/148‑1 (MQ=255)
                                                                                                               ccGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCgg                                  <  1:45542/148‑1 (MQ=255)
                                                                                                                cGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCt                                <  2:316255/149‑1 (MQ=255)
                                                                                                                cGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTACGTCAGCGGTCTGCTGAGACTTCTGAGCGGCt                                >  6:227052/1‑149 (MQ=255)
                                                                                                                 ggTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAg                                                                                          <  3:71245/90‑1 (MQ=255)
                                                                                                                 ggTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAg                                                                                          <  3:75592/90‑1 (MQ=255)
                                                                                                                 ggTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAg                                                                                          >  4:71245/1‑90 (MQ=255)
                                                                                                                 ggTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAg                                                                                          >  4:75592/1‑90 (MQ=255)
                                                                                                                   tttCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCt                                                                 >  8:53257/1‑113 (MQ=255)
                                                                                                                   tttCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCt                                                                 <  7:53257/113‑1 (MQ=255)
                                                                                                                   tttCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGCGCCTTCTGAGCGGCTGc                              >  2:100257/1‑148 (MQ=255)
                                                                                                                       aaaGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTg                         >  2:219081/1‑149 (MQ=255)
                                                                                                                         aGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGcg                       >  2:131392/1‑149 (MQ=255)
                                                                                                                             aaTAAAACAGCCTCCCCACACGTGCGTCGGCTTCTTAACCATCCGTGAGTATTATTACTTGCGGCTGGCGTTGTCCAGCATACGCCGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGcgtcgt                   <  4:182809/149‑1 (MQ=255)
                                                                                                                               ttaagaaGCCGACCCACTTGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGcgtcgt                   >  6:308949/6‑147 (MQ=255)
                                                                                                                               tAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCa                 >  8:235598/1‑149 (MQ=255)
                                                                                                                                aaaaaaGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAg                >  1:248546/1‑149 (MQ=255)
                                                                                                                                 aaaaaGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAGc               <  1:219081/149‑1 (MQ=255)
                                                                                                                                  agaaGCCGACCCACTTGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAGcc              >  3:217014/3‑149 (MQ=255)
                                                                                                                                  aaaaGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAGcc              >  7:234756/1‑149 (MQ=255)
                                                                                                                                       ccGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAGCCTTACg         <  7:275955/149‑1 (MQ=255)
                                                                                                                                            ccACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTACTGAGCGGCTGCCATTGCGTCGTCAGCCTGAAGGTAgg    >  4:68648/1‑149 (MQ=255)
                                                                                                                                             cacaAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAGCCTTACGGTAg     <  6:314456/147‑1 (MQ=255)
                                                                                                                                             cacaAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAGCCTTACGGTAg     >  5:314456/1‑147 (MQ=255)
                                                                                                                                             cacaAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAGCCTTACGGTAGGc   <  2:248546/149‑1 (MQ=255)
                                                                                                                                              acaAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAGCCTTACGGTAGGCt  >  5:297305/1‑149 (MQ=255)
                                                                                                                                              acaAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAGCCTTACGGTAGGCt  <  6:191784/149‑1 (MQ=255)
                                                                                                                                               caAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGcgt                      >  7:293876/1‑128 (MQ=255)
                                                                                                                                               caAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGcgt                      <  8:293876/128‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GGGGGCCGCACCGATGGTGACGAGCATTCCGCCCGAACTGCAGTAAGGTAATTTCGGCGTACCCAGGCCTGCCCTTAGTGCAGGCCTTTTCGTACCTGCTGACCGCGTTTGCCGGTTTCAAAGGCAATAAAAAAGCCGACCCACAAGTGGGTCGGCTTCTTAACAATCCGTGAGGATTATTACTTGCGGCTGGCTTTGTCCAGCATACGCAGGGCGCGCTCGTTGGCTTCGTCAGCGGTCTGCTGAGCCTTCTGAGCGGCTGCCATTGCGTCGTCAGCCTTACGGTAGGCT  >  NC_002947/2650016‑2650306

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