Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,504,957 G→C 11.7% intergenic (+261/+143) USA300HOU_2365 → / ← gltT hypothetical protein/DAACS family dicarboxylate/amino acid:sodium (Na+) or proton (H+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,504,9570GC11.7% 177.9 / 14.4 78intergenic (+261/+143)USA300HOU_2365/gltThypothetical protein/DAACS family dicarboxylate/amino acid:sodium (Na+) or proton (H+) symporter
Reads supporting (aligned to +/- strand):  ref base G (25/43);  new base C (4/6);  total (29/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

TTGAATAGCTCCACTTGCCGTTCGCCTGCGCTTTGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTAT  >  CP000730/2504818‑2505096
                                                                                                                                           |                                                                                                                                           
ttGAATAGCTCCACTTGCCGTTCGCCTGCGCTTTACACATGCATAGAAGCCCCTAACGACATGAGGTCACTGCGCTTCTGTTCGCCTGCGCTTTCGGTATTGCTTAAAACACACTAACAACGTGAGGTCACTATGCTTCgg                                                                                                                                            >  2:123365/1‑141 (MQ=255)
  gAATAGCTCCACTTGCCGTTCGCCTGCGCTTTGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTCCGCCTGCGCGTTAGCGCATGCATAAAAGCCCCTAACACCCTGAGGTCACTACGCTTCGGtt                                                                                                                                          <  1:122448/141‑1 (MQ=255)
  gAATAGCTCCACTTGCCGTTCGCCTGCGCTTTGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTGCGCCCCGGTTCGCCTGCGCTTTAACGCATGCATAAAAGCCCCTAACAACCGGGGGTCGCTGCGCTTCggat                                                                                                                                          >  2:98448/1‑139 (MQ=18)
    aTAGCTCCACTTGCCGTTCGCCTGCGCTTTGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCCGTTCGCCGGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCg                                                                                                                                        <  2:116252/141‑1 (MQ=255)
    aTAGCTCCACTTGCCGTTCGCCTGCGCTTTGCGCATGCATAAAAGCCCCCAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCg                                                                                                                                        <  2:285122/141‑1 (MQ=255)
     tAGCTCCACTTGCCGTTCGCCTGCGCTTTGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGGGCTTTAGCGCCTGCATAAAAGCCCCTAACAACCTGAGGTCACTCCGCTTCGTTTCGc                                                                                                                                       >  1:135653/1‑141 (MQ=255)
         cccACTTGCCGTCCGCCTGCGCTTTGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTgc                                                                                                                                   <  2:15605/140‑1 (MQ=255)
            aCTTGCCGTTCGCCTGCGCTTTGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCAGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGACTGTGCt                                                                                                                                >  2:213235/1‑141 (MQ=255)
             cTTGCCGTTCGCCTGCGCTTTGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCta                                                                                                                               <  2:35576/141‑2 (MQ=255)
                 ccGTTCGCCTGCGCTTTGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCCTGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCCGGGCTTTatc                                                                                                                           >  2:210402/1‑139 (MQ=255)
                 ccGTTCGCCTGCGCTTTGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAgc                                                                                                                           >  1:238157/1‑141 (MQ=255)
                    ttCGCCTGCGCTTTGCGCATGCATAAAAGCCCCTAACAACCTTAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCgca                                                                                                                        <  2:153646/141‑1 (MQ=255)
                       gCCTGCGCTTTGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCgcatgc                                                                                                                     >  1:119988/1‑141 (MQ=255)
                        ccTGCGCTTTGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATCAAAGCCCCTAACAACCTGAGGTCACTACGATTCCGTTCGCCTGCGCGTTAgcaaaagca                                                                                                                    >  2:223591/1‑134 (MQ=14)
                         cTGCGCTTTGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCgcatgcat                                                                                                                   <  1:89008/141‑1 (MQ=255)
                                  gcggATGCATAAAGGCGCCTAACAACCTGAGGTCACTCCGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTccc                                                                                                          <  2:123257/141‑1 (MQ=11)
                                  gcgcATGCATAAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTAGGGGTCGGTTCGCCTAAGATTTAGCGGATGCGTAAAAGTgcc                                                                                                          >  2:80425/1‑141 (MQ=14)
                                       tgcatAAAAGTCCCTAACCACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTaaca                                                                                                     <  1:275045/141‑1 (MQ=11)
                                       tgcatAAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTaaca                                                                                                     <  1:237681/141‑1 (MQ=255)
                                          atAAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAAcc                                                                                                  <  2:134595/141‑1 (MQ=255)
                                           tAAACGCCCCTAACCACCTGGGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCt                                                                                                 <  2:209029/141‑1 (MQ=255)
                                                 tcccTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCCAACAACCTGAGGTCACAAGGCGTCCAATCGCCGACGCTGTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGtt                                                                                           >  2:102410/2‑140 (MQ=11)
                                                 ccccTAACAACCGGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTc                                                                                           <  1:220136/141‑1 (MQ=255)
                                                      aacaacCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTAc                                                                                      >  1:276983/1‑141 (MQ=255)
                                                        caacCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCGGAGGTTACTAGGc                                                                                    >  2:150695/1‑141 (MQ=255)
                                                          acCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCtt                                                                                  >  1:196260/1‑141 (MQ=255)
                                                            cTGAGGTCACTACGCCTCCGTTCCCCTGCGCTTTAGCGCATGCATAAAAGCCCCCAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTcc                                                                                <  2:107481/141‑1 (MQ=14)
                                                             tGAGGTCACTACGCCTCCGTGCGCCTGCGTTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTCTAGCGAATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCg                                                                               <  1:223591/141‑1 (MQ=14)
                                                                ggTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAAGAACCTGAGGTCACTACGCTTCCGTTc                                                                            >  1:175146/1‑141 (MQ=255)
                                                                ggTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCCGAGGTTACTGCGGTTCCGTTc                                                                            >  2:241142/1‑141 (MQ=14)
                                                                 gTCACTCCGCTTCGGTTCGGCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCg                                                                           <  1:123365/141‑1 (MQ=14)
                                                                 gTCACTACGCTCCCGGTCGCCGGCGCTTTTGCGCATGCATAAAAGCCCCTAACAACCTAAGATCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCg                                                                           <  1:158766/141‑1 (MQ=14)
                                                                  tCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGc                                                                          >  2:147151/1‑141 (MQ=255)
                                                                   cACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTGGcc                                                                         >  1:241569/1‑141 (MQ=255)
                                                                    aCTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCt                                                                        >  1:227438/1‑141 (MQ=255)
                                                                      cccgcTTCGGGTCGGCCGCGCTTTAGCGCATGCATAAAAGCCCCTACCAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTgc                                                                      <  1:11775/139‑1 (MQ=11)
                                                                       cacttcgcgTTCGCCGTGCTTTTGCGCATGCATAAAAGCCCCTAACCACCTGAGGTCACTGCGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTgcgc                                                                    <  2:46019/133‑1 (MQ=11)
                                                                        cgcTTCCGTTCCCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTCAGCTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTgcgc                                                                    <  1:222289/141‑1 (MQ=14)
                                                                        cgcTCCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTgcgc                                                                    <  2:189178/141‑1 (MQ=255)
                                                                         gcTTCGGTCCGCCTGCGTTTCTGCGCGTGCATAAAGGCTCCTAACAACCTGGGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCt                                                                   <  2:232160/141‑1 (MQ=14)
                                                                          ttgctgtgccccTGCGCTTTGCGTATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCttt                                                                 <  1:51625/131‑1 (MQ=11)
                                                                          cttgggTTCGCCCTGGCCTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCtt                                                                  <  1:163596/137‑1 (MQ=14)
                                                                          cTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCtt                                                                  <  1:130605/141‑1 (MQ=255)
                                                                           tCCGGTTCGCCCGCGCTGTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCttt                                                                 <  1:106333/141‑1 (MQ=255)
                                                                             ccgTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAg                                                               <  2:227438/139‑1 (MQ=255)
                                                                             cGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCGCTACGCGTCGGATCGCCTGCGCTTTat                                                               >  2:274366/1‑140 (MQ=11)
                                                                                 tCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCgca                                                           >  1:153708/1‑141 (MQ=255)
                                                                                 tCGCCTGCGCTTTAGCGCATGCAGAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCgca                                                           <  1:3307/141‑1 (MQ=11)
                                                                                 gcGCCTGCGCCTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCgca                                                           <  1:165115/140‑1 (MQ=11)
                                                                                    ccTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCgcatgc                                                        >  1:76440/1‑141 (MQ=255)
                                                                                      tctgcTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTCCGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCgcatgcat                                                      <  2:155351/138‑1 (MQ=255)
                                                                                      tGCGCTTTTGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCgcatgcat                                                      <  2:205030/141‑1 (MQ=255)
                                                                                      tGCGCTTATGCGCTTGCATAAAAGCCCCGAACCACCTGAGCGCACTACCCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAGGTCCCTAACAACCCGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCgcatgcat                                                      <  2:121890/141‑1 (MQ=14)
                                                                                          cTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATaaaa                                                  >  1:159394/1‑141 (MQ=255)
                                                                                                  catgcatAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCgcatgca                                                                                                                    <  2:254705/67‑1 (MQ=40)
                                                                                                  catgcatAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCgcatgca                                                                                                                    <  1:190888/67‑1 (MQ=255)
                                                                                                  cagtcatAAAAGCCCCGAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTaa                                          <  1:46652/137‑1 (MQ=255)
                                                                                                    tgcatAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTGCTGTTCGCCGGCGCGTTAGCGCAAGCAGGAAAGCCCCTaaca                                        >  2:256950/1‑141 (MQ=11)
                                                                                                     gcatAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTaacaa                                       <  2:226384/141‑1 (MQ=255)
                                                                                                      catAAACGCCCCTAACAACCTGAGGTCCCTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTCCGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTaacaac                                      <  1:98448/141‑1 (MQ=255)
                                                                                                       atAAAAGCCCCTAACAACCTGGGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCTCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAAcc                                     <  1:38081/141‑1 (MQ=255)
                                                                                                        tAAAAGCCCCTAACAACCTGAGGTTTCTGCGCCTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCCTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAAcct                                    >  1:68694/1‑140 (MQ=11)
                                                                                                         aaaaGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTCCGCTTCCGTTCGCCTTCGCTTTAGCCCATGCATAAAAGCCCCTAACAAcctg                                   >  1:12740/1‑139 (MQ=255)
                                                                                                           aggCCCCCAACAACCTGAGGTCACTCCGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCTTGCATAAAAGCCCCTTACAAcctttg                                 <  2:177398/139‑5 (MQ=14)
                                                                                                           aaGCCCCTCACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTTACAAcctttg                                 <  1:256950/141‑5 (MQ=255)
                                                                                                           aaGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATa                                                                                                                  <  1:122063/60‑1 (MQ=255)
                                                                                                              ccccTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAAcctgaggtc                              <  1:155779/141‑8 (MQ=255)
                                                                                                                ccTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGt                            <  1:80425/141‑1 (MQ=255)
                                                                                                                ccTAACAACCCGCGGCCACTACGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGt                            <  2:255126/141‑1 (MQ=14)
                                                                                                                 cTACCAACCTGGGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGGCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTa                           <  2:175146/141‑1 (MQ=255)
                                                                                                                   aacaacCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCCGAGGTCGCTACGCTTCCGTTCGGCTGCACATTAGCGCATCAATAGACGCCGCTGACAACCCACGGTTGTAAg                         >  2:287374/1‑141 (MQ=17)
                                                                                                                   aacaacCTGAGGGCACTACGCATCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAg                         >  1:148558/1‑141 (MQ=255)
                                                                                                                       acCTTACGTCACTTCCCTCCGCTTCGCCTGCGCTTTAGCGCATGCAGAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGc                     <  2:135653/141‑1 (MQ=14)
                                                                                                                         cggAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCtt                   <  1:241142/139‑1 (MQ=255)
                                                                                                                            aGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTACGGTTGTCACGGGCtcgtt                >  2:87998/1‑137 (MQ=14)
                                                                                                                             ggTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAGGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTAAGGTTgctacgggtttgttt               >  2:116740/1‑127 (MQ=14)
                                                                                                                             ggTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTaaa               <  1:213235/141‑1 (MQ=255)
                                                                                                                                 aCTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTc           >  1:54799/1‑141 (MQ=255)
                                                                                                                                    aCGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGa        <  2:76440/141‑1 (MQ=255)
                                                                                                                                          cGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGACCCGCTTTTCTTGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTAt  >  1:66435/1‑141 (MQ=21)
                                                                                                                                           |                                                                                                                                           
TTGAATAGCTCCACTTGCCGTTCGCCTGCGCTTTGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTGCGCTCCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCTGAGGTCACTACGCTTCGGTTCGCCTGCGCTTTAGCGCATGCATAAAAGTCCCTAACAACCTGAGGTCACTACGCTTCCGTTCGCCTGCGCTTTAGCGCATGCATAAAAGCCCCTAACAACCAAAGGTTGTAAGGGGCTTTAAATGTCAGATAGTAT  >  CP000730/2504818‑2505096

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