Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,707,330:1 +A 8.0% intergenic (+106/+36) PP_1501 → / ← PP_1502 hypothetical protein/OmpA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,707,3301.A 8.0% 338.8 / 20.8 88intergenic (+106/+36)PP_1501/PP_1502hypothetical protein/OmpA family protein
Reads supporting (aligned to +/- strand):  ref base . (37/44);  new base A (3/4);  total (40/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGTCCG  >  NC_002947/1707181‑1707475
                                                                                                                                                       |                                                                                                                                                 
cgGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGACTTTGTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCcgcg                                                                                                                                                    >  2:2133940/1‑150 (MQ=255)
          cATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCgg                                                                                                                                         <  2:814252/150‑1 (MQ=255)
             gggCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGATCGACGCGGGGc                                                                                                                                      <  2:887075/150‑1 (MQ=255)
               gCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCtt                                                                                                                                    <  2:2123189/150‑1 (MQ=255)
                 gAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGt                                                                                                                                  >  1:1632301/1‑150 (MQ=255)
                 gAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGt                                                                                                                                  <  1:1926671/150‑1 (MQ=255)
                 gAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGt                                                                                                                                  >  1:794097/1‑150 (MQ=255)
                   ggAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTcc                                                                                                                                >  2:2100269/1‑150 (MQ=255)
                    gAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCCGAATTTGTGGGGCAGGCCAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGATCGACGCGGGGCttttctt                                                                                                                               <  1:2341869/150‑6 (MQ=37)
                    gAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCt                                                                                                                               <  2:1781924/150‑1 (MQ=255)
                     aaaCACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGATCGACGCGGGGCttttcttt                                                                                                                              <  2:2109385/150‑7 (MQ=38)
                       acacCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTc                                                                                                                            <  1:1111422/150‑1 (MQ=255)
                        cacCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGCGCTTGCCCTTCCa                                                                                                                           >  2:1111422/1‑150 (MQ=255)
                         acCGCGGCAGGTGCAGTGGGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAg                                                                                                                          >  2:1474271/1‑150 (MQ=255)
                         acCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAg                                                                                                                          >  1:549862/1‑150 (MQ=255)
                         acCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAg                                                                                                                          >  1:105257/1‑150 (MQ=255)
                          ccGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGc                                                                                                                         >  2:268288/1‑150 (MQ=255)
                                cAGGTGCAGTGAGGTCTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTagtca                                                                                                                   <  2:1678677/150‑1 (MQ=255)
                                  ggTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTagccagc                                                                                                                 <  2:626895/150‑1 (MQ=255)
                                  ggTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTagccagc                                                                                                                 <  1:2100269/150‑1 (MQ=255)
                                  ggTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTagccagc                                                                                                                 <  2:533697/150‑1 (MQ=255)
                                   gTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCATAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTagccagcc                                                                                                                <  1:581026/150‑1 (MQ=255)
                                   gTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTc                                                                                                                            >  1:141769/1‑138 (MQ=255)
                                   gTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTc                                                                                                                            <  2:141769/138‑1 (MQ=255)
                                        gTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCgggg                                                                                                                                       <  2:36328/122‑1 (MQ=255)
                                        gTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCgggg                                                                                                                                       >  1:36328/1‑122 (MQ=255)
                                           aGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCAGCAAGTGGCTTTGTCGGTCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCt                                                                                                        <  2:1872736/150‑1 (MQ=255)
                                                   gcgcAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTgg                                                                                                <  2:1632301/150‑1 (MQ=255)
                                                   gcgcAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTgg                                                                                                <  1:300592/150‑1 (MQ=255)
                                                   gcgcAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTgg                                                                                                <  2:549862/150‑1 (MQ=255)
                                                   gcgcAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTgg                                                                                                <  1:268288/150‑1 (MQ=255)
                                                     gcagcaTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCt                                                                                              >  1:1097561/1‑150 (MQ=255)
                                                         caTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTagcca                                                                                                                   >  2:1334652/1‑125 (MQ=255)
                                                         caTTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTagcca                                                                                                                   <  1:1334652/125‑1 (MQ=255)
                                                                  cgcCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTGAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCgtggt                                                                                 <  2:1748961/150‑1 (MQ=255)
                                                                     ccGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTagccagc                                                                                                                 >  1:742809/1‑115 (MQ=255)
                                                                     ccGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTagccagc                                                                                                                 <  2:742809/115‑1 (MQ=255)
                                                                       gtgagaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCAc                                                                            <  2:2214376/148‑1 (MQ=255)
                                                                               gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGt                                                                    >  1:1025658/1‑150 (MQ=255)
                                                                                aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCt                                                                                                                               <  2:1005346/90‑1 (MQ=255)
                                                                                aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCt                                                                                                                               >  1:1005346/1‑90 (MQ=255)
                                                                                aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGtt                                                                   <  1:784344/150‑1 (MQ=255)
                                                                                 gcgcCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGtt                                                                   >  1:1829964/1‑149 (MQ=255)
                                                                                 gcgcCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGtt                                                                   <  2:1829964/149‑1 (MQ=255)
                                                                                  cgccgcGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCAc                                                                                           <  1:120181/124‑1 (MQ=255)
                                                                                  cgccgcGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCAc                                                                                           >  2:120181/1‑124 (MQ=255)
                                                                                  cgccgcGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGt                                                                 <  2:1018594/150‑1 (MQ=255)
                                                                                       cgcgGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTc                                                                                                                            >  2:664468/1‑86 (MQ=255)
                                                                                       cgcgGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTc                                                                                                                            <  1:664468/86‑1 (MQ=255)
                                                                                        gcgggcggCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTgg                                                           <  1:105578/150‑1 (MQ=255)
                                                                                              cggcgCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATg                                                      <  2:1097561/150‑1 (MQ=255)
                                                                                                 cgcTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCCCCGGCTGGTTGTCGTTGGAGATGGTc                                                   <  1:2133940/150‑1 (MQ=255)
                                                                                                 cgcTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTc                                                   <  2:1025658/150‑1 (MQ=255)
                                                                                                  gcTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCCCCGCGCGTGGTGCCATCGGCTGGTTGTGGTTGGAGAtcggaa                                                  >  2:1054045/1‑145 (MQ=255)
                                                                                                      gATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTg                                                            >  1:1291764/1‑136 (MQ=255)
                                                                                                      gATCTCACAGGCGGCACAACACTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTg                                                            <  2:1291764/136‑1 (MQ=255)
                                                                                                                   gcACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTc                                 >  2:1358160/1‑150 (MQ=255)
                                                                                                                   gcACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTc                                 >  1:912601/1‑150 (MQ=255)
                                                                                                                    cacaACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCt                                >  1:1171609/1‑150 (MQ=255)
                                                                                                                     acaacaTTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTg                               <  1:1315069/150‑1 (MQ=255)
                                                                                                                       aacaTTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCt                                               <  1:1007141/132‑1 (MQ=255)
                                                                                                                       aacaTTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCt                                               >  2:1007141/1‑132 (MQ=255)
                                                                                                                          aTTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGtt                                                             <  1:2169207/115‑1 (MQ=255)
                                                                                                                          aTTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGtt                                                             >  2:2169207/1‑115 (MQ=255)
                                                                                                                           ttATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTg                                                                                                 >  2:216562/1‑78 (MQ=255)
                                                                                                                           ttATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTg                                                                                                 <  1:216562/78‑1 (MQ=255)
                                                                                                                            tATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGt                        <  2:173899/150‑1 (MQ=255)
                                                                                                                            tATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCACTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGt                        >  1:968348/1‑150 (MQ=255)
                                                                                                                               cGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTccc                     <  1:1474271/150‑1 (MQ=255)
                                                                                                                                 gCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGa                                       >  1:900512/1‑130 (MQ=255)
                                                                                                                                 gCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGa                                       <  2:900512/130‑1 (MQ=255)
                                                                                                                                   gAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGt                 >  1:848922/1‑150 (MQ=255)
                                                                                                                                      cAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGACTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGc              >  1:1095084/1‑150 (MQ=255)
                                                                                                                                       aaaGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTgc                              <  1:240512/133‑1 (MQ=255)
                                                                                                                                       aaaGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTgc                              >  2:240512/1‑133 (MQ=255)
                                                                                                                                       aaaGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCt             >  2:497605/1‑150 (MQ=255)
                                                                                                                                       aaaGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCAACCGGCGGTTCTGCGCACGTCCCGAGTTGCt             >  2:796638/1‑150 (MQ=255)
                                                                                                                                          gAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGa          <  1:1358160/150‑1 (MQ=255)
                                                                                                                                           aaaaGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAg         >  1:1395132/1‑150 (MQ=255)
                                                                                                                                             aaGCCCCGCGATCGACGCGGGGCTTGTCCTTTCAGCTAGCGAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGga                                                          >  1:985465/1‑99 (MQ=39)
                                                                                                                                             aaGCCCCGCGATCGACGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGga                                                          <  2:985465/99‑1 (MQ=39)
                                                                                                                                             aaGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCAc                                                 <  1:2309066/108‑1 (MQ=255)
                                                                                                                                             aaGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCAc                                                 >  2:2309066/1‑108 (MQ=255)
                                                                                                                                               gCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATg                                                      >  2:579441/1‑101 (MQ=255)
                                                                                                                                               gCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATg                                                      <  1:579441/101‑1 (MQ=255)
                                                                                                                                                ccccGCGATCGACGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGTCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGTc    >  1:234795/1‑150 (MQ=255)
                                                                                                                                                ccccGCGATCGACGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGTc    >  2:1705803/1‑150 (MQ=255)
                                                                                                                                                ccccGCGCCGATCGCCGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAGCCTCGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTGGACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGTc    <  2:794097/150‑1 (MQ=255)
                                                                                                                                                  ccGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCATACGCGCTTAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGCCGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGTCCg  >  1:422942/1‑150 (MQ=255)
                                                                                                                                                   cgcgTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGTGCACGcccc                     <  1:1175468/130‑1 (MQ=255)
                                                                                                                                                   cgcgTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGTGCACGcccc                     >  2:1175468/1‑130 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
CGGTGCGCTGCATGGGCGAGGAAACACCGCGGCAGGTGCAGTGAGGTTTCGGCGCAGCATTTTTGGCGCCCGCGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGGCACAACATTATCGGCGAACAAAGAAAAGCCCCGCGTCGATCGCGGGGCTTGTCCTTTCAGCTAGCCAGCCTTAGCGCTGAACCTGGCTCACCGAGCGTGGTGCCACCGGCTGGTTGTCGTTGGAGATGGTCACCTCCACCCGACGGTTCTGCGCACGTCCCGAGTTGCTCGAGTTGTCCG  >  NC_002947/1707181‑1707475

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