Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,386,712 A→T 6.7% C493S (TGT→AGT)  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,7120AT 6.7% 207.1 / 4.3 69C493S (TGT→AGT) aaeBp‑hydroxybenzoic acid efflux system component
Reads supporting (aligned to +/- strand):  ref base A (35/29);  new base T (2/3);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

GAGGTGGTTCTCTTTACGACGTGCCACATTGGTAGTCATCGCGGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATC  >  NC_000913/3386572‑3386858
                                                                                                                                            |                                                                                                                                                  
taGGTGGTTCTCTTTACGACGTGCCACATTGGTAGTCATCGCGGAAACAGCGGCAGAAACAAAATGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGACCGCGAGCACACAGCcgacg                                                                                                                                             <  6:145640/147‑1 (MQ=255)
  ggtggtTCTCTTTACGACGTGCCACATTGGTAGTCATCGCGGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGAt                                                                                                                                           <  2:142551/148‑1 (MQ=255)
  gggggtGCTCTTTACGACGTGCCACATTGGTAGGCATCGCGGAAACAGCGGCAGAAACAACCTGATAAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGAt                                                                                                                                           <  2:212126/145‑1 (MQ=255)
      gtTCTCTTTACGAGGTGCCACATTGGTAGTCATCTCGGAAACAGCGGGAGAAACAAAATGATTAAGCAGTACACGTCGGGTTCTGTCGCGCGATTTATCCCGCACCAGCAAAATACCGGTGAACGCTAGCACACTGCCGACGATTTgg                                                                                                                                       <  8:160555/148‑2 (MQ=255)
      gtTCTCTTTACGACGTGCCACATTGGTAGTCATCGCGGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACTGCCGACGATTTGc                                                                                                                                       >  1:45819/1‑148 (MQ=255)
      gtTCTCTTTACGACGTGCCACATTGGTAGTCATCGCGGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACTGCCGACGATTTGc                                                                                                                                       >  1:50474/1‑148 (MQ=255)
                  acgTGCCACATTGGTAGTCATCGCGGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTg                                                                                                                           >  5:116964/1‑148 (MQ=255)
                  acgTGCCACATTGGTAGTCATCGCGGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTg                                                                                                                           >  1:113044/1‑148 (MQ=255)
                       ccACATTGGTAGTCATCGCGGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCgag                                                                                                                      <  4:274366/148‑1 (MQ=255)
                          caTTGGTAGTCATCGCGGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGcacaca                                                                                                                                                    <  5:118460/115‑1 (MQ=255)
                          caTTGGTAGTCATCGCGGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGcacaca                                                                                                                                                    >  6:118460/1‑115 (MQ=255)
                            ttGGTAGTCATCGCGGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAact                                                                                                                 >  3:127538/1‑148 (MQ=255)
                              ggTAGACATCGCGGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGGCCTGGCGCGCGATTAATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAactg                                                                                                                <  6:210849/147‑1 (MQ=255)
                               gTAGTCATCGCGGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAactgac                                                                                                              <  4:291315/148‑1 (MQ=255)
                               gTAGTCATCGCGGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAactga                                                                                                               <  3:22755/147‑1 (MQ=255)
                                   tCATCGCGGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGaa                                                                                                          <  4:127538/148‑1 (MQ=255)
                                    cATCGCGGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGaaa                                                                                                         >  5:160232/1‑148 (MQ=255)
                                    cATCGCGGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGAGTTGCACTAATGCGCTGTCGAGAAACTGACTGaaa                                                                                                         >  6:213765/1‑148 (MQ=255)
                                       cgcgGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATg                                                                                                       >  4:144065/1‑147 (MQ=255)
                                         cgGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACCCGTCCGGTCCTCTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGCACGCGAGCACACAGCCGACGATTTGCGCTAATGCGCTGTCGAGAAACTGACTGAAATGGa                                                                                                     >  1:123286/1‑147 (MQ=255)
                                         cgGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGCCTGAAATGGa                                                                                                     >  2:142830/1‑147 (MQ=255)
                                         cgGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGa                                                                                                     >  3:291628/1‑147 (MQ=255)
                                              aCAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGACCGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCa                                                                                               <  5:197948/148‑1 (MQ=255)
                                               cAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCGTGTGGCGCTATTTATCCCGCACCTGCAAAATAAGGGGGAACGGGAGGTCACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCAt                                                                                              <  2:123286/148‑1 (MQ=255)
                                                 gCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCg                                                                                            <  2:113044/148‑1 (MQ=255)
                                                  cGGCAGAAACAAACTGATTAAGCAGTACACGGCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCgg                                                                                           <  3:114238/148‑1 (MQ=255)
                                                        aaacaaacTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATc                                                                                     >  4:188655/1‑148 (MQ=255)
                                                        aaacaaacTGATTAAGCAGTACACGTCCGGCCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTATCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATc                                                                                     >  3:93048/1‑148 (MQ=255)
                                                         aacaaacTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATc                                                                                             <  1:115402/139‑1 (MQ=255)
                                                          acaaacTGATAAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTCAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCACACGAGAAACAGCCTGAAAAGGAAAGTCAACg                                                                                            >  2:115402/1‑139 (MQ=255)
                                                           caaacTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGTTGAACGCGAGCACACAGCCGACGATTTGCCCTaa                                                                                                                                  <  5:295275/100‑1 (MQ=255)
                                                           caaacTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTaa                                                                                                                                  >  6:295275/1‑100 (MQ=255)
                                                            aaacTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGc                                                                                 <  4:291628/148‑1 (MQ=255)
                                                                       gCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATAtt                                                                      >  1:18016/1‑148 (MQ=255)
                                                                        cAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATAtt                                                                      >  8:39908/1‑147 (MQ=255)
                                                                            acacGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATTGAAAGTCATCGGGTTATCCAGCACGATAATATTTATgg                                                                 >  4:209437/1‑148 (MQ=255)
                                                                            acacGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAatat                                                                       >  8:3282/1‑142 (MQ=255)
                                                                            acacGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAatat                                                                       <  7:3282/142‑1 (MQ=255)
                                                                              acGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTg                                                               <  6:160232/148‑1 (MQ=255)
                                                                              acGGCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAACTGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTg                                                               <  4:220422/148‑1 (MQ=255)
                                                                                gTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCa                                                                                               <  8:171451/114‑1 (MQ=255)
                                                                                gTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCa                                                                                               >  7:171451/1‑114 (MQ=255)
                                                                                 tCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTg                                                            >  6:243071/1‑148 (MQ=255)
                                                                                 tCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTg                                                            >  5:27018/1‑148 (MQ=255)
                                                                                   cGGTCCTGTCTCTCGATTTATCCAGCACCAGCAAAATAACGGTGAAAGCGAGCACACAGCAGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTgg                                                           <  5:213765/147‑1 (MQ=255)
                                                                                    ggTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACTGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGcc                                                         <  2:50474/148‑1 (MQ=255)
                                                                                    ggTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACTGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGcc                                                         <  2:45819/148‑1 (MQ=255)
                                                                                    ggTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATTGTGCTGGcc                                                         <  3:209437/148‑1 (MQ=255)
                                                                                    ggTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGcc                                                         <  3:188655/148‑1 (MQ=255)
                                                                                    ggTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGcc                                                         <  5:243071/148‑1 (MQ=255)
                                                                                     gTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGc                                                          >  3:125155/1‑146 (MQ=255)
                                                                                     gTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCa                                                        <  3:144065/148‑1 (MQ=255)
                                                                                      tCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAg                                                       >  2:270146/1‑148 (MQ=255)
                                                                                       ccTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGt                                                      >  7:249637/1‑148 (MQ=255)
                                                                                         tGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTa                                                                                                                                   >  4:89097/1‑69 (MQ=255)
                                                                                         tGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTa                                                                                                                                   <  3:89097/69‑1 (MQ=255)
                                                                                         tGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGc                                                    >  1:304645/1‑148 (MQ=255)
                                                                                         tGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGc                                                    >  4:273061/1‑148 (MQ=255)
                                                                                            cgcgcgATTTATCCCGAACCAGCAAAATAACGGTGAAAGCGAGCACACAGCCGACGAATTGCCCTAATGCGCTGTCGAGAAAATGAATGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCAAGTGcacc                                                 >  5:266945/1‑148 (MQ=255)
                                                                                                             aCCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGcc                                <  1:272069/148‑1 (MQ=255)
                                                                                                                aGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGtt                             >  4:295516/1‑148 (MQ=255)
                                                                                                                 gCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGttt                            >  5:242862/1‑148 (MQ=255)
                                                                                                                 gCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGttt                            >  5:192699/1‑148 (MQ=255)
                                                                                                                   aaaaTAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCt                          <  6:192699/148‑1 (MQ=255)
                                                                                                                    aaaTAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGATTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTg                         >  3:262873/1‑148 (MQ=255)
                                                                                                                         aCGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACtt                    <  6:116964/148‑1 (MQ=255)
                                                                                                                              gAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCtat                >  7:178067/1‑147 (MQ=255)
                                                                                                                                   cgAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAg          >  4:296199/1‑148 (MQ=255)
                                                                                                                                     aGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGc                                                              >  3:149944/1‑94 (MQ=255)
                                                                                                                                     aGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGc                                                              <  4:149944/94‑1 (MQ=255)
                                                                                                                                     aGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCTAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAgg         >  3:208947/1‑147 (MQ=255)
                                                                                                                                        acacaGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGgaaga     <  8:178067/148‑1 (MQ=255)
                                                                                                                                           caGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCTGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATc  >  4:96813/1‑148 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
GAGGTGGTTCTCTTTACGACGTGCCACATTGGTAGTCATCGCGGAAACAGCGGCAGAAACAAACTGATTAAGCAGTACACGTCCGGTCCTGTCGCGCGATTTATCCCGCACCAGCAAAATAACGGTGAACGCGAGCACACAGCCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATC  >  NC_000913/3386572‑3386858

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