Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 242,088 A→C 5.5% intergenic (+20/‑203) PP_0189 → / → PP_0190 heme biosynthesis protein/disulfide bond formation protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947242,0880AC 5.5% 223.9 / 2.2 72intergenic (+20/‑203)PP_0189/PP_0190heme biosynthesis protein/disulfide bond formation protein B
Reads supporting (aligned to +/- strand):  ref base A (33/35);  new base C (2/2);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

AGCTGGCCCGCCTGCTCGCGGGGCTGGGCGAGACCGAGCGCAGCAACCAGCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGATTCAGTGACTTGTCGTTGCTGTCGTGCCTTGAAACAAAGC  >  NC_002947/241944‑242234
                                                                                                                                                |                                                                                                                                                  
aGCTGGCCCGCCTGCTCGCGGGGCTGGGCGAGACCGAGCGCAGCAACCAGCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCaaa                                                                                                                                                  <  5:162632/147‑1 (MQ=255)
aGCTGGCCCGCCTGCTCGCGGGGCTGGGCGAGACCGAGCGCAGCAACCAGCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCaaa                                                                                                                                                  <  8:237330/147‑1 (MQ=255)
       ccGCCTGCTCGCGGGGCTGGGCGAGACCGAGCGCAGCAACCAGCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGc                                                                                                                                         >  5:69667/1‑149 (MQ=255)
           cTGCTCGCGTGGCTGGGCGAGACCGAGCGCAGCAACCAGCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTg                                                                                                                                     <  3:83578/149‑1 (MQ=255)
           cTGCTCGCGGGGCTGGGCGAGACCGAGCGCAGCAACCAGCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTg                                                                                                                                     <  6:75101/149‑1 (MQ=255)
              cTCGCGGGGCTGGGCGAGACCGAGCGCAGCAACCAGCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTgggg                                                                                                                                  >  3:226351/1‑149 (MQ=255)
                        tGGGCGAGACCGAGCGCAGCAACCAGCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAAtt                                                                                                                        >  5:293534/1‑149 (MQ=255)
                        tGGGCGAGACCGAGCGCAGCAACCAGCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAAtt                                                                                                                        >  5:296112/1‑149 (MQ=255)
                          ggCGAGACCGAGCGCAGCAACCAGCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTa                                                                                                                      >  4:166422/1‑149 (MQ=255)
                            cGAGACCGAGCGCAGCAACCAGCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCACCTTTGCGCGGGCCTTggggtaggt                                                                                                                             <  1:289031/140‑6 (MQ=38)
                            cGAGACCGAGCGCAGCAACCAGCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCACCTTTGCGCGGGCCTTggggtaggt                                                                                                                             >  2:289031/1‑135 (MQ=38)
                             gagaCCGAGCGCAGCAACCAGCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCt                                                                                                                   >  4:296430/1‑149 (MQ=255)
                             gagaCCGAGCGCAGCAACCAGCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCt                                                                                                                   <  6:190090/149‑1 (MQ=255)
                             gagaCCGAGCGCAGCAACCAGCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCt                                                                                                                   >  2:136502/1‑149 (MQ=255)
                              agaCCGAGCGCAGCAACCAGCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCtt                                                                                                                  <  3:296430/149‑1 (MQ=255)
                                 ccGAGCGCAGCAACCAGCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCg                                                                                                               <  8:242827/149‑1 (MQ=255)
                                    aGCGCAGCAACCAGCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTg                                                                                                            >  2:996/1‑149 (MQ=255)
                                       gcagcaACCAGCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCa                                                                                                         <  6:217567/149‑1 (MQ=255)
                                        cagcaACCAGCTGTTCCAGGAGGGCCTTGGCCTTGTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAAGCCCGCACCTTTGCGCGGGCCTTggggtaggtc                                                                                                                            <  7:43299/129‑7 (MQ=25)
                                        cagaaACCACCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCCCATGATGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCACAAAGCCCGCACCAAACCGCGGGCCATggggtaggtc                                                                                                                            >  8:43299/1‑123 (MQ=255)
                                          gcaACCAGCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCg                                                                                                               >  4:76000/1‑140 (MQ=255)
                                             aCCAGCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGa                                                                                                   <  6:69667/149‑1 (MQ=255)
                                               cAGGTGTTCCAGGAGGGCCTTGGCATGCTGGACGATCGCCTGCTGGCGCTGCCTTTGCCAGTAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATc                                                                                                 <  7:125693/149‑1 (MQ=255)
                                                 gCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCt                                                                                                                   >  8:275163/1‑129 (MQ=255)
                                                 gCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCt                                                                                                                   <  7:275163/129‑1 (MQ=255)
                                                   tGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGaa                                                                                             <  3:61886/149‑1 (MQ=255)
                                                     ttCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGaa                                                                                             <  4:226351/147‑1 (MQ=255)
                                                     ttCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGt                                                                                           >  2:204469/1‑149 (MQ=255)
                                                     ttCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGt                                                                                           >  4:183899/1‑149 (MQ=255)
                                                        cAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCa                                                                                        <  6:296112/149‑1 (MQ=255)
                                                               gCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATttcttc                                                                                 >  5:197874/1‑149 (MQ=255)
                                                                ccTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTCATCTGAAGTCCATTTCTTcc                                                                                >  7:243173/1‑149 (MQ=255)
                                                                    ggCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCg                                                                            >  5:76381/1‑149 (MQ=255)
                                                                     gCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGc                                                                           <  6:76381/149‑1 (MQ=255)
                                                                       ctgctgGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGcct                                                                         >  8:234066/1‑149 (MQ=255)
                                                                          ctgGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGcctc                                                                        <  1:996/147‑1 (MQ=255)
                                                                           tgGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTa                                                                     >  8:141538/1‑149 (MQ=255)
                                                                            ggACGACCGCCTGCTGGCGCTGCCTTTGCCAGAACGCGTGCGCGCCTGACCTAACCCAAGCCCCGCGCAAAGGAGCGCGCTTTGGGTGCAGGTAATTTATCTATCGCTGTCAGTTTCATATCAAGTCCATTTCTTCCGTCGCCTCCTAc                                                                    >  6:194001/1‑149 (MQ=255)
                                                                            gcaCGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGACTGACCTACCCCAAGCCCCGCGCAAAGGTGCGGGCTTTGCGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTAc                                                                    >  4:118164/3‑149 (MQ=255)
                                                                             gACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTAc                                                                    <  6:197874/148‑1 (MQ=255)
                                                                              aCGAGCGCCTGCTGGCGCTGCCTTTGCCTGCAGGCGTGCGAGCCTGAGCTACCCCAAGGCGCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACgg                                                                  <  8:243173/149‑1 (MQ=255)
                                                                                gAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCa                                                                >  1:236500/1‑149 (MQ=255)
                                                                                      ctgctgGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTg                                                          >  8:183285/1‑149 (MQ=255)
                                                                                               cTGCCTTTGGCAGATGGCGTGCTCGCCTGACCTAACGCAAGGCCCGGGGATATGTGCGGGCTTTGGGTGCAGGTAATTTATCTTTCTCTGTCAGTTTTATCTGATGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTcc                                                 <  5:194001/149‑1 (MQ=255)
                                                                                                  ccTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTAc                                              >  3:3245/1‑149 (MQ=255)
                                                                                                    tttGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTa                                                                                                                      >  8:146272/1‑75 (MQ=255)
                                                                                                    tttGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTa                                                                                                                      <  7:146272/75‑1 (MQ=255)
                                                                                                      tGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTAcgcg                                           <  3:166422/148‑1 (MQ=255)
                                                                                                        ccAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAACTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGAtt                                        <  5:125813/149‑1 (MQ=255)
                                                                                                         cAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGAtt                                        >  4:227553/1‑148 (MQ=255)
                                                                                                         cAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGAtt                                        >  3:289928/1‑148 (MQ=255)
                                                                                                          aGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCat                                                               <  4:130972/124‑1 (MQ=255)
                                                                                                          aGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGAtt                                        >  1:126889/1‑147 (MQ=255)
                                                                                                             aGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGATTCAGt                                    <  4:289928/148‑1 (MQ=255)
                                                                                                              ggCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGATTCAGTGa                                  <  2:236500/149‑1 (MQ=255)
                                                                                                                       gcCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGATTCAGTGACTTGTCGtt                         >  2:191022/1‑149 (MQ=255)
                                                                                                                         cTGACCTCCCCCAAGGCCCGCGCATAGGTGCTGGCTTTGGGGGCAGGTATTTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGATTCAGTGAc                                 <  8:108433/139‑1 (MQ=255)
                                                                                                                         cTGACCTACCCCAAGGCCAGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGATTCAGTGAc                                 >  7:108433/1‑139 (MQ=255)
                                                                                                                                 ccccAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGATTCAGTGACTTTTCGTTGCTGTCGTGc               <  3:118164/149‑1 (MQ=255)
                                                                                                                                 ccccAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGATTCAGTGACTTGTCGTTGCTGTCGTGc               >  8:85309/1‑149 (MQ=255)
                                                                                                                                  cccAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGATTCAGTGACTTGTCGTTGCTGTCGTGcc              <  7:234066/149‑1 (MQ=255)
                                                                                                                                   ccAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGATTCAGTGACTTGTCGTTGCTGTCGTGCCt             <  3:227553/149‑1 (MQ=255)
                                                                                                                                    cAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGATTCAGTGACTTGTCGTTGCTGTCGTGCCtt            >  6:102704/1‑149 (MQ=255)
                                                                                                                                      aGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGATTCAGTGACTTGTCGtt                         >  2:132317/1‑134 (MQ=255)
                                                                                                                                      aGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGATTCAGTGACTTGTCGtt                         <  1:132317/134‑1 (MQ=255)
                                                                                                                                      aGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGATTCAGTGACTTGTCGTTGCTGTCGTGCCTTGa          <  6:293534/149‑1 (MQ=255)
                                                                                                                                      aGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGATTCAGTGACTTGTCGTTGCTGTCGTGCCTTGa          <  3:183899/149‑1 (MQ=255)
                                                                                                                                      aGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGATTCAGTGACTTGTCGTTGCTGTCGTGCCTTGa          <  1:204469/149‑1 (MQ=255)
                                                                                                                                        gCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGATTCAGTGACTTGTCGTTGCTGTCGTGCCTTGaaa        >  2:188490/1‑149 (MQ=255)
                                                                                                                                        gCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGATTCAGTGACTTGTCGTTGCTGTCGTGCCTTGaaa        <  1:136502/149‑1 (MQ=255)
                                                                                                                                        gCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGATTCAGTGACTTGTCGTTGCTGTCGTGCCTTGaa         >  7:219815/1‑148 (MQ=255)
                                                                                                                                         cccGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGATTCAGTGACTTGTCGTTGCTGTCGTGCCTTGAAAc       <  2:148390/149‑1 (MQ=255)
                                                                                                                                            gcgcAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGATTCAGTGACTTGTCGTTGCTGTCGTGCCTTGAAACaaa    <  5:102704/149‑1 (MQ=255)
                                                                                                                                              gcAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGATTCAGTGACTTGTCGTTGCTGTCGTGCCTTGAAACAAAGc  <  7:85309/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
AGCTGGCCCGCCTGCTCGCGGGGCTGGGCGAGACCGAGCGCAGCAACCAGCTGTTCCAGGAGGGCCTTGGCCTGCTGGACGAGCGCCTGCTGGCGCTGCCTTTGCCAGAAGGCGTGCGCGCCTGACCTACCCCAAGGCCCGCGCAAAGGTGCGGGCTTTGGGGGCAGGTAATTTATCTTTCGCTGTCAGTTTGATCTGAAGTCCATTTCTTCCGTCGCCTCCTACGGCATATCTGGAAATTTCCTACGCGATTCAGTGACTTGTCGTTGCTGTCGTGCCTTGAAACAAAGC  >  NC_002947/241944‑242234

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: