Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,387,553 T→A 8.2% R212R (CGA→CGT 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,387,5530TA 8.2% 320.7 / 16.1 97R212R (CGA→CGTaaeBp‑hydroxybenzoic acid efflux system component
Reads supporting (aligned to +/- strand):  ref base T (44/45);  new base A (5/3);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.27e-01

GCGTATTCTGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAAGAG  >  NC_000913/3387405‑3387693
                                                                                                                                                    |                                                                                                                                            
gcgTATTCTGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGt                                                                                                                                              >  2:577076/1‑149 (MQ=255)
gcgTATTCTGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACATCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGt                                                                                                                                              <  4:352600/149‑1 (MQ=255)
 cgTATTCTGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTc                                                                                                                                             <  1:757004/149‑1 (MQ=255)
  gTATTCTGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCg                                                                                                                                            >  1:512078/1‑149 (MQ=255)
   tATTCTGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGc                                                                                                                                           <  3:292124/149‑1 (MQ=255)
    aTTCTGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCa                                                                                                                                          <  1:514794/149‑1 (MQ=255)
       cTGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCa                                                                                                                                       >  2:449303/1‑149 (MQ=255)
        tGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCTCACCAg                                                                                                                                      >  4:580464/1‑149 (MQ=255)
        tGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAg                                                                                                                                      <  1:574623/149‑1 (MQ=255)
        tGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAg                                                                                                                                      >  1:565197/1‑149 (MQ=255)
          aaTAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGt                                                                                                                                    >  2:764634/1‑149 (MQ=255)
                attagTTTGGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGTGCCGTGGTGCGTCGCACCAGGTCGcccc                                                                                                                              >  4:307116/4‑149 (MQ=255)
                ataagTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGcccc                                                                                                                              >  4:269635/1‑149 (MQ=255)
                ataagTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGcccc                                                                                                                              >  3:703194/1‑149 (MQ=255)
                ataagTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGcccc                                                                                                                              >  1:464597/1‑149 (MQ=255)
                  aagTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAg                                                                                                                            <  3:186014/149‑1 (MQ=255)
                    gTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGc                                                                                                                          <  3:809365/149‑1 (MQ=255)
                     tttCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCt                                                                                                                         <  1:464406/149‑1 (MQ=255)
                       tCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCttt                                                                                                                       <  4:300733/149‑1 (MQ=255)
                         gCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTAt                                                                                                                     <  1:661263/149‑1 (MQ=255)
                         gCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTAt                                                                                                                     >  2:661263/1‑149 (MQ=255)
                         gCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTAt                                                                                                                     >  3:432865/1‑149 (MQ=255)
                         gCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGCCGCCCCAGGCTTTAt                                                                                                                     <  4:845398/149‑1 (MQ=255)
                          cAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATc                                                                                                                    <  1:14921/149‑1 (MQ=255)
                               ttGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACaa                                                                                                               >  3:647386/1‑149 (MQ=255)
                                 gggTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACt                                                                                                             >  3:627700/1‑149 (MQ=255)
                                   gTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTc                                                                                                           >  2:340628/1‑149 (MQ=255)
                                     aaTCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCAc                                                                                                         <  2:340315/149‑1 (MQ=255)
                                     aaTCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCAc                                                                                                         <  1:788596/149‑1 (MQ=255)
                                       tcagcgtcagcgATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCa                                                                                                          <  4:713524/146‑1 (MQ=255)
                                       tcagcgtcagcgATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCa                                                                                                          >  3:713524/1‑146 (MQ=255)
                                             tcagcgATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCa                                                                                                 <  3:248015/149‑1 (MQ=255)
                                             tcagcgATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCa                                                                                                 >  3:214861/1‑149 (MQ=255)
                                                gcgATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGc                                                                                              <  1:764634/149‑1 (MQ=255)
                                                 cgATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCt                                                                                             <  1:833993/149‑1 (MQ=255)
                                                     aGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGAt                                                                                         <  2:368546/149‑1 (MQ=255)
                                                        gTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACa                                                                                      <  1:449303/149‑1 (MQ=255)
                                                        gTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCATCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACa                                                                                      >  1:624522/1‑149 (MQ=255)
                                                          aTTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACaga                                                                                    <  4:279714/149‑1 (MQ=255)
                                                          aTTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACaga                                                                                    <  2:464597/149‑1 (MQ=255)
                                                           ttGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACagag                                                                                   <  1:836176/149‑1 (MQ=255)
                                                             gATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGtt                                                                                 >  2:295213/1‑149 (MQ=255)
                                                              aTCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTg                                                                                >  4:225294/1‑149 (MQ=255)
                                                                  cTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCAtt                                                                            >  3:167484/1‑149 (MQ=255)
                                                                   ttttAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCAttt                                                                           >  1:792211/1‑148 (MQ=255)
                                                                     ttAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAAt                                                                         <  1:600287/149‑1 (MQ=255)
                                                                     ttAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAAt                                                                         <  1:388510/149‑1 (MQ=255)
                                                                        aaCGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCAt                                                                                                <  2:547846/123‑1 (MQ=255)
                                                                        aaCGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCAt                                                                                                >  1:547846/1‑123 (MQ=255)
                                                                           gTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGatat                                                                   >  2:97178/1‑149 (MQ=255)
                                                                              gATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTgc                                                                <  2:512078/149‑1 (MQ=255)
                                                                              gATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTgc                                                                >  2:446207/1‑149 (MQ=255)
                                                                                   gCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGAc                                                           <  1:856690/149‑1 (MQ=255)
                                                                                     ccGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCa                                                         >  3:643393/1‑149 (MQ=255)
                                                                                          cccAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCaaa                                                    <  1:446207/149‑1 (MQ=255)
                                                                                           ccAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACagatt                                                                                  <  2:149664/118‑3 (MQ=255)
                                                                                           ccAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGt                                                                                  >  1:149664/1‑118 (MQ=255)
                                                                                           ccAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATAtt                                                                  >  2:30490/1‑134 (MQ=255)
                                                                                           ccAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATAtt                                                                  <  1:30490/134‑1 (MQ=255)
                                                                                            cAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACt                                                  >  3:866327/1‑149 (MQ=255)
                                                                                                gggAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTcc                                              >  3:170532/1‑149 (MQ=255)
                                                                                                gggAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTcc                                              <  2:792211/149‑1 (MQ=255)
                                                                                                 ggAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCa                                             <  4:646206/149‑1 (MQ=255)
                                                                                                   aaGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGc                                           >  2:494205/1‑149 (MQ=255)
                                                                                                       ttCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCtctc                                       >  4:234798/1‑149 (MQ=255)
                                                                                                       ttCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCtctc                                       >  4:325796/1‑149 (MQ=255)
                                                                                                       ttCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCtctc                                       >  3:524254/1‑149 (MQ=255)
                                                                                                       ttCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCtctc                                       >  4:321469/1‑149 (MQ=255)
                                                                                                          cATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGAt                                    >  1:278958/1‑149 (MQ=255)
                                                                                                          cATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGAt                                    >  2:630086/1‑149 (MQ=255)
                                                                                                          cATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCCCGACCAGCAAACTTTCCAGCTCTCGAt                                    >  4:579454/1‑149 (MQ=255)
                                                                                                           atatTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATc                                   <  1:519205/149‑1 (MQ=255)
                                                                                                              ttCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCAc                                >  2:265946/1‑149 (MQ=255)
                                                                                                                   gTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCActtctt                           <  1:186800/149‑1 (MQ=255)
                                                                                                                     tgctgcGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGt                         <  2:125210/149‑1 (MQ=255)
                                                                                                                           gcATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTgatc                   <  3:321469/149‑1 (MQ=255)
                                                                                                                             aTGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTgatcga                 <  4:524254/149‑1 (MQ=255)
                                                                                                                              tGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTgatcgat                >  2:88020/1‑149 (MQ=255)
                                                                                                                                 cTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATcgc             <  4:385902/149‑1 (MQ=255)
                                                                                                                                 cTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATcgc             >  4:806493/1‑149 (MQ=255)
                                                                                                                                    gTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGga          <  3:69742/149‑1 (MQ=255)
                                                                                                                                    gTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTAACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCAACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGga          <  3:672097/149‑1 (MQ=255)
                                                                                                                                       gcgcCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGaa       >  3:166097/1‑149 (MQ=255)
                                                                                                                                         gcCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGATCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGaaaa     <  1:577076/149‑1 (MQ=255)
                                                                                                                                         gcCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGaaaa     <  3:829471/149‑1 (MQ=255)
                                                                                                                                         gcCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGaaaa     <  3:334169/149‑1 (MQ=255)
                                                                                                                                          ccGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATTCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAag    >  4:156485/1‑149 (MQ=255)
                                                                                                                                          ccGTGATGCGTCGCACCAGGTCGCCCCAGGCTTAATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCTGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAag    >  1:254715/1‑149 (MQ=255)
                                                                                                                                          acctggtgCGACGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAag    >  3:436354/4‑149 (MQ=255)
                                                                                                                                          acctggtgCGACGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAag    >  4:54316/4‑149 (MQ=255)
                                                                                                                                           cGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAaga   <  1:340628/149‑1 (MQ=255)
                                                                                                                                           cGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAaga   <  4:627700/149‑1 (MQ=255)
                                                                                                                                            gtggtgCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAagag  >  2:856434/1‑149 (MQ=255)
                                                                                                                                            ctggtgCGACGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCa                                                         <  2:736165/93‑1 (MQ=255)
                                                                                                                                            ctggtgCGACGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCa                                                         >  1:736165/2‑94 (MQ=255)
                                                                                                                                            ctggtgCGACGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAagag  >  1:125076/2‑149 (MQ=255)
                                                                                                                                            ctggtgCGACGCACCACGGCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAagag  <  4:87100/148‑1 (MQ=255)
                                                                                                                                                tgCGACGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAaga   <  2:658654/144‑1 (MQ=255)
                                                                                                                                                tgCGACGCACCAGGTCGCCCCAGGCTTTATAGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAaga   >  1:658654/1‑144 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
GCGTATTCTGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGCCCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAAGAG  >  NC_000913/3387405‑3387693

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: