Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 531,520 C→A 18.7% intergenic (+1091/‑4634) USA300HOU_0494 → / → rplY hypothetical protein/ribosomal protein L25

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730531,5200CA18.7% 133.6 / 12.0 68intergenic (+1091/‑4634)USA300HOU_0494/rplYhypothetical protein/ribosomal protein L25
Reads supporting (aligned to +/- strand):  ref base C (31/23);  new base A (2/11);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

TCGATTAGAGCCAATTTCTGTTTCAAAACACCCACAAATCGATAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGGATAGAAGGGTTGAAAAGATGAGATATAAACGAAGCGAGAGAATTGTTTTTATGACGCAATATTTGATGAACC  >  CP000730/531380‑531655
                                                                                                                                             |                                                                                                                                       
tCGATTAGAGCCAATTTCTGTTTCAAAACACCCACAAATCGATAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGc                                                                                                                                         <  1:709470/141‑1 (MQ=255)
ccGATTAGAGCCAATTTCTGTTTCAAAACACCCACAAATCGACAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGc                                                                                                                                         <  1:281757/140‑1 (MQ=255)
  gATTAGAGCCAATTTCTGTTTCAAAACACCCACAAATCGATAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCaa                                                                                                                                       <  1:50882/141‑1 (MQ=255)
    ttAGAGTCAATTTCTGTTTCAAAACACCCACAAATCGATAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCaaaa                                                                                                                                     <  1:64885/141‑1 (MQ=255)
      agagCCAATTTCTGTTTCAAAACACCCACAAATCGATAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCaaaaaa                                                                                                                                   >  2:412090/1‑141 (MQ=255)
       gagCCAATTTCTGTTTCAAAACACCCTCAAATCGATAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCaaaaaaa                                                                                                                                  >  2:474136/1‑140 (MQ=255)
        agCCAATTTCTGTTTCAAAACACCCACAAATCGATAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCaaaaaaaa                                                                                                                                 >  1:433724/1‑139 (MQ=255)
        agCCAATTTCTGTTTCAAAACACCCACAAATCGATAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAata                                                                                                                                 >  2:411207/1‑141 (MQ=255)
        acccAATTTCTGTTTCAAAACACCCACAAAACGATAAAAAAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAAAAAAAAAAAAAata                                                                                                                                 <  2:121340/139‑1 (MQ=255)
         gCCAATTTCTGTTTCAAAACACCCACAAATCGATAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAatat                                                                                                                                >  2:671645/1‑141 (MQ=255)
         gCCAATTTCTGTTTCAAAACACCCACAAATCGATAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAatat                                                                                                                                >  2:470789/1‑141 (MQ=255)
         gCCAATTTCTGTTTCAAAACACCCACAAATCGATAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAatat                                                                                                                                >  2:599098/1‑141 (MQ=255)
          ccAATTTCTGTTTCAAAACACCCACAAATCGATAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATAtt                                                                                                                               >  2:766488/1‑141 (MQ=255)
          ccAATTTCTGTTTCAAAACACCCACAAATCGATAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATAtt                                                                                                                               <  1:360283/141‑1 (MQ=255)
            aaTTTCTGTTTAAAAACACCCACAAATCGATAAATTAAAACATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAACAAAAAAAATAttt                                                                                                                              >  1:550640/1‑141 (MQ=255)
                  tGTTTCAAAACAACCAAAAATCGATAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTCTATGGGCTAGCACGAAAAGAAAAACAAAAAAAAAATATTTATAATGc                                                                                                                       <  2:510320/141‑1 (MQ=255)
                  tGTTTAAAAAAAACAAAAAAACAATAAATTAAAAAATAATAAGTGGAAAAGTGGAGCAGATGGTGCGTTAATGAGAGGAAGCTGACCTACTGTGTATGGGCTAGCACGAAAAGAAAAAAAAAAAAAAAATATTTATAATGc                                                                                                                       <  2:281905/141‑1 (MQ=255)
                      tCAAAACACCCACAAATCGATAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGAACGAGAAGAAAGCCAAGCAAAAAAAATTAATAATGAAGtt                                                                                                                   >  2:553396/1‑141 (MQ=255)
                       cAAAACACCCACAAATCGATAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTa                                                                                                                  <  1:722267/141‑1 (MQ=255)
                        aaacaacCCACAAATTGATAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGTGGACTTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTaa                                                                                                                 <  2:14451/136‑1 (MQ=255)
                          aaaacaaaaaaaTCGATAAATTAAAAAATAATATGTTAAAAAGTGGTGCAGATGGTGCGTTAATCAGTGGAAGCGTACCTACTGTGTATGGGCTAGCACGAAAAGAAAAAAAAAAAAAAAATATTTATAATGCAGTTAACg                                                                                                               <  2:626334/132‑1 (MQ=255)
                            cacCCACAAATCGATAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGgt                                                                                                             >  1:523501/1‑141 (MQ=255)
                              cccACAAATCTATAAATTAAGAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAAAAAAAAAAAAAATATTTATAATGCAGTTAACGgttg                                                                                                           <  2:537278/141‑1 (MQ=255)
                               ccACAAATCGATAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGgttgt                                                                                                          <  1:711691/141‑1 (MQ=255)
                               ccACAAATCGATAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGgttgt                                                                                                          >  1:657658/1‑141 (MQ=255)
                                  caAATCGATAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTCGCACGAAAAGAAAGCCAAGCAACAAATATCTATAATGCAGTTAACGGTtgttgt                                                                                                       >  1:299306/1‑141 (MQ=255)
                                   aaaTCGATAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGATGCGTTGAGGAGTGGGAGCGGACCTACTGGGGATGGGCTGGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACTGTAGTTGTa                                                                                                      >  2:593160/1‑141 (MQ=255)
                                    aaTCGATAAATTAAAAAATAATATGTGGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGAAAAACAAAAAATATTTATAATGCAGTTAACGGTTGTTGTaa                                                                                                     <  2:454558/141‑1 (MQ=255)
                                      tCGATTGATTAAAAAATAAGAAGTTGAAAAGGGGAGCAGATGGTGCGTTAATGAGTGGAAGATGACCTACTGTGTATGGGCTAGCACGAAAAAAAAAAAAAAAAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATg                                                                                                   <  2:438124/141‑1 (MQ=255)
                                          tAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGTGGTTGGAATGAAgg                                                                                               >  1:630800/1‑140 (MQ=255)
                                          tAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGATGTTGTAATGAAgt                                                                                               >  2:516887/1‑141 (MQ=255)
                                          tAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGAAAAAAATCTTTCTAATGACGTTCAGGGTCATGGTGATGAAgg                                                                                               >  2:346889/1‑140 (MQ=255)
                                           aaaataaaaaataaTATGTTGAAAAGGGGAGCAGATGGGGAGTAAATGAGGGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAAAAAAAAAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAgtg                                                                                              <  2:50229/141‑1 (MQ=255)
                                              ttAAAATATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTAc                                                                                           <  1:369899/141‑1 (MQ=255)
                                                    agtaataTGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTCTACttagtt                                                                                     <  1:578687/139‑1 (MQ=255)
                                                     ataacaTGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTATAATGAAGTGTACttagtta                                                                                    >  1:536825/1‑141 (MQ=255)
                                                                 aaaGTGGTGCAGATGGTGCGTTAATGAGGGGAAGCGTACCGACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGGAt                                                                        <  1:694251/141‑1 (MQ=255)
                                                                      ggtgCAGGTGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGAGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATAAAGTGTACTTAGTTAGACTATTAGGAGAGAAg                                                                   >  2:603331/1‑141 (MQ=255)
                                                                       gtgCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTGACGGTTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGGATAGAAgg                                                                  >  2:235125/1‑141 (MQ=255)
                                                                        tgCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATAAAGTGTACTTAGTTAGACTATTAGGACAGAAggg                                                                 >  1:471496/1‑141 (MQ=255)
                                                                         gCAGATGGGGGCTTAAAGAATTGAAGCAGACATACGTTGTATGGGCTAGCACTAAAAGAAAAAAAAAAAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGGATAGAAGGGt                                                                <  2:663929/141‑1 (MQ=255)
                                                                               ggtgggtgaagtggTGGACGCGGACCTAATGCTTATTGGCTAGCACGAAAAGAAAGCCACGCAAAAAATAGTTACAAGTCATTTAACGGTGGTTGTAATGGAGTGTACTTAGTTAGACTATTAGGATAGAAGGGTGGaaaa                                                          <  2:259824/128‑1 (MQ=255)
                                                                               ggTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGGATAGAAGGGTTGaaaa                                                          <  1:445367/141‑1 (MQ=255)
                                                                               ggTGAGTTAATGAGTGGAAGCGGAACTACTGTGTATGGGCTAGCACGAAAAGAAAGAAAAAAAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGGATAGAAGGGTTGaaaa                                                          <  2:99982/141‑1 (MQ=255)
                                                                                   cGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGGATAGAAGGGTTGAAAAGATg                                                      >  1:294148/1‑141 (MQ=255)
                                                                                       agtGAGGGGAAGCGGAACTACTGTGGATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGGAGTTAACGGGTGTTGTAATGAAGTGTATTTAGTTCGACTGTGAGGATAGAAGGGTTGAAAAGATGAGat                                                  <  1:382623/139‑1 (MQ=255)
                                                                                       aaTGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGGATAGAAGGGTTGAAAAGATGAGat                                                  <  1:750617/141‑1 (MQ=255)
                                                                                       aaTGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGGATAGAAGGGTTGAAAAGATGAGat                                                  <  1:673413/141‑1 (MQ=255)
                                                                                           actGGAGGCGGACCTCATGTGTAAGTGATAACACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTTTACTTAGTTAGACTATTAGGATAGAAGTGTTGAAAAGATGAGATATaa                                              <  2:37643/139‑1 (MQ=255)
                                                                                            gTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCCGATAACGGTTGTTGTAATGAAGTGTACTTTGTTAGACTATTAGGATAGAAGGTTTGAAAAGAGTAGATATaaa                                             >  1:503380/1‑141 (MQ=255)
                                                                                             tGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTAGTTAGTTAGACTATTAGGATAGATGGGTTGAAAAGGTGAGATATaaag                                            >  1:401759/1‑140 (MQ=255)
                                                                                             tGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGGATAGAAGGGTTGAAAAGATGAGATATAAAc                                            <  2:379025/141‑1 (MQ=255)
                                                                                                 aGCGGGCCTACTGTGTATGGGCTAGCACGAAAAGAATGCCACACAAAAAATTTTTATCATGGAGTTAACGGTTGTTGTAATGAAATGTACTTGCTAAGACTATATGGATAGAAGGGTAGAAAAGCTGAGATATAAACGaaa                                        >  2:308401/1‑140 (MQ=255)
                                                                                                       ccTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGGATAGAAGGGTTGAAAAGATGAGATATAAACGAAGCgagag                                  <  2:369965/141‑1 (MQ=255)
                                                                                                       ccTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGGATAGAAGGGTTGAAAAGATGAGATATAAACGAAGCgagag                                  <  2:638753/141‑1 (MQ=255)
                                                                                                        cTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGGATAGAAGGGTTGAAAAGATGAGATATAAACGAAGCgagaga                                 <  1:194032/141‑1 (MQ=255)
                                                                                                              tgtATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGGATAGAAGGGTTGAAAAGATGAGATATAAACGAAGTGAGAGAATTGtt                           <  1:410147/141‑1 (MQ=255)
                                                                                                                    ggCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGGATAGAAGGGTTGAAAAGATGAGATATAAACGAAGCGAGAGAATTGTTTTTATg                     >  2:708802/1‑141 (MQ=255)
                                                                                                                    ggCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGATGTTAACGGTTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGGATAGAAGAGTTGAAAAGATGAGATATAAACGAAGCGAGAGAATTGTTTTTATg                     >  2:720149/1‑141 (MQ=255)
                                                                                                                     gCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGAGTACTTAGTTAGACTATTAGGATAGAAGGGTTGAAAAGATGAGATATAAACGAAGCGAGAGAATTGTTTTTATGa                    >  2:775264/1‑141 (MQ=255)
                                                                                                                          cACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGATAGAAGGGTTGAAAAGATGAGATATAAACGAAGCGAGAGAATTGTTTTTATGACGCAat              >  2:49738/1‑141 (MQ=255)
                                                                                                                            cGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGGATAGAAGGGTTGAAAAGATGAGATATAAACGAAGCGAGAGAATTGTTTTTATGACGCAata             >  1:505835/1‑141 (MQ=255)
                                                                                                                            cGAAAAAAAAAAAAAGAAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGGATAGAAGGGTTGAAAAGATGAGATATAAACGAAGCGAGAGAATTGTTTTTATGACGCAata             <  2:359106/141‑1 (MQ=255)
                                                                                                                            agAAAAAAAAACAAAAAAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGGATAGAAGGGTTGAAAAGATGAGATATAAACGAAGCGAGAGAATTGTTTTTATGACGCAata             <  2:3153/140‑1 (MQ=255)
                                                                                                                                aagaaagCCAAGCAAAAAAGATTTATAAGGCAGTTAACGGTGGGTGTAATGAAGTGTACTTAGTTGGCATAATAGGATCGAAGGGATGCACAGATGAGAAAGAAACGAAGCAAGAGAAATGTTTTTCTGACGTAACAttta         >  2:267740/1‑140 (MQ=255)
                                                                                                                                aagaaagCCAACCAAAAAAGATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGGATAGAAGGGTTGAAAAGATGAGATATAAACGAAGCGAGAGAATTGTTTTTATGACGCAATATTtg         <  1:177312/141‑1 (MQ=255)
                                                                                                                                 agaaagCCAAGGAAAAGATATTTACAACGCAGGTAACGGGTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGGATAGAAGGGTTGAAAAGATGAGATATAAACGAAGCGAGAGAATTGTTTTTATGACGCAATATTtga        <  2:25901/141‑1 (MQ=255)
                                                                                                                                 agaaagCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGGATAGAAGGGTTGAAAAGATGAGATATAAACGAAGCGAGAGAATTGTTTTTATGACGCAATATTtga        <  2:722451/141‑1 (MQ=255)
                                                                                                                                    aagCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGGATAGAAGGGTTGAAAAGATGAGATATAAACGAAGCGAGAGAATTGTTTTTATGACGCAATATTtgatga     >  2:705679/1‑141 (MQ=255)
                                                                                                                                       ccAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTACTTAGTTAGAATATTAGGATAGAAGGGTTGAAAAGATGAGATATAGACGAAGCGAGAGAATTGTTTTTATGACGCAATATTTGATGAAcc  >  2:242522/1‑141 (MQ=255)
                                                                                                                                             |                                                                                                                                       
TCGATTAGAGCCAATTTCTGTTTCAAAACACCCACAAATCGATAAATTAAAAAATAATATGTTGAAAAGTGGTGCAGATGGTGCGTTAATGAGTGGAAGCGGACCTACTGTGTATGGGCTAGCACGAAAAGAAAGCCAAGCAAAAAATATTTATAATGCAGTTAACGGTTGTTGTAATGAAGTGTACTTAGTTAGACTATTAGGATAGAAGGGTTGAAAAGATGAGATATAAACGAAGCGAGAGAATTGTTTTTATGACGCAATATTTGATGAACC  >  CP000730/531380‑531655

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: