Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,279,708 A→T 7.3% L18Q (CTG→CAG)  bcr ← bicyclomycin/cysteine/sulfonamide efflux transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,279,7080AT 7.3% 324.6 / 12.3 95L18Q (CTG→CAG) bcrbicyclomycin/cysteine/sulfonamide efflux transporter
Reads supporting (aligned to +/- strand):  ref base A (32/56);  new base T (2/5);  total (34/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

CCGTAGATTAACTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAATCAGCATC  >  NC_000913/2279564‑2279853
                                                                                                                                                |                                                                                                                                                 
ccGTAGATTAACTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCcagcag                                                                                                                                               >  3:454866/1‑149 (MQ=255)
  gTAGATTAACTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGc                                                                                                                                             >  2:575642/1‑149 (MQ=255)
     gATTAACTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCaa                                                                                                                                          <  2:168122/149‑1 (MQ=255)
     gATTAACTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCaa                                                                                                                                          >  2:252079/1‑149 (MQ=255)
      aTTAACTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAg                                                                                                                                         <  3:791981/149‑1 (MQ=255)
      aTTAACTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAg                                                                                                                                         <  2:562359/149‑1 (MQ=255)
       ttAACTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCTGCTGGCCatgt                                                                                                                                        <  1:407508/149‑4 (MQ=255)
       ttAACTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAgg                                                                                                                                        >  2:778651/1‑149 (MQ=255)
       ttAACTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTAATCTGCGTACTGCCCGCCGGTACGCCAAACGGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCTGGTGGCCatgt                                                                                                                                        <  3:323314/149‑4 (MQ=255)
        tAACTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCTGCTGGCCatgtt                                                                                                                                       <  1:841487/149‑5 (MQ=255)
        tAACTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGa                                                                                                                                       <  3:47942/149‑1 (MQ=255)
           cTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCTGCTGGCCatgttgat                                                                                                                                    <  1:728157/149‑8 (MQ=255)
            tGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATaaa                                                                                                                                   <  1:778651/149‑1 (MQ=255)
              ccccAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATaaaaa                                                                                                                                 <  1:62552/149‑1 (MQ=255)
              ccccAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATaaaaa                                                                                                                                 <  4:237176/149‑1 (MQ=255)
              ccccAACGCAAAGCCCAGAATATAAGTACTGAGGCTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATaaaaa                                                                                                                                 <  4:201421/149‑1 (MQ=255)
                ccAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACa                                                                                                                               >  1:623602/1‑149 (MQ=255)
                 cAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACaa                                                                                                                              <  1:479350/149‑1 (MQ=255)
                 cAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACaa                                                                                                                              >  1:68521/1‑149 (MQ=255)
                  aaCGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAAt                                                                                                                             >  1:493013/1‑149 (MQ=255)
                  aaCGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAAt                                                                                                                             >  3:839839/1‑149 (MQ=255)
                   aCGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATa                                                                                                                            <  4:449843/149‑1 (MQ=255)
                   aCGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATa                                                                                                                            <  4:657749/149‑1 (MQ=255)
                   aCGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATa                                                                                                                            <  3:369894/149‑1 (MQ=255)
                           cccAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGTAGGCCAAGGATAAAAACAATAGCAAacga                                                                                                                    <  1:117826/149‑1 (MQ=255)
                           cccAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAacga                                                                                                                    <  1:369431/149‑1 (MQ=255)
                            ccAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAacgac                                                                                                                   <  2:184164/149‑1 (MQ=255)
                             cAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCTGCTGGCCatgt                                                                                                                                        <  4:660604/127‑4 (MQ=255)
                             cAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCTGCTGGCCatgt                                                                                                                                        >  3:660604/1‑124 (MQ=255)
                             cAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAacgacg                                                                                                                  >  3:777789/1‑149 (MQ=255)
                             cAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAacgacg                                                                                                                  >  2:618285/1‑149 (MQ=255)
                             cAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAacgacg                                                                                                                  >  4:608465/1‑149 (MQ=255)
                                aaTATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAAt                                                                                                               >  4:144907/1‑149 (MQ=255)
                                   ataAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCt                                                                                                            <  1:512168/149‑1 (MQ=255)
                                    taAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTg                                                                                                           >  1:292594/1‑149 (MQ=255)
                                        tACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATTGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCgg                                                                                                       <  3:452432/149‑1 (MQ=255)
                                           tGAGGGTCATCTGCGGACTGCCCGCCGGTAGGCCAAACGGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGagtg                                                                                                    <  3:116544/149‑1 (MQ=255)
                                                gTCATCTGCGTACTACCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTcac                                                                                               >  3:67487/1‑149 (MQ=255)
                                                   aTCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAAc                                                                                            >  2:696870/1‑149 (MQ=255)
                                                      tGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACggg                                                                                         >  4:760322/1‑149 (MQ=255)
                                                       gCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGGAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGc                                                                                        <  3:39453/149‑1 (MQ=255)
                                                        cGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCgg                                                                                                       >  2:422642/1‑133 (MQ=255)
                                                        cGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCgg                                                                                                       <  1:422642/133‑1 (MQ=255)
                                                         gTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTc                                                                                      <  4:13217/149‑1 (MQ=255)
                                                          tACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTcc                                                                                     <  2:753918/149‑1 (MQ=255)
                                                          tACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGACAAGGATAAAAACAATAGCAAACGACGCATGCTGTCGGGTGGTCACAACGGGCTcc                                                                                     <  3:810499/149‑1 (MQ=255)
                                                             tGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGa                                                                                  <  4:189505/149‑1 (MQ=255)
                                                              gCCCGCCGATACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGaa                                                                                 >  1:118567/1‑149 (MQ=255)
                                                               cccGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGaaa                                                                                <  2:426202/149‑1 (MQ=255)
                                                               cccGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGaaa                                                                                <  1:429224/149‑1 (MQ=255)
                                                               cccGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGaaa                                                                                <  2:266337/149‑1 (MQ=255)
                                                               cccGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGaaa                                                                                <  2:493013/149‑1 (MQ=255)
                                                                     ggTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCAtt                                                                          <  2:717562/149‑1 (MQ=255)
                                                                          gCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGaaaa                                                                     <  4:533996/149‑1 (MQ=255)
                                                                          gCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGaaaa                                                                     <  2:86127/149‑1 (MQ=255)
                                                                          gCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGaaaa                                                                     >  2:403042/1‑149 (MQ=255)
                                                                          gCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGaaaa                                                                     <  1:27390/149‑1 (MQ=255)
                                                                           ccAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCt                                                                                                            >  3:826719/1‑109 (MQ=255)
                                                                           ccAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCt                                                                                                            <  4:826719/109‑1 (MQ=255)
                                                                             aaaCTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGtt                                                                  <  1:750599/149‑1 (MQ=255)
                                                                               aCTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAg                                                                >  1:73111/1‑149 (MQ=255)
                                                                                cTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAga                                                               >  1:617337/1‑149 (MQ=255)
                                                                                  gcgcTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAgacg                                                             >  1:752669/1‑149 (MQ=255)
                                                                                                  tAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATttctt                                             <  4:67487/149‑1 (MQ=255)
                                                                                                   aGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATttcttc                                            <  3:372001/149‑1 (MQ=255)
                                                                                                       cgGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCa                                        <  4:77342/149‑1 (MQ=255)
                                                                                                       cgGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCa                                        <  4:26305/149‑1 (MQ=255)
                                                                                                       cgGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCa                                        <  2:73111/149‑1 (MQ=255)
                                                                                                        gggCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAg                                       >  3:71043/1‑149 (MQ=255)
                                                                                                          gCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGtt                                     <  4:777789/149‑1 (MQ=255)
                                                                                                            aGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTaa                                   >  1:835559/1‑149 (MQ=255)
                                                                                                                aCATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGa                               <  2:848656/149‑1 (MQ=255)
                                                                                                                aCATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGa                               <  1:98526/149‑1 (MQ=255)
                                                                                                                 cATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGAc                              <  1:437064/149‑1 (MQ=255)
                                                                                                                   tatCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGa                            <  2:623602/149‑1 (MQ=255)
                                                                                                                     tCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGata                          >  1:71023/1‑149 (MQ=255)
                                                                                                                       aaTCGACAGCGGCATCAACATGGCCAGCAGGGCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATAtt                        >  2:884552/1‑149 (MQ=255)
                                                                                                                       aaTCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATAtt                        <  1:403042/149‑1 (MQ=255)
                                                                                                                        aTCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTc                       >  4:875006/1‑149 (MQ=255)
                                                                                                                           gACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCAcc                    >  4:357529/1‑149 (MQ=255)
                                                                                                                            aCAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCg                   <  4:454866/149‑1 (MQ=255)
                                                                                                                              aGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCggg                 <  2:426920/149‑1 (MQ=255)
                                                                                                                               gCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCgggg                <  1:575642/149‑1 (MQ=255)
                                                                                                                               gCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAAAGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATGGCAAAGTTAGACGACGCAGGCAATTTCTTCTTCAGTTAACGTACGTTATTCACCgggg                >  2:13213/1‑149 (MQ=255)
                                                                                                                                   cATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGc                                                         <  4:375057/104‑1 (MQ=255)
                                                                                                                                   cATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGc                                                         >  3:375057/1‑104 (MQ=255)
                                                                                                                                     tCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTaaa          <  1:645912/149‑1 (MQ=255)
                                                                                                                                     tCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTaaa          <  2:835559/149‑1 (MQ=255)
                                                                                                                                     tCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTaaa          <  3:86925/149‑1 (MQ=255)
                                                                                                                                        aCATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAATCa       >  2:105329/1‑149 (MQ=255)
                                                                                                                                        aCATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAATCa       <  4:189740/149‑1 (MQ=255)
                                                                                                                                         cttGGCCTGCTGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAATCAg      >  1:608898/3‑149 (MQ=255)
                                                                                                                                          aTGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAATCAGc     <  1:696870/149‑1 (MQ=255)
                                                                                                                                           tGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAATCAGCa    >  4:328615/1‑149 (MQ=255)
                                                                                                                                           tGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAATCAGCa    <  1:327051/149‑1 (MQ=255)
                                                                                                                                             gCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAATCAGCATc  <  3:156511/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CCGTAGATTAACTGCCCCAACGCAAAGCCCAGAATATAAGTACTGAGGGTCATCTGCGTACTGCCCGCCGGTACGCCAAACTGCGCTGAAATTACCGGTAGCGCGGGCAGATACATATCAATCGACAGCGGCATCAACATGGCCAGCAGGCCAAGGATAAAAACAATAGCAAACGACGAATGCTGTCGGGTGGTCACAACGGGCTCCTGAAAGTCATTGAAAAGTTAGACGACGCTGGCAATTTCTTCTTCAGTTAACGGACGATATTCACCGGGGGCTAAATCAGCATC  >  NC_000913/2279564‑2279853

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