Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,386,794 G→C 6.3% T465T (ACC→ACG)  aaeB ← p‑hydroxybenzoic acid efflux system component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,386,7940GC 6.3% 257.9 / 4.3 77T465T (ACC→ACG) aaeBp‑hydroxybenzoic acid efflux system component
Reads supporting (aligned to +/- strand):  ref base G (30/41);  new base C (2/3);  total (32/45)
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

TACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAA  >  NC_000913/3386647‑3386940
                                                                                                                                                   |                                                                                                                                                  
tACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATgg                                                                                                                                                   >  3:10280/1‑149 (MQ=255)
   acGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGc                                                                                                                                                >  3:67203/1‑149 (MQ=255)
    cGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCt                                                                                                                                               <  4:67203/149‑1 (MQ=255)
      tCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTgg                                                                                                                                             <  7:209927/149‑1 (MQ=255)
        cGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGc                                                                                                                                            >  1:366803/1‑148 (MQ=255)
         ggTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCa                                                                                                                                          <  2:65434/149‑1 (MQ=255)
                    gcgATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATc                                                                                                                               <  4:11851/149‑1 (MQ=255)
                      gATTTATCCCGCCCCAGCAAAATAACGGGGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCCCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGGTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGa                                                                                                                             <  6:183424/149‑1 (MQ=255)
                      gATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGa                                                                                                                             <  7:306382/149‑1 (MQ=255)
                       aTTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAg                                                                                                                            >  5:166739/1‑149 (MQ=255)
                           aTCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACCATAAATATTATCGTGCTggat                                                                                                                                           <  3:28686/130‑3 (MQ=255)
                           aTCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACCATAAATATTATCGTGCTggat                                                                                                                                           >  4:28686/1‑128 (MQ=255)
                            tCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAg                                                                                                                       >  6:293692/1‑149 (MQ=255)
                            tCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAg                                                                                                                       <  8:139137/149‑1 (MQ=255)
                                gCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGc                                                                                                                   >  5:8489/1‑149 (MQ=255)
                                gCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGc                                                                                                                   >  5:3485/1‑149 (MQ=255)
                                  aCCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCg                                                                                                                 <  5:145802/149‑1 (MQ=255)
                                     aGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGttt                                                                                                              <  4:10280/149‑1 (MQ=255)
                                     aGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGttt                                                                                                              >  1:220582/1‑149 (MQ=255)
                                     aGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGttt                                                                                                              >  7:333663/1‑149 (MQ=255)
                                          aaTAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTa                                                                                                         <  8:327849/149‑1 (MQ=255)
                                              aCGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACtt                                                                                                      <  2:274397/148‑1 (MQ=255)
                                              aCGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTc                                                                                                     <  3:103273/149‑1 (MQ=255)
                                                ggTGAACGCGAGCACACAGCCGACGATTTACCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCta                                                                                                   >  3:308547/1‑149 (MQ=255)
                                                 gTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCta                                                                                                   >  7:19397/1‑148 (MQ=255)
                                                    aaCGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATAcc                                                                                               >  4:210251/1‑149 (MQ=255)
                                                      cgcgAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCGGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGa                                                                                             <  3:51403/149‑1 (MQ=255)
                                                          aGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGgaa                                                                                         <  8:19397/149‑1 (MQ=255)
                                                            ctcacaGCCGACGAATTGCCCTAATGCGCAGTCTCGAAACTCACTGAAATGGAAAGTCATCCGGTTATACAGCACGATAATAATTATGGCGCAGGCCTCTGCCCCCATCGAGCTCAGGCGCCGGTTCTGTACTTCTATACCGAGgaa                                                                                         <  5:105280/145‑1 (MQ=255)
                                                             acacaGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGGTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGgaagaa                                                                                      <  2:366803/149‑1 (MQ=255)
                                                              cacaGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAAt                                                                                     <  3:154489/149‑1 (MQ=255)
                                                              cacaGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTAATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCGGTTTCTGTACTTCTATACCGAGGAAGAAt                                                                                     <  2:231249/149‑1 (MQ=255)
                                                               acaGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATc                                                                                    >  5:47231/1‑149 (MQ=255)
                                                                  gCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCa                                                                                 <  3:210251/149‑1 (MQ=255)
                                                                   ccGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAATGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAg                                                                                <  6:353621/149‑1 (MQ=255)
                                                                   ccGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAg                                                                                >  2:245244/1‑149 (MQ=255)
                                                                    cgacgaTTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAg                                                                                >  5:333573/1‑148 (MQ=255)
                                                                     gacgaTTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCa                                                                              >  4:202494/1‑149 (MQ=255)
                                                                     gacgaTTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCa                                                                              <  2:220582/149‑1 (MQ=255)
                                                                     gacgaTTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCa                                                                              <  1:245244/149‑1 (MQ=255)
                                                                              cccTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGc                                                                     >  4:85402/1‑149 (MQ=255)
                                                                                  aaTGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAAt                                                                 >  3:33445/1‑149 (MQ=255)
                                                                                   ctGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAAt                                                                 <  1:271825/147‑1 (MQ=255)
                                                                                   aTGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTaa                                                                  <  6:3485/147‑1 (MQ=255)
                                                                                   aTGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTaa                                                                  <  6:8489/147‑1 (MQ=255)
                                                                                      cgcTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGcac                                                             <  6:166739/149‑1 (MQ=255)
                                                                                       gcTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGCGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGca                                                              >  7:130684/1‑147 (MQ=255)
                                                                                             gagaAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGcc                                                                                                                          >  4:178734/1‑81 (MQ=255)
                                                                                             gagaAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGcc                                                                                                                          <  3:178734/81‑1 (MQ=255)
                                                                                             gagaAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGgaag                                                                                        <  1:350837/115‑1 (MQ=255)
                                                                                             gagaAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGgaag                                                                                        >  2:350837/1‑115 (MQ=255)
                                                                                                aaaCTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACa                                                   <  7:135621/149‑1 (MQ=255)
                                                                                                    tgactgaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCCAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGc                                                >  5:12429/1‑148 (MQ=255)
                                                                                                        tgaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCCGTCGCCGTTTGTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTATTGCACAGCAGCAACATGCTCTGt                                           >  1:160632/1‑149 (MQ=255)
                                                                                                        tgaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGt                                           >  5:220717/1‑149 (MQ=255)
                                                                                                        tgaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGt                                           >  4:12512/1‑149 (MQ=255)
                                                                                                        tgaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTGATGGTGCTGGCCAGTGCCCCCCTCGAGCCCAGGCGCCGTTGCGGTACTTCTATACCGAGGATGAATCCCAGCACTGCCCGGCTAATGCACCGCAGCAACATGCTCTGt                                           >  6:175673/1‑149 (MQ=255)
                                                                                                          acaTGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATTGTGCTGGCCAGTGCCCCCATCGTGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGtt                                          <  2:275802/146‑1 (MQ=255)
                                                                                                          aaaTGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGtt                                          <  1:214749/148‑1 (MQ=255)
                                                                                                           aaTGGAAAGTCATCGGGTTATCCAGCACGATAATATTAATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTgg                                        <  8:354563/149‑1 (MQ=255)
                                                                                                                aaaGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTAtt                                   <  5:293692/149‑1 (MQ=255)
                                                                                                                  aGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAg                                 <  8:333663/149‑1 (MQ=255)
                                                                                                                  aGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAg                                 <  6:47231/149‑1 (MQ=255)
                                                                                                                         gggTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATaa                           <  3:85402/148‑1 (MQ=255)
                                                                                                                         gcgTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGGCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAAt                          <  6:353995/147‑1 (MQ=255)
                                                                                                                              aTCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCAAGTCGCCATTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCAcc                      >  1:300825/1‑148 (MQ=255)
                                                                                                                               gccAGCACCATAAATATTATCGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCa                     >  5:288633/2‑148 (MQ=255)
                                                                                                                                  aGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCaaa                   >  8:8286/1‑147 (MQ=255)
                                                                                                                                    cACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAg                >  6:114134/1‑148 (MQ=255)
                                                                                                                                    cACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAg                >  4:269319/1‑148 (MQ=255)
                                                                                                                                       gATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTaga            >  1:272236/1‑149 (MQ=255)
                                                                                                                                        ataataTTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTagag           <  2:160394/149‑1 (MQ=255)
                                                                                                                                        ataataTTTATCGTGCTGTCAAGTGACCCCATCGACCCCAGTCGCAGTTTCTGTACATCTATACCGAGGAAGAATCCCAGCACTTCAATGCTAAAGCACAGCAGCAACATGCTCTGTAGGGAATTAGGGATAATCACCAAAAAGTagag           <  1:93859/149‑1 (MQ=255)
                                                                                                                                           atatTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGagca         <  4:33445/148‑1 (MQ=255)
                                                                                                                                            atattATCGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGagcagc       <  6:288633/146‑1 (MQ=255)
                                                                                                                                              tttATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGcc      <  4:308547/148‑1 (MQ=255)
                                                                                                                                              tttATGGGGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGccc     <  2:160632/149‑1 (MQ=255)
                                                                                                                                                  tGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTaa  <  8:140055/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAA  >  NC_000913/3386647‑3386940

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