Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,064,576 G→C 10.0% intergenic (+17/+17) ubiA → / ← PP_5319 4‑hydroxybenzoate octaprenyltransferase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,064,5760GC10.0% 157.8 / 8.2 60intergenic (+17/+17)ubiA/PP_53194‑hydroxybenzoate octaprenyltransferase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/26);  new base C (3/3);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

TTTGCCTGGGAATACTGGTCAACGCGCAAGCTGGACCGGGAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTGTCCTTGATGAAGGTATACAGCGGCTTGCCGTCGTAAGCCCACTGCATCGAACCGTCGTCACGC  >  NC_002947/6064434‑6064719
                                                                                                                                              |                                                                                                                                                
tttGCCTGGGAATACTGGTCAACGCGCAAGCTGGACCGGGAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGc                                                                                                                                            <  3:79554/149‑1 (MQ=255)
  tGCCTGGGAATACTGGTCAACGCGCAAGCTGGACCGGGAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGcc                                                                                                                                           <  5:87514/148‑1 (MQ=255)
    ccTGGGAATACTGGTCAACGCGCAAGCTGGACCGGGAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGcccc                                                                                                                                         <  4:108971/148‑1 (MQ=255)
     cTGGGAATACTGGTCAACGCGCAAGCTGGACCGGGAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTCGGGGCCGCTCTGCGGcccca                                                                                                                                        >  1:271752/1‑147 (MQ=255)
      tGGGAATACTGGTCAACGCGCAAGCTGGACCGGGAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGc                                                                                                                                      >  4:87253/1‑149 (MQ=255)
       gggAATTCTGGTCAACGCTCAAGGTTGACCGGGAATCGTTCTTCAAGGCATTTTTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTCGGGGCCGCTCTGCGGccccaaa                                                                                                                                      <  2:271752/148‑4 (MQ=255)
            tACTGGTCAACGCGCAAGCTGGACCGGGAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCTCTGCGGccccaaac                                                                                                                                     <  3:193014/144‑5 (MQ=255)
            tACTGGTCAACGCGCAAGCTGGACCGGGAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCTCTGCGGccccaaac                                                                                                                                     >  4:193014/1‑140 (MQ=255)
             aCTGGTCAACGCGCAAGCTGGACCGGGAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCa                                                                                                                               >  2:187308/1‑149 (MQ=255)
               tGGTCAACGCGCAAGCTGGACCGGGAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATg                                                                                                                             >  1:83153/1‑149 (MQ=255)
                 gTCAACGCGCAAGCTGGACCGGGAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGAGGGGATGCTGGAGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGCGCGGCCCCGGCTATGACATGGc                                                                                                                           >  4:262203/1‑149 (MQ=255)
                    aaCGCGCAAGCTGGACCGGGAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTg                                                                                                                       <  1:231730/149‑1 (MQ=255)
                    aaCGCGCAAGCTGGACCGGGAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTGGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTg                                                                                                                       <  2:229951/149‑1 (MQ=255)
                     aCGCGCAAGCTGGACCGGGAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTgg                                                                                                                      <  1:187308/149‑1 (MQ=255)
                     aCGAGCAAGCTGGACGGGGAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGCGGGTTTGGTGGTGTTTATCGGAGTGGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGTCCGCAGAGCGGCCCCGGCTATCTCATGGCTTgg                                                                                                                      <  2:253409/149‑1 (MQ=255)
                       gcgcAAGCTGGACCGGGAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGc                                                                                                                     >  6:47811/1‑148 (MQ=255)
                             gCTGGACCGGGAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGc                                                                                                              <  3:87253/149‑1 (MQ=255)
                               tGGACCGGGAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCa                                                                                                            >  3:291217/1‑149 (MQ=255)
                                 gACGGGTAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGGGGGTATGGTGGTGTTTATCGGAGTTGTGCTGTATTACGGGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAg                                                                                                           <  4:150456/148‑1 (MQ=255)
                                 gACCGGGAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCa                                                                                                            <  2:65066/147‑1 (MQ=255)
                                  aCCGGGAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGAc                                                                                                         <  4:136438/149‑1 (MQ=255)
                                   ccGGGAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACg                                                                                                        >  7:259896/1‑149 (MQ=255)
                                   ccGGGAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACg                                                                                                        >  3:175271/1‑149 (MQ=255)
                                       gAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCtt                                                                                                    <  8:259896/149‑1 (MQ=255)
                                                ttCAAGGCATTTCTGCGCAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGt                                                                                           >  4:113832/1‑149 (MQ=255)
                                                    aGGCATTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCAcc                                                                                       <  4:175271/149‑1 (MQ=255)
                                                        aTTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGt                                                                                    >  3:208471/1‑148 (MQ=255)
                                                        aTTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGt                                                                                    >  7:232220/1‑148 (MQ=255)
                                                          ttCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCAt                                                                                 >  3:15777/1‑149 (MQ=255)
                                                             tGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGt                                                                               >  4:131665/1‑148 (MQ=255)
                                                               cacaACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCAc                                                                            >  2:299039/1‑149 (MQ=255)
                                                                acaGCCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCAcc                                                                           >  6:275025/1‑149 (MQ=255)
                                                                acaACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCAcc                                                                           <  1:299039/149‑1 (MQ=255)
                                                                  aaCCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCgg                                                                         >  2:287022/1‑149 (MQ=255)
                                                                          ggCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTg                                                                 >  4:290009/1‑149 (MQ=255)
                                                                              gggATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTGTCCt                                                             <  4:291217/149‑1 (MQ=255)
                                                                              gggATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTGTCCt                                                             <  3:113832/149‑1 (MQ=255)
                                                                              gggAGGCGGGGGTTTATCGGAGTTGTGCGGTATTACGCGTGGCGCTGAGTAGTTTGGGGCCGCCGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTGTCCt                                                             <  4:15777/149‑1 (MQ=255)
                                                                                 aTGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTGTCCTtga                                                          <  2:83153/149‑1 (MQ=255)
                                                                                 aTGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTGTCCTtga                                                          <  3:131665/149‑1 (MQ=255)
                                                                                     tGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTGTCCTtgatg                                                        >  3:94800/1‑147 (MQ=255)
                                                                                     tGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTGTCCTtgatg                                                        >  5:171778/1‑147 (MQ=255)
                                                                                       gtgtTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTGTCCTTGATGAAgg                                                    >  7:270438/1‑149 (MQ=255)
                                                                                         gtTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCAt                                                                                 >  8:71909/1‑118 (MQ=255)
                                                                                         gtTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCAt                                                                                 <  7:71909/118‑1 (MQ=255)
                                                                                          tttATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTGTCCTTGATGAAGGta                                                  >  4:276709/1‑148 (MQ=255)
                                                                                                    ttGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTGTCCTTGATGAAGGTATACAGCGGCtt                                       <  6:171778/149‑1 (MQ=255)
                                                                                                         cTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCcttc                                                                    >  4:175449/1‑115 (MQ=255)
                                                                                                         cTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCcttc                                                                    <  3:175449/115‑1 (MQ=255)
                                                                                                            gATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTGTCCTTGATGAAGGTATACAGCGGCTTGCcgtcg                                >  2:238783/1‑148 (MQ=255)
                                                                                                            gATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTGTCCTTGATGAAGGTATACAGCGGCTTGCcgtc                                 <  5:47811/147‑1 (MQ=255)
                                                                                                              ttACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTGTCCTTGATGAAGGTATACAGCGGCTTGCcgtcgt                               >  5:68306/1‑147 (MQ=255)
                                                                                                              ttACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTGTCCTTGATGAAGGTATACAGCGGCTTGCcgtcgt                               >  7:159261/1‑147 (MQ=255)
                                                                                                                    gTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTGTCCTTGATGAAGGTATACAGCGGCTTGCCGTCGTAAGCCCAc                       <  5:275025/149‑1 (MQ=255)
                                                                                                                              gatagccggggCCGCTCTGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTGTCCTTGATGAAGGTATACAGCGGCTTGCCGTCGTAAGCCCACTGCATCGAAc             <  3:191024/142‑1 (MQ=255)
                                                                                                                                   ccggggCCGCTCTGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTGTCCTTGATGAAGGTATACAGCGGCTTGCCGTCGTAAGCCCACTGCATCGAACcgtcg        >  3:197576/3‑149 (MQ=255)
                                                                                                                                       ggCCGCAGAGCGGCCCCGGGTTTCTCATGGCTTGGCAACGTGCCCGCCGTCTTTGATGGCGTCACCGGTGATGTGACCGGCCTTCTTGTGCTTGGTGAAGGTATACAGCGGCTTGGCGTCGTGAGGCCAGAGCATGGGACCgtcgtcca    >  8:265615/1‑147 (MQ=255)
                                                                                                                                         ccGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTGTCCTTGATGAAGGTATACAGCGGCTTGCCGTCGTAAGCCCACTGCATCGAACcgtc         <  5:303839/142‑1 (MQ=255)
                                                                                                                                         ccGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTGTCCTTGATGAAGGTATACAGCGGCTTGCCGTCGTAAGCCCACTGCATCGAACcgtc         >  6:303839/1‑142 (MQ=255)
                                                                                                                                         ccGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTGTCCTTGATGAAGGTATACAGCGGCTTGCCGTCGTAAGCCCACTGCATCGAACCGTCGTCAcgc  <  1:287022/149‑1 (MQ=255)
                                                                                                                                             agagCGGCCCCGGCTATCTCATGGCTGTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTGTCCTTGATGAAg                                                     >  1:3592/1‑95 (MQ=255)
                                                                                                                                             agagCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTGTCCTTGATGAAg                                                     <  2:3592/94‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                
TTTGCCTGGGAATACTGGTCAACGCGCAAGCTGGACCGGGAATCGTGCTTCAAGGCATTTCTGCACAACCACTGGGCGGGGATGCTGGTGTTTATCGGAGTTGTGCTGGATTACGCGTTGCGCTGAGTAGTTTGGGGCCGCAGAGCGGCCCCGGCTATCTCATGGCTTGGCAACGTGCCAGACGTCTTTCATGCCGTCACCGGTCATGTCACCGGCCTTCTTGTCCTTGATGAAGGTATACAGCGGCTTGCCGTCGTAAGCCCACTGCATCGAACCGTCGTCACGC  >  NC_002947/6064434‑6064719

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: