Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,387,553 T→A 6.8% 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 6.8% 296.4 / 6.3 85R212R (CGA→CGT) aaeBp‑hydroxybenzoic acid efflux system component
Reads supporting (aligned to +/- strand):  ref base T (43/36);  new base A (3/3);  total (46/39)
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

CGTATTCTGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAAGAG  >  NC_000913/3387406‑3387693
                                                                                                                                                   |                                                                                                                                             
cgTATTCTGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTc                                                                                                                                              <  3:136705/149‑1 (MQ=255)
     tCTGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCAc                                                                                                                                          >  2:612992/1‑148 (MQ=255)
       tGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCa                                                                                                                                        <  2:284917/148‑1 (MQ=255)
       tGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAg                                                                                                                                       >  3:703888/1‑149 (MQ=255)
        gAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGACGCACCacg                                                                                                                                      <  3:888158/149‑3 (MQ=255)
             agataagTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGACCTGGTGCGACGCACCACGgcgc                                                                                                                                  <  1:632569/148‑1 (MQ=255)
              gataagTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGc‑cc                                                                                                                               >  4:405073/1‑149 (MQ=255)
               ataagTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCa                                                                                                                                        <  4:340666/140‑1 (MQ=255)
               ataagTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCa                                                                                                                                        >  3:340666/1‑140 (MQ=255)
                  agTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAgg                                                                                                                           >  4:777918/1‑149 (MQ=255)
                   gTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGc                                                                                                                          >  4:324000/1‑149 (MQ=255)
                         cAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTAt                                                                                                                     >  1:446537/1‑148 (MQ=255)
                            gATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGAc                                                                                                                 >  3:672833/1‑149 (MQ=255)
                              ttGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACa                                                                                                                <  2:382886/148‑1 (MQ=255)
                                    aaTCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCa                                                                                                          >  2:342248/1‑148 (MQ=255)
                                        agcgtcagcgATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATc                                                                                                     >  3:771003/1‑149 (MQ=255)
                                          cgtcagcgATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCg                                                                                                    >  1:186579/1‑148 (MQ=255)
                                            tcagcgATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATc                                                                                                                    >  3:268612/1‑130 (MQ=255)
                                            tcagcgATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATc                                                                                                                    <  4:268612/130‑1 (MQ=255)
                                            tcagcgATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCa                                                                                                 <  3:777918/149‑1 (MQ=255)
                                             cagcgATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCa                                                                                                 >  1:740457/1‑148 (MQ=255)
                                                   tAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTg                                                                                           >  2:678540/1‑148 (MQ=255)
                                                   tAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGa                                                                                          <  4:179562/149‑1 (MQ=255)
                                                     gCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATa                                                                                        >  4:375953/1‑149 (MQ=255)
                                                      cGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGAt                                                                                         <  1:859668/147‑1 (MQ=255)
                                                      cGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGAt                                                                                         >  2:859668/1‑147 (MQ=255)
                                                             aTCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTg                                                                                >  4:355454/1‑149 (MQ=255)
                                                              tCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGc‑cc                                                                                                                               >  3:681213/1‑101 (MQ=255)
                                                              tCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGc‑cc                                                                                                                               <  4:681213/101‑1 (MQ=255)
                                                              tCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑TCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGc                                                                               >  4:758558/1‑149 (MQ=255)
                                                              tCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTg                                                                                >  2:250961/1‑148 (MQ=255)
                                                              tCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTg                                                                                >  1:493333/1‑148 (MQ=255)
                                                              tCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTg                                                                                >  2:710576/1‑148 (MQ=255)
                                                               cGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGc                                                                               <  2:493333/148‑1 (MQ=255)
                                                                gCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCa                                                                              <  2:793263/148‑1 (MQ=255)
                                                                   tttAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTaa                                                                          <  3:405073/149‑1 (MQ=255)
                                                                       aaCGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGa                                                                      <  3:375953/149‑1 (MQ=255)
                                                                        aCGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGat                                                                     <  3:544228/149‑1 (MQ=255)
                                                                        aCGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGCTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACa                                                                                      >  4:489373/1‑132 (MQ=255)
                                                                        aCGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGCTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACa                                                                                      <  3:489373/132‑1 (MQ=255)
                                                                           tCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGatat                                                                   >  2:92129/1‑148 (MQ=255)
                                                                             gATTGGCCCGCGCCGAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTgc                                                                >  3:63882/1‑149 (MQ=255)
                                                                              aTTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTgc                                                                >  2:267342/1‑148 (MQ=255)
                                                                               ttGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTgcgc                                                              >  4:743809/1‑149 (MQ=255)
                                                                                    ccGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGatat                                                                   <  1:806282/139‑1 (MQ=255)
                                                                                    ccGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGatat                                                                   >  2:806282/1‑139 (MQ=255)
                                                                                     cgcgCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCa                                                         >  1:758947/1‑148 (MQ=255)
                                                                                      gcgcCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGc                                                       <  3:599551/149‑1 (MQ=255)
                                                                                         cccAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCaaa                                                    <  3:303234/149‑1 (MQ=255)
                                                                                          ccAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAAc                                                   >  4:721126/1‑149 (MQ=255)
                                                                                           cAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTgcgc                                                              <  3:681001/137‑1 (MQ=255)
                                                                                           cAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTgcgc                                                              >  4:681001/1‑137 (MQ=255)
                                                                                           cAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTCGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAAc                                                   >  2:155445/1‑148 (MQ=255)
                                                                                              cGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACAAGCAAACTTTc                                               <  4:771003/149‑1 (MQ=255)
                                                                                                ggAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCa                                             >  3:184120/1‑149 (MQ=255)
                                                                                                 gAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCa                                             <  2:273431/148‑1 (MQ=255)
                                                                                                   aGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGct                                          >  4:331221/1‑149 (MQ=255)
                                                                                                    gATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGct                                          <  1:92129/148‑1 (MQ=255)
                                                                                                      ttCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCtctc                                       >  3:20173/1‑149 (MQ=255)
                                                                                                           tatTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATcc                                  >  3:58745/1‑149 (MQ=255)
                                                                                                           tatTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATcc                                  >  4:549772/1‑149 (MQ=255)
                                                                                                              tCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTcc                                              >  2:80299/1‑134 (MQ=255)
                                                                                                              tCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTcc                                              <  1:80299/134‑1 (MQ=255)
                                                                                                                 ggTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCActtct                            <  3:324000/149‑1 (MQ=255)
                                                                                                                 ggTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCActtc                             <  1:710576/148‑1 (MQ=255)
                                                                                                                    tgctgcGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGt                         <  4:703888/149‑1 (MQ=255)
                                                                                                                     gctgcGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGt                         <  1:612992/148‑1 (MQ=255)
                                                                                                                         cgcATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTgat                    <  3:331221/149‑1 (MQ=255)
                                                                                                                               ccTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGttt                       <  4:645686/140‑1 (MQ=255)
                                                                                                                               ccTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGttt                       >  3:645686/1‑140 (MQ=255)
                                                                                                                                 ttGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATcgc             >  1:455964/1‑148 (MQ=255)
                                                                                                                                     ggcgacctggtgCGACGCACCACGGCGCTACAAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGga          >  2:423396/8‑148 (MQ=38)
                                                                                                                                      gcgcCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGaa       >  3:274917/1‑149 (MQ=255)
                                                                                                                                        gcCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGaaaa     <  3:721126/149‑1 (MQ=255)
                                                                                                                                        gcCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTCGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGaaa      <  1:155445/148‑1 (MQ=255)
                                                                                                                                         ccGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGaaaa     >  2:148610/1‑148 (MQ=255)
                                                                                                                                         acctggtgCGACGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCaaa                                                    <  3:633553/98‑1 (MQ=255)
                                                                                                                                         acctggtgCGACGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCaaa                                                    >  4:633553/4‑101 (MQ=255)
                                                                                                                                          cctggtgCGACGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAaga   >  4:431760/3‑149 (MQ=255)
                                                                                                                                          cGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAaga   <  4:672833/149‑1 (MQ=255)
                                                                                                                                          cGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAaga   <  4:274917/149‑1 (MQ=255)
                                                                                                                                          cGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAag    <  1:771842/148‑1 (MQ=255)
                                                                                                                                           gtggtgCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAagag  <  3:355454/149‑1 (MQ=255)
                                                                                                                                                gCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAaga   <  3:504640/143‑1 (MQ=255)
                                                                                                                                                gCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAaga   >  4:504640/1‑143 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
CGTATTCTGAATAAGATAAGTTTCGCAGGATTGGGTAATCAGCGTCAGCGATAGCGTATTGATCGCTTTTAAACGTCGATTGGCCCGCGCCCAGCGGGAAGATTCCATATTCAGGTTGCTGCGCATGCCTTGTAGCGCCGTGGTGCGTCGCACCAGGTCGC‑CCCAGGCTTTATCGACAACTTCACCATCGCCATGCTTGATACAGAGTTGCATTAATTGATATTGCGCGACCAGCAAACTTTCCAGCTCTCGATCCACTTCTTGTTTGATCGATCGCGGAGAAAAGAG  >  NC_000913/3387406‑3387693

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