Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,657,804 A→T 7.8% intergenic (‑193/+16) 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,8040AT 7.8% 195.4 / 5.3 64intergenic (‑193/+16)pabB/PP_2330aminodeoxychorismate synthase/para‑aminobenzoate synthase multi‑enzyme complex/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base T (2/3);  total (32/32)
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.44e-01

TGCAGGAAAAGCGCCTAGCGCTGGGTACGACATTCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTGCGCCTCGATGGCATGGGTGCGCTGAGTCAGGCCGCA  >  NC_002947/2657658‑2657951
                                                                                                                                                  |                                                                                                                                                   
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)
                                                                                                                                                  |                                                                                                                                                   
TGCAGGAAAAGCGCCTAGCGCTGGGTACGACATTCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTGCGCCTCGATGGCATGGGTGCGCTGAGTCAGGCCGCA  >  NC_002947/2657658‑2657951

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: ‑