Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,415,076 T→C 7.2% intergenic (+141/+27) PP_5644 → / ← PP_3899 membrane protein/hicB family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,415,0760TC 7.2% 295.9 / 9.0 97intergenic (+141/+27)PP_5644/PP_3899membrane protein/hicB family protein
Reads supporting (aligned to +/- strand):  ref base T (48/42);  new base C (4/3);  total (52/45)
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.44e-01

AGTAGATGGCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCT  >  NC_002947/4414931‑4415216
                                                                                                                                                 |                                                                                                                                            
aGTAGATGGCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCa                                                                                                                                          <  2:14932/150‑1 (MQ=255)
aGTAGATGGCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCa                                                                                                                                          <  2:2368229/150‑1 (MQ=255)
        gCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGc                                                                                                                                  >  2:991312/1‑150 (MQ=255)
          tATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCtt                                                                                                                                >  1:377905/1‑150 (MQ=255)
              aaGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCt                                                                                                                            <  1:1138549/150‑1 (MQ=255)
                gCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCttt                                                                                                                          <  2:578159/150‑1 (MQ=255)
                 cACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCtttt                                                                                                                         >  2:856860/1‑149 (MQ=255)
                   ctcctcTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTggg                                                                                                                       >  2:724390/1‑150 (MQ=255)
                   ctcctcTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTGTCACGCCCTCCCCTGCACAGTGAAGGCTTATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTggg                                                                                                                       <  1:157589/150‑1 (MQ=255)
                    tcctcTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGggc                                                                                                                      >  1:2215614/1‑150 (MQ=255)
                     cctcTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGggcg                                                                                                                     >  1:1073382/1‑150 (MQ=255)
                      ctctTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTggg                                                                                                                       <  2:393133/147‑1 (MQ=255)
                      ctctTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTggg                                                                                                                       >  1:393133/1‑147 (MQ=255)
                      ctctTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGggcgg                                                                                                                    >  2:2013844/1‑150 (MQ=255)
                      ctctTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGggcgg                                                                                                                    >  2:1550859/1‑150 (MQ=255)
                        ctTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCt                                                                                                                  >  1:1996382/1‑150 (MQ=255)
                        ctTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCACTGCAATGCGGGc                                                                                                                                  >  1:1337143/1‑134 (MQ=255)
                        ctTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCACTGCAATGCGGGc                                                                                                                                  <  2:1337143/134‑1 (MQ=255)
                              gccgccATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGt                                                                                                            >  1:157853/1‑150 (MQ=255)
                               ccgccATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGtg                                                                                                           <  2:2187925/150‑1 (MQ=255)
                                cgccATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGtga                                                                                                          <  1:1550859/150‑1 (MQ=255)
                                 gccATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGtgat                                                                                                         <  2:1073382/150‑1 (MQ=255)
                                   cATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGtgatga                                                                                                       <  2:157853/150‑1 (MQ=255)
                                    aTCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCACTGCAATGCGGGCtaa                                                                                                                               >  1:979081/1‑123 (MQ=255)
                                    aTCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCACTGCAATGCGGGCtaa                                                                                                                               <  2:979081/125‑3 (MQ=255)
                                          tGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGt                                                                                                                                        <  1:1864745/110‑1 (MQ=255)
                                          tGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGt                                                                                                                                        >  2:1864745/1‑110 (MQ=255)
                                              cccACAGTCCGCCTCGGCGGGCTTCACCATGCATTAACGATTTCTTCAAGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTATGGGCGGCGACGCGTGATGACCCAGGTACTc                                                                                            >  1:2223999/1‑150 (MQ=255)
                                               ccACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGt                                                                                                            >  2:1093560/1‑133 (MQ=255)
                                               ccACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGt                                                                                                            <  1:1093560/133‑1 (MQ=255)
                                                cacaGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCa                                                                                                                                          <  2:938506/102‑1 (MQ=255)
                                                cacaGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCa                                                                                                                                          <  2:938300/102‑1 (MQ=255)
                                                cacaGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCa                                                                                                                                          >  1:938506/1‑102 (MQ=255)
                                                cacaGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCa                                                                                                                                          >  1:938300/1‑102 (MQ=255)
                                                cacaGTCCGCCTCGGCGGGCTTCAACGCATATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCAc                                                                                          >  1:1270714/1‑150 (MQ=255)
                                                     tCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTg                                                                                                                                       >  1:231714/1‑100 (MQ=255)
                                                     tCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTg                                                                                                                                       <  2:231714/100‑1 (MQ=255)
                                                     tCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTTGGCGGCTACGCGTGATGACCCAGGTACTCACTAAgc                                                                                     >  2:212020/1‑150 (MQ=255)
                                                        gCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCAGGCCCTTAACTGCAAAATGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTACGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCt                                                                                  >  1:2102283/1‑150 (MQ=255)
                                                        gCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCa                                                                                                   <  1:2028796/133‑1 (MQ=255)
                                                        gCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCa                                                                                                   >  2:2028796/1‑133 (MQ=255)
                                                         ccTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCtt                                                                                 >  2:1615450/1‑150 (MQ=255)
                                                           tCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTg                                                                               >  2:1927299/1‑150 (MQ=255)
                                                            cggcggGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGc                                                                              <  2:2223999/150‑1 (MQ=255)
                                                            cggcggGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGc                                                                              <  2:377905/150‑1 (MQ=255)
                                                              gcggGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCgt                                                                            <  1:724390/150‑1 (MQ=255)
                                                                gggCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGATCCAGGTACTCACTAAGCGCTTTGCgtgt                                                                          <  1:665838/150‑1 (MQ=255)
                                                                gggCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCa                                                                                           <  1:2396336/133‑1 (MQ=255)
                                                                gggCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCa                                                                                           >  2:2396336/1‑133 (MQ=255)
                                                                      cACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATc                                                                    <  2:2215614/150‑1 (MQ=255)
                                                                       aCCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCa                                                                   >  2:1877846/1‑150 (MQ=255)
                                                                         cATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCAtt                                                                 <  1:2013844/150‑1 (MQ=255)
                                                                          aTGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTc                                                                >  1:1435702/1‑150 (MQ=255)
                                                                            gCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCaa                                                              <  1:1877846/150‑1 (MQ=255)
                                                                             taTTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAg                                                             <  2:2102283/149‑1 (MQ=255)
                                                                             cATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAg                                                             <  1:1927299/150‑1 (MQ=255)
                                                                              aTTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCtttt                                                                                                                              <  1:2283729/84‑1 (MQ=255)
                                                                              aTTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCtttt                                                                                                                              >  2:2283729/1‑84 (MQ=255)
                                                                              aTTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGTGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCgt                                                                            >  2:115241/1‑134 (MQ=255)
                                                                              aTTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGTGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCgt                                                                            <  1:115241/134‑1 (MQ=255)
                                                                                  aCGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGa                                                           <  1:1783488/147‑1 (MQ=255)
                                                                                  aCGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGa                                                           >  2:1783488/1‑147 (MQ=255)
                                                                                      ttttttCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTACCTATTGCTGTACGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACtt                                                    <  2:1270714/150‑1 (MQ=255)
                                                                                       tttttCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACttt                                                   <  2:1435702/150‑1 (MQ=255)
                                                                                        ttttCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTg                                                  <  2:223692/150‑1 (MQ=255)
                                                                                              cGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTc                                                                                                                             >  1:2018369/1‑69 (MQ=255)
                                                                                              cGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTc                                                                                                                             <  2:2018369/69‑1 (MQ=255)
                                                                                               gCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTtggtg                                           >  2:1492584/1‑150 (MQ=255)
                                                                                                  cTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCt                                        >  1:1708470/1‑150 (MQ=255)
                                                                                                     cccTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGcc                                     >  1:788316/1‑150 (MQ=255)
                                                                                                      ccTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCa                                    >  1:199376/1‑150 (MQ=255)
                                                                                                       cTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTGTGCCAGAGGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCAc                                   >  2:1891381/1‑150 (MQ=255)
                                                                                                         tcacaGTGAAGCCGCATCCGCTTCCCTACGTTAGCCAGCCTTGCATTGCGGCCTTTTGTTGGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATGATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGc                                 >  2:1289592/2‑150 (MQ=255)
                                                                                                          cacaGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACTCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCt                                <  1:991312/150‑1 (MQ=255)
                                                                                                             aGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCAAGCACGCTGGc                             >  1:1575858/1‑150 (MQ=255)
                                                                                                               tGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACttt                                                   >  1:1672126/1‑126 (MQ=255)
                                                                                                               tGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACttt                                                   <  2:1672126/126‑1 (MQ=255)
                                                                                                               tGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCg                           >  1:1525070/1‑150 (MQ=255)
                                                                                                               tGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCg                           >  1:1524898/1‑150 (MQ=255)
                                                                                                                   ggCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAg                       <  1:2204224/150‑1 (MQ=255)
                                                                                                                     cTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCAGAGACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGt                     >  2:1417998/1‑150 (MQ=255)
                                                                                                                       cATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTc                                                                                            >  2:2122905/1‑77 (MQ=255)
                                                                                                                       cATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTc                                                                                            <  1:2122905/77‑1 (MQ=255)
                                                                                                                        aTCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCa                  >  1:208210/1‑150 (MQ=255)
                                                                                                                           cGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCtt                                                                                 <  1:1705258/84‑1 (MQ=255)
                                                                                                                           cGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCtt                                                                                 >  2:1705258/1‑84 (MQ=255)
                                                                                                                           cGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGt               >  2:724786/1‑150 (MQ=255)
                                                                                                                                    aCGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCCGATCATTCAAg                                                             >  1:1956648/1‑95 (MQ=255)
                                                                                                                                    aCGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAg                                                             <  2:1956648/95‑1 (MQ=255)
                                                                                                                                        tAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCt  <  2:1708470/150‑1 (MQ=255)
                                                                                                                                        aaGCCCGCACTGCAATGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCt  >  2:378646/2‑150 (MQ=255)
                                                                                                                                           cccGCACTGCAATGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAg                 <  2:305147/132‑1 (MQ=255)
                                                                                                                                           cccGCACTGCAATGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAg                 >  1:305147/1‑132 (MQ=255)
                                                                                                                                            ccGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTc              <  1:1608669/134‑1 (MQ=255)
                                                                                                                                            ccGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTc              >  2:1608669/1‑134 (MQ=255)
                                                                                                                                             cGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCt                                                <  2:212648/99‑1 (MQ=255)
                                                                                                                                             cGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCt                                                >  1:212648/1‑99 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
AGTAGATGGCTATGAAGCACTCCTCTTGAAGCCGCCATCAAGTGATCCCACAGTCCGCCTCGGCGGGCTTCACCATGCATTCACGCTTTTTTCACGCCCTCCCCTGCACAGTGAAGGCTCATCCGCTTCCCTACGTTAGCCCGCATTGCAGTGCGGGCTTTTCTTTGGGCGGCTACGCGTGATGACCCAGGTACTCACTAAGCGCTTTGCGTGTCAGATCATTCAAGGAAATACTTTGCTTGGTGGCTGCCACGCTGGCCGCTAGGTGCAGGTCGTGCCCTACCCT  >  NC_002947/4414931‑4415216

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: