Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,657,804 A→T 12.6% 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,8040AT12.6% 163.3 / 11.7 63intergenic (‑193/+16)pabB/PP_2330aminodeoxychorismate synthase/para‑aminobenzoate synthase multi‑enzyme complex/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/27);  new base T (4/4);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGAAAAGCGCCTAGCGCTGGGTACGACATTCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTGCGCCTCGATGGCATGGGTGCGCTGAGTCAGGCCG  >  NC_002947/2657660‑2657949
                                                                                                                                                |                                                                                                                                                 
cAGGAAAAGCGCCTAGCGCTGGGTACGACATTCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCa                                                                                                                                               <  5:163479/149‑1 (MQ=255)
   gAAAAGCGCCTAGCGCTGGGTACGACATTCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGcc                                                                                                                                            >  7:220823/1‑149 (MQ=255)
    aaaaGCGCCTAGCGCTGGGTACGACATTCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGccc                                                                                                                                           <  1:138949/149‑1 (MQ=255)
          gcCTAGCGCTGGGTACGACATTCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGACCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTAt                                                                                                                                     <  3:52334/149‑1 (MQ=255)
           ccTAGCGCTGGGTACGACATTCAACATTGTGTTTGGTCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCCTTGCAGcccattcg                                                                                                                                      <  8:86362/147‑6 (MQ=255)
           ccTAGCGCTGGGTACGACATTCAACATTGTGTTTGGTCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCCTTGCAGcccattcg                                                                                                                                      >  7:86362/1‑142 (MQ=255)
            cTAGCGCTGGGTACGACATTCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTAtt                                                                                                                                    >  4:19868/1‑148 (MQ=255)
                        aCGACATTCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTggg                                                                                                                        >  4:267208/1‑148 (MQ=255)
                         cGACATTCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGca                                                                                                                      <  8:220823/149‑1 (MQ=255)
                          gACATTCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGca                                                                                                                      >  6:201619/1‑148 (MQ=255)
                              ttCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGc                                                                                                                 >  8:180209/1‑149 (MQ=255)
                              ttCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGc                                                                                                                 <  5:201619/149‑1 (MQ=255)
                              ttCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAAGGGGCTGCCTTGCAGcccatt                                                                                                                                        >  5:190844/1‑123 (MQ=255)
                              ttCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAAGGGGCTGCCTTGCAGcccatt                                                                                                                                        <  6:190844/126‑4 (MQ=255)
                               tCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGcc                                                                                                                >  6:290433/1‑149 (MQ=255)
                                 aaCATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGccc                                                                                                               <  8:10375/148‑1 (MQ=255)
                                    aTTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGaa                                                                                                            >  7:113344/1‑148 (MQ=255)
                                    aTTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAAc                                                                                                           <  3:19868/149‑1 (MQ=255)
                                         gtTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCg                                                                                                              <  2:196691/141‑1 (MQ=255)
                                         gtTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCg                                                                                                              >  1:196691/1‑141 (MQ=255)
                                          tttGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGAccc                                                                                                     >  8:132567/1‑149 (MQ=255)
                                          tttGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTCAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGAccc                                                                                                     >  6:289057/1‑149 (MQ=255)
                                               cccATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGc                                                                                                >  4:184680/1‑149 (MQ=255)
                                               cccATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGc                                                                                                <  8:113344/149‑1 (MQ=255)
                                                ccATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGc                                                                                                >  2:22070/1‑148 (MQ=255)
                                                 cATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCta                                                                                              <  7:132567/149‑1 (MQ=255)
                                                      gcagcaCAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACAt                                                                                                                   <  1:227868/123‑1 (MQ=255)
                                                      gcagcaCAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACAt                                                                                                                   >  2:227868/1‑123 (MQ=255)
                                                         gcaCAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCtacccgtacc                                                                                      <  7:180209/149‑1 (MQ=255)
                                                          cacaAGGCAGTACCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCtacccgtac                                                                                       >  7:228420/1‑147 (MQ=255)
                                                                  cgtgccTGCAGGGCGAGCGACGCGCTTGAATCAACTGCAGTGCTTGCTAGCGCACCCTTGTGCTGCGAAGGGGCGGCAAGGCCGGCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGTCCCTGTGCTACCCGTACCCGCTGcgcg                                                                             <  6:309571/145‑1 (MQ=255)
                                                                     tCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGc                                                                           <  5:290433/148‑1 (MQ=255)
                                                                        tGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCt                                                                       >  4:281837/1‑149 (MQ=255)
                                                                         gCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCtac                                                                                             <  1:299545/126‑1 (MQ=255)
                                                                         gCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCtac                                                                                             >  2:299545/1‑126 (MQ=255)
                                                                              tCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTc                                                                 <  3:267208/149‑1 (MQ=255)
                                                                                          cTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTg                                                     <  3:184680/149‑1 (MQ=255)
                                                                                          cTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTg                                                     <  7:39809/149‑1 (MQ=255)
                                                                                                tCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGGAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGt                                               >  5:215721/1‑149 (MQ=255)
                                                                                                tCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGt                                               >  8:10569/1‑149 (MQ=255)
                                                                                                tCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGt                                               >  2:278470/1‑149 (MQ=255)
                                                                                                tCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGt                                               <  3:281837/149‑1 (MQ=255)
                                                                                                        aGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGAt                                       <  5:289057/149‑1 (MQ=255)
                                                                                                         gTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGAtg                                      <  1:278470/149‑1 (MQ=255)
                                                                                                             ttGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTgcg                                  >  3:250426/1‑149 (MQ=255)
                                                                                                              tGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTgcgc                                 >  1:140872/1‑149 (MQ=255)
                                                                                                                   gagCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCacac                                                             >  3:74839/1‑116 (MQ=255)
                                                                                                                   gagCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCacac                                                             <  4:74839/116‑1 (MQ=255)
                                                                                                                   gagCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGc                                             >  7:249974/1‑132 (MQ=255)
                                                                                                                   gagCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGc                                             <  8:249974/132‑1 (MQ=255)
                                                                                                                                   cGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCtt                                           >  6:300786/1‑118 (MQ=255)
                                                                                                                                   cGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCtt                                           <  5:300786/118‑1 (MQ=255)
                                                                                                                                   cGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTGCGCCTCGATGGCATGGGTgcgc              <  2:140872/147‑1 (MQ=255)
                                                                                                                                    gAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTGCGCCTCGATGGCATGGGTGCGCt             >  6:306009/1‑147 (MQ=255)
                                                                                                                                     caggggCTGCCTTGCAGCCCCTGTATTTCAGCGGTTGGGCATATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCtt                                           <  3:289227/113‑1 (MQ=255)
                                                                                                                                     caggggCTGCCTTGCAGCCCCTGTATTTCAGCGGTTGGGCATATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCtt                                           >  4:289227/4‑116 (MQ=255)
                                                                                                                                      aggggCTGCCTTGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTGCGCCTCGATGGCATGGGTGCGCTGAg          <  1:307754/146‑1 (MQ=255)
                                                                                                                                      aggggCTGCCTTGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTGCGCCTCGATGGCATGGGTGCGCTGAg          >  2:307754/3‑148 (MQ=255)
                                                                                                                                      aTGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTGCGCCTCGATGGCATGGGTGCGCTGAGt         >  2:168886/1‑149 (MQ=255)
                                                                                                                                         gtcTGCAAGTCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGAGCCTGTGCTACCCGTACGCGCTGCTCGGCGTCTTCGGTCACTCCCTCCTTGGCCAGTGCTTCGTTGTGCGCCTCGATGGCATGGGTGCGCTGAGTCAg      <  8:228420/147‑1 (MQ=255)
                                                                                                                                          gCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTgc                                   >  7:205486/1‑119 (MQ=255)
                                                                                                                                          gCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTgc                                   <  8:205486/119‑1 (MQ=255)
                                                                                                                                              cAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTGCGCCTCGATGGCATGGGTGCGCTGAGTCAGGCCg  >  3:167143/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CAGGAAAAGCGCCTAGCGCTGGGTACGACATTCAACATTGTGTTTGGCCCATACGCAGCACAAGGCAGTTCCTGCAGGTCGAGCGAAGCGCTTGAATCAACTGCAGTGCTTGCAAGAGCAACCTTGTGCTGCGAATGGGCTGCAAGGCAGCCCCTGTATTTCAGCGGTTGGGCACATGCCCGAACAGACCCTGTGCTACCCGTACCCGCTGCTCGGCGTCTTCGGTCACACCCTCCTTGGCCAGTGCTTCGATGTGCGCCTCGATGGCATGGGTGCGCTGAGTCAGGCCG  >  NC_002947/2657660‑2657949

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