Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,657,805 G→T 8.9% intergenic (‑194/+15) pabB ← / ← PP_2330 aminodeoxychorismate synthase/para‑aminobenzoate synthase multi‑enzyme complex/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,657,8050GT 8.9% 200.8 / 6.5 64intergenic (‑194/+15)pabB/PP_2330aminodeoxychorismate synthase/para‑aminobenzoate synthase multi‑enzyme complex/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/28);  new base T (2/4);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCAGGAAAAGCGCCTAGCGCTGGGTACGACATTCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTGCGCCTCGATGGCATGGGTGCGCTGAGTCAGGCCGCAGT  >  NC_002947/2657658‑2657953
                                                                                                                                                   |                                                                                                                                                    
tgCAGGAAAAGCGCCTAGCGCTGGGTACGACATTCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAg                                                                                                                                                      <  3:201268/148‑1 (MQ=255)
  cAGGAAAAGCGCCTAGCGCTGGGTACGACATTCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAgg                                                                                                                                                     >  3:270423/1‑147 (MQ=255)
         aGCGCCTAGCGCTGGGTACGACATTCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCCTTGCAGcccatt                                                                                                                                            <  8:384351/149‑4 (MQ=255)
             ccTAGCGCTGGGTACGACATTCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTAtt                                                                                                                                        >  8:312570/1‑149 (MQ=255)
               tAGCGCTGGGTACGACATTCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAgg                                                                                                                                                     <  6:196878/134‑1 (MQ=255)
                aGCTCTGGGTACTACTTTCAACATTGTGTTTGGACCATACTCAGCACTAGTCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTTCTGCGAATGGGCTGCAATGCAGCCCCTGTATTTc                                                                                                                                      <  2:407875/148‑1 (MQ=255)
                   gcTGGGTACGACATTCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCg                                                                                                                                  <  4:270423/149‑1 (MQ=255)
                        gTACGACATTCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTgg                                                                                                                             <  7:312570/149‑1 (MQ=255)
                         tACGACATTCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTggg                                                                                                                            >  5:35846/1‑149 (MQ=255)
                            gACATTCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGca                                                                                                                          >  3:153536/1‑148 (MQ=255)
                                ttCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGc                                                                                                                     >  1:366288/1‑149 (MQ=255)
                                ttCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGc                                                                                                                     >  5:95310/1‑149 (MQ=255)
                                   aaCATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTgg                                                                                                                             <  3:275883/138‑1 (MQ=255)
                                   aaCATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTgg                                                                                                                             >  4:275883/1‑138 (MQ=255)
                                   aaCATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGccc                                                                                                                   <  2:284993/148‑1 (MQ=255)
                                   aaCATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGccc                                                                                                                   <  2:366288/148‑1 (MQ=255)
                                      aTTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGaa                                                                                                                >  8:315850/1‑148 (MQ=255)
                                          tgtTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGccc                                                                                                                                               <  8:362890/113‑1 (MQ=255)
                                          tgtTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGccc                                                                                                                                               >  7:362890/1‑113 (MQ=255)
                                                gCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCtgtg                                                                                                     <  6:95310/149‑1 (MQ=255)
                                                gCCCATACGCAGCACAAGGCAGGTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCtgt                                                                                                      <  4:77684/148‑1 (MQ=255)
                                                    aTACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCCTTGCAGcccatt                                                                                                                                            >  2:59369/1‑103 (MQ=255)
                                                    aTACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCCTTGCAGcccatt                                                                                                                                            <  1:59369/106‑4 (MQ=255)
                                                     tACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCtacc                                                                                                <  6:35846/149‑1 (MQ=255)
                                                      aCGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCtaccc                                                                                               >  2:154313/1‑149 (MQ=255)
                                                      aCGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCtaccc                                                                                               >  1:167965/1‑149 (MQ=255)
                                                              caAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCgc                                                                                       <  4:153536/149‑1 (MQ=255)
                                                                 ggCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCt                                                                                                                                             >  3:95632/1‑92 (MQ=255)
                                                                 ggCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCt                                                                                                                                             <  4:95632/92‑1 (MQ=255)
                                                                   cAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTc                                                                                  <  7:315850/149‑1 (MQ=255)
                                                                     gTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGtt                                                                                                                               >  6:221315/1‑102 (MQ=255)
                                                                     gTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGtt                                                                                                                               <  5:221315/102‑1 (MQ=255)
                                                                          tGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCt                                                                           >  2:304279/1‑149 (MQ=255)
                                                                              ggTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCgg                                                                       <  2:167965/149‑1 (MQ=255)
                                                                               gTCGTGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGAc                                                                                                           <  5:357944/112‑1 (MQ=255)
                                                                               gTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGAc                                                                                                           >  6:357944/1‑112 (MQ=255)
                                                                               gTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGt                                                                      >  6:110047/1‑149 (MQ=255)
                                                                                  gAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAg                                                                                                                                    <  1:169151/84‑1 (MQ=255)
                                                                                  gAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAg                                                                                                                                    >  2:169151/1‑84 (MQ=255)
                                                                                      gAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGt                                                                             <  4:124332/135‑1 (MQ=255)
                                                                                      gAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGt                                                                             >  3:124332/1‑135 (MQ=255)
                                                                                         gcgcTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGACCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACcctcc                                                            >  6:333133/1‑149 (MQ=255)
                                                                                           gcTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCtt                                                          >  3:191896/1‑149 (MQ=255)
                                                                                               gAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGc                                                       <  1:154313/148‑1 (MQ=255)
                                                                                                  tCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGAcc                                                                                                          <  4:406523/94‑1 (MQ=255)
                                                                                                  tCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGAcc                                                                                                          >  3:406523/1‑94 (MQ=255)
                                                                                                  tCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGt                                                   >  5:352174/1‑149 (MQ=255)
                                                                                                    aaCTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTg                                                  <  6:352174/148‑1 (MQ=255)
                                                                                                         cAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCGTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTc                                              <  5:333133/147‑1 (MQ=255)
                                                                                                            tGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGTTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGt                                                   >  7:358629/1‑139 (MQ=255)
                                                                                                            tGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGTTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGt                                                   <  8:358629/139‑1 (MQ=255)
                                                                                                               ttGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGTCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACCCCATCCTTGGCCAGTGCTTCGAtgtgtg                                      >  2:144270/1‑147 (MQ=255)
                                                                                                               ttGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTg                                                         >  8:374281/1‑130 (MQ=255)
                                                                                                               ttGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTg                                                         <  7:374281/130‑1 (MQ=255)
                                                                                                               ttGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTgcg                                      >  3:195128/1‑149 (MQ=255)
                                                                                                                 gCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTGCGcc                                    >  1:156827/1‑149 (MQ=255)
                                                                                                                     gagCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTGCGCCTCGa                                >  4:114548/1‑149 (MQ=255)
                                                                                                                      agCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTGCGCCTCGAt                               <  2:126708/149‑1 (MQ=255)
                                                                                                                          aCCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTGCGCCTCGATGGCa                           >  5:320902/1‑149 (MQ=255)
                                                                                                                                         ggggCTGCCTTGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTGCGCCTCGATGGCATgg                        >  2:239738/2‑137 (MQ=255)
                                                                                                                                         ggggCTGCCTTGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTGCGCCTCGATGGCATgg                        <  1:239738/136‑1 (MQ=255)
                                                                                                                                            gCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCtaccc                                                                                               >  6:331408/1‑63 (MQ=255)
                                                                                                                                            gCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCtaccc                                                                                               <  5:331408/63‑1 (MQ=255)
                                                                                                                                            gCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTGCGCCTCGATGGCATGGGTGCGCTGAGTCAgg         <  4:191896/149‑1 (MQ=255)
                                                                                                                                                 aaGGGAGCCCCTGTTTTTCAGCGGTTGGGCAGATGCCCGAACAGACCCTGTGCTACCCGTAGCCGCTGCTCGGGGTCTTCGGTCACACCCTCCTTGGGCAGTGCTTCGATGTGCGCCTCGATGGCATGGGTGCGCTGAGTCAGGCCGCa    <  4:203609/149‑1 (MQ=255)
                                                                                                                                                 aaGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTGCGCCTCGATGGCATGGGTGCGCTGAGTCAGGCCGCa    <  2:156827/149‑1 (MQ=255)
                                                                                                                                                 aaGGCAGCCCCTGTATTTAAGCGGTTGGCCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGACTTCGGTCACACCCTCCTTGGCCAGAGCTTCGATGTGCGCCTCGATGGCATGGGTGCGCTGAGTCAGGCCGCa    >  3:203609/1‑149 (MQ=255)
                                                                                                                                                   ggCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTACGATGTGCGCCTCGATGGCATGGGTCCGCAGAGTCAGGCCGCAGt  >  8:73184/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGCAGGAAAAGCGCCTAGCGCTGGGTACGACATTCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTGCGCCTCGATGGCATGGGTGCGCTGAGTCAGGCCGCAGT  >  NC_002947/2657658‑2657953

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑