Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,218,711 Δ1 bp 5.1% intergenic (+51/+32) oprJ → / ← nfxB outer membrane protein OprJ/HTH‑type transcriptional regulator NfxB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,218,7110G. 5.1% 259.0 / 7.3 79intergenic (+51/+32)oprJ/nfxBouter membrane protein OprJ/HTH‑type transcriptional regulator NfxB
Reads supporting (aligned to +/- strand):  ref base G (37/38);  new base . (2/2);  total (39/40)
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.79e-01

CGAAGCAGCCTTCATCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCATCGACGATGCCGTACACCATGCTCAG  >  NC_002947/3218567‑3218859
                                                                                                                                                 |                                                                                                                                                    
cGAAGCAGCCTTCATCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTg                                                                                                                                                 >  1:679303/1‑150 (MQ=255)
cGAAGCAGCCTTCATCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTg                                                                                                                                                 <  2:679303/150‑1 (MQ=255)
 gAAGCAGCCTTCATCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTgg                                                                                                                                                <  2:723356/150‑1 (MQ=255)
 gAAGCAGCCTTCATCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTgg                                                                                                                                                >  2:682012/1‑150 (MQ=255)
  aaGCAGCCTTCATCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGCATCGATCCCTggg                                                                                                                                               <  2:464186/150‑1 (MQ=255)
  aaGCAGCCTTCATCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCGCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTggg                                                                                                                                               <  2:239363/150‑1 (MQ=255)
      agcCTTCATCGAAGCCAGTACCCAGCGACAACTGGCGTTGGCCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGttt                                                                                                                                           >  2:823567/1‑150 (MQ=255)
          ttCATCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGttttttt                                                                                                                                       <  1:759489/150‑1 (MQ=255)
            cATCGAAGCCAGTACCCAGTGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTgg                                                                                                                                     <  1:612880/150‑1 (MQ=255)
            cATCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTgg                                                                                                                                     <  1:750782/150‑1 (MQ=255)
             aTCGAAGCCAGTCCCCAGCGACAACTGGCGGTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCGGGCCGAGCAGGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTGGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTggg                                                                                                                                    >  2:832804/1‑149 (MQ=255)
             aTCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGc                                                                                                                                    <  2:1234561/150‑1 (MQ=255)
               cGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTg                                                                                                                                  >  2:1160388/1‑150 (MQ=255)
                    ccAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGaa                                                                                                                             >  2:1248581/1‑150 (MQ=255)
                     cAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCACTGGGGTTTTTTTGGCTGGTGAAg                                                                                                                            >  2:57929/1‑150 (MQ=255)
                         aCCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAg                                                                                                                        <  1:375127/150‑1 (MQ=255)
                           ccAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGt                                                                                                                      >  1:758788/1‑150 (MQ=255)
                           ccAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGt                                                                                                                      <  1:682012/150‑1 (MQ=255)
                            cAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCa                                                                                                                         <  1:959030/146‑1 (MQ=255)
                            cAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCa                                                                                                                         >  2:959030/1‑146 (MQ=255)
                             aGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCa                                                                                                                    <  2:845267/150‑1 (MQ=255)
                                   aaCTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCagcgagc                                                                                                              <  2:993002/150‑1 (MQ=255)
                                     cTGGCGTTGGCCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCgg                                                                                                            <  1:823567/150‑1 (MQ=255)
                                        gCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTg                                                                                                         >  1:1002096/1‑150 (MQ=255)
                                         cGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGc                                                                                                        <  1:177153/150‑1 (MQ=255)
                                                 aaCCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCAc                                                                                                >  1:498475/1‑150 (MQ=255)
                                                   ccTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGtttt                                                                                                                                          <  2:273897/106‑1 (MQ=255)
                                                   ccTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGtttt                                                                                                                                          >  1:273897/1‑106 (MQ=255)
                                                    cTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGt                                                                                             >  1:260168/1‑150 (MQ=255)
                                                     tGTTCCGCGCACTCGTCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTc                                                                                            >  2:502779/1‑150 (MQ=255)
                                                     tGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAACCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGc                                                                                                                                                     <  2:949472/93‑1 (MQ=255)
                                                     tGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAACCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGc                                                                                                                                                     >  1:949472/1‑93 (MQ=255)
                                                       ttCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGa                                                                                          >  1:541011/1‑150 (MQ=255)
                                                       ttCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGa                                                                                          <  1:1160388/150‑1 (MQ=255)
                                                        tCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAg                                                                                         <  2:498475/150‑1 (MQ=255)
                                                          cgcgCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCg                                                                                                      >  2:356694/1‑135 (MQ=255)
                                                          cgcgCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCg                                                                                                      <  1:356694/135‑1 (MQ=255)
                                                              cACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAag                                                                                   <  1:502779/150‑1 (MQ=255)
                                                              cACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCATGGGTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAag                                                                                   <  1:654913/150‑1 (MQ=255)
                                                                 tCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTggg                                                                                                                                               <  1:734148/87‑1 (MQ=255)
                                                                 tCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTggg                                                                                                                                               >  2:734148/1‑87 (MQ=255)
                                                                          gcTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAagca                                                                                 <  1:753470/140‑1 (MQ=255)
                                                                          gcTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAagca                                                                                 >  2:753470/1‑140 (MQ=255)
                                                                                cGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGCATCGATCCCTggg                                                                                                                                               <  1:91354/72‑1 (MQ=38)
                                                                                cGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGCATCGATCCCTggg                                                                                                                                               >  2:91354/1‑72 (MQ=38)
                                                                                 gAGCAAGGGCTGACACCGCCTCCGTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACg                                                                >  2:1277695/1‑150 (MQ=255)
                                                                                   gCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGa                                                              <  1:1277695/150‑1 (MQ=255)
                                                                                     aaGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACg                                                            >  1:1212061/1‑150 (MQ=255)
                                                                                       gggCTGAAACAGTCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCAc                                                                 <  2:842764/143‑1 (MQ=255)
                                                                                       gggCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCAc                                                                 >  1:842764/1‑143 (MQ=255)
                                                                                        ggCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGa                                                                                          >  1:228790/1‑117 (MQ=255)
                                                                                        ggCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGa                                                                                          <  2:228790/117‑1 (MQ=255)
                                                                                          cTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTcgcg                                                       >  1:305037/1‑150 (MQ=255)
                                                                                             aaaCAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGa                                                                                          <  2:413342/112‑1 (MQ=255)
                                                                                             aaaCAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGa                                                                                          >  1:413342/1‑112 (MQ=255)
                                                                                               aCAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCgcgc                                                      >  2:1040568/1‑146 (MQ=255)
                                                                                               aCAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCCAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCgcgc                                                      <  1:1040568/146‑1 (MQ=255)
                                                                                                cAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGgcc                                                 <  2:1163017/150‑1 (MQ=255)
                                                                                                   cctcctTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGaaa                                                                                     <  2:436901/111‑1 (MQ=255)
                                                                                                   cctcctTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGaaa                                                                                     >  1:436901/1‑111 (MQ=255)
                                                                                                        tttGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGccgcccgccc                                         >  1:1086691/1‑150 (MQ=255)
                                                                                                        tttGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGccgcccgccc                                         >  2:829094/1‑150 (MQ=255)
                                                                                                             ttggCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGcgtc                                    <  1:344268/150‑1 (MQ=255)
                                                                                                                gCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGtt                                 >  2:348521/1‑150 (MQ=255)
                                                                                                                 cTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTc                                <  2:541011/150‑1 (MQ=255)
                                                                                                                  ttGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCg                               >  1:387314/1‑150 (MQ=255)
                                                                                                                  ttGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCg                               <  2:1212061/150‑1 (MQ=255)
                                                                                                                  ttGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCg                               >  2:1285095/1‑150 (MQ=255)
                                                                                                                     tCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCa                            <  1:1221929/150‑1 (MQ=255)
                                                                                                                       cAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAATGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGAAGCTCCAGCGTTTGCACGGACGTCGCGCTTGCCGCCCGCCCGCGTCGTTCGGCATc                          >  2:47677/1‑150 (MQ=255)
                                                                                                                       cAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCATc                          >  2:708659/1‑150 (MQ=255)
                                                                                                                         gCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCTAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGttt                                                                     <  2:1263468/105‑1 (MQ=255)
                                                                                                                         gCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCTAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGttt                                                                     >  1:1263468/1‑105 (MQ=255)
                                                                                                                            aCTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGCGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCATcgacga                     <  1:166748/150‑1 (MQ=255)
                                                                                                                                   cccAGGGCATCGATCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGGGCCGCAGCACCGAAGAAGGAGGGCTGGTCCAGAGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCATCGACGATGCCGTa              >  2:304856/1‑150 (MQ=255)
                                                                                                                                    ccAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCATCGACGATGCCGTac             >  1:306935/1‑150 (MQ=255)
                                                                                                                                        ggATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCATCGACGATGCCGTACACCa         >  1:177715/1‑150 (MQ=255)
                                                                                                                                              gATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCCCCGCGACGCAGCTGCATCTGGTGGCCAGGGTTTGCGCTGACGTAGCGATGGCCGCACGCCCGCGTCTTGCGGCATCGACGATGCCGTACACCATGCTc    <  1:1254752/150‑1 (MQ=255)
                                                                                                                                                tGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCATCGACGATGCCGTACACCATGCTCAg  <  2:758788/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CGAAGCAGCCTTCATCGAAGCCAGTACCCAGCGACAACTGGCGTTGGTCAACCTGTTCCGCGCACTCGGCGGCGGCTGGCCGAGCAAGGGCTGAAACAGCCTCCTTTGGTTGGCTTGTCCAGCCACTTGGCCCCAGGGATCGATGCCCTGGGGTTTTTTTGGCTGGTGAAGGCAGGTCAGCGAGCGGGTGCCGCAGCACCGTCGAGGAAAAGCAGCTCCAGCGTTTGCACGGACGTCGCGCTGGCCGCCCGCCCGCGTCGTTCGGCATCGACGATGCCGTACACCATGCTCAG  >  NC_002947/3218567‑3218859

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