Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 607,884 G→A 89.8% R209* (CGA→TGA)  USA300HOU_0550 ← deoxynucleoside kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730607,8840GA89.8% 193.9 / 13.1 80R209* (CGA→TGA) USA300HOU_0550deoxynucleoside kinase
Reads supporting (aligned to +/- strand):  ref base G (3/5);  new base A (26/44);  total (29/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

ATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTCCAG  >  CP000730/607750‑608020
                                                                                                                                      |                                                                                                                                        
atGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAAt                                                                                                                                    <  2:725874/141‑1 (MQ=255)
  gtgAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAAttt                                                                                                                                  >  1:298691/1‑141 (MQ=255)
           tGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCa                                                                                                                         >  2:1012705/1‑141 (MQ=255)
           tGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTCTGTTTATCa                                                                                                                         >  1:362068/1‑141 (MQ=255)
            gTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTTTTTATCAt                                                                                                                        >  1:372008/1‑141 (MQ=255)
            gTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCAt                                                                                                                        >  1:732736/1‑141 (MQ=255)
             tATTCGATAGAATTGATAAAATAATCGATAATATCTGTGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCAATATATTTGAATATATAAAGAAATTTTGTTTAtataa                                                                                                                       >  1:829824/1‑137 (MQ=255)
             tATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAGTTTTGTTTTTCATa                                                                                                                       >  1:207657/1‑141 (MQ=255)
             tATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATa                                                                                                                       <  2:718130/141‑1 (MQ=255)
              aTTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAg                                                                                                                      >  2:416185/1‑141 (MQ=255)
              aTTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAg                                                                                                                      <  2:593450/141‑1 (MQ=255)
                    tAGAATTGATAAAATAATCGATTATATCTATGGAATATATATTGAATATTAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGGCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTa                                                                                                                <  2:393355/141‑1 (MQ=255)
                    tAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTa                                                                                                                >  1:536341/1‑141 (MQ=255)
                    tAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTa                                                                                                                <  1:842394/141‑1 (MQ=255)
                    tAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGAGTTa                                                                                                                >  1:1020184/1‑141 (MQ=255)
                     aGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTaa                                                                                                               <  1:230655/141‑1 (MQ=255)
                     aGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTaa                                                                                                               <  2:104513/141‑1 (MQ=255)
                      gAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCCAGCAATTTTGTTTATCATAGGATTTaaa                                                                                                              <  1:887371/141‑1 (MQ=255)
                        aTTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGa                                                                                                            <  2:557283/141‑1 (MQ=255)
                          tGATAACATAATCGATAATATCTATTGAATATACATTGAATATACACATATATGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAACCAATTTTGTTTATCATAGTATTTAAAGAtt                                                                                                          <  2:964279/141‑1 (MQ=255)
                          tGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGAtt                                                                                                          <  2:52365/141‑1 (MQ=255)
                           gATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTc                                                                                                         >  2:56287/1‑141 (MQ=255)
                             tAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCgg                                                                                                       <  2:585544/141‑1 (MQ=255)
                             aaaaaTAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCgg                                                                                                       <  2:310317/140‑1 (MQ=255)
                               aaaTAATCTATAATAGCTATGGAATGTACATTGAATAAAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTTTCTACTTGTCGATATGTTTCAATAATTCCACCAATTTTGTCTATCATATGATTTAAAGATTCgggg                                                                                                     <  2:658402/141‑1 (MQ=255)
                                   aatCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCt                                                                                                 <  1:708852/141‑1 (MQ=255)
                                    atCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTCTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGTGGTCCtt                                                                                                >  2:325390/1‑141 (MQ=255)
                                     tCGATAATAACTATGGAATATACCTTGGATAAAAACATGTACGCTATGTCTTTAGTCTTTTATCGTATATCTACTTGTCGATGTGTTTGAATAATTCAAGCAATTTTGTTTATCGTAGGATTTAAAGATTCGGGGTCCttt                                                                                               >  1:779542/1‑140 (MQ=255)
                                          aataTTTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGCCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGat                                                                                          <  2:974698/141‑1 (MQ=255)
                                          aataTCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGat                                                                                          <  2:298691/141‑1 (MQ=255)
                                             atCTATCGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATc                                                                                       <  2:1037744/141‑1 (MQ=255)
                                              tCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCa                                                                                      <  2:464510/141‑1 (MQ=255)
                                                tATGGAATATACATTGAATATAAACATATACGTTATGTCTTTAGTCTTTTACCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCata                                                                                    <  1:26815/141‑1 (MQ=255)
                                                 aTGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCatat                                                                                   <  1:163900/141‑1 (MQ=255)
                                                    gAATATAAATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCa                                                                                <  1:403426/141‑1 (MQ=255)
                                                      atatACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCATGATGGATATCATATTCAtt                                                                              >  2:264486/1‑141 (MQ=255)
                                                      atatACATTGAATATAAACATAGACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAAGTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCAtt                                                                              <  1:23537/141‑1 (MQ=255)
                                                        ataCATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTaa                                                                            <  1:967470/141‑1 (MQ=255)
                                                          aCATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAata                                                                          <  2:268683/141‑1 (MQ=255)
                                                          aCATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAata                                                                          <  2:390849/141‑1 (MQ=255)
                                                           cATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTAGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAatat                                                                         >  1:1055448/1‑141 (MQ=255)
                                                               gAATATTACCCTATACGTTATGTTTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGAt                                                                     <  2:372008/141‑1 (MQ=255)
                                                                 atatTAACATATACGCTATGTCTTTAGTCTTTTATCGTTTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATAc                                                                   <  1:116915/141‑1 (MQ=255)
                                                                     aaaCATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTAc                                                               >  2:824332/1‑141 (MQ=255)
                                                                      aaCATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACa                                                              >  1:1043588/1‑141 (MQ=255)
                                                                      aaCATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACa                                                              >  1:25039/1‑141 (MQ=255)
                                                                       aCATATACGCTATGTCGTTAGTGCTTTATCGCGTATCAACTTGTCGATATGTTTGAATAATACAAACAATTTTGTTTATCATAGGATGTAAAGATTCGGGGTCTATATGGATATCATATTCATTAGTATTGATACGTAcca                                                             >  2:672/1‑139 (MQ=255)
                                                                             aCGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGAc                                                       <  1:796499/141‑1 (MQ=255)
                                                                             aCGCTATGACTTTAGTCTTTTATCGTGTATGTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTCCAATTGGAc                                                       >  2:380718/1‑141 (MQ=255)
                                                                              cgccgtggcgttAGGCTTTTTTCGTGTATTTACTTATCGATATGTTGGAATAATTCAAGCAATTTTGTTTATCATAGGAATTTAAGTTTCGGGGTCCTTATGGATATCGTATTCATTAATATTGATACGTTGAACTGGACa                                                      <  2:238501/131‑1 (MQ=255)
                                                                              cGCTATGTCTTTAGTTGTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACa                                                      <  2:328966/141‑1 (MQ=255)
                                                                              cGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACa                                                      <  2:1041716/141‑1 (MQ=255)
                                                                              cGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACa                                                      <  1:10515/141‑1 (MQ=255)
                                                                              cGCTATGTCTTTAGTCTTTTATCGGGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACa                                                      <  1:466358/141‑1 (MQ=255)
                                                                                    gTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCAtt                                                >  1:9911/1‑141 (MQ=255)
                                                                                       tttAGTTTTTTATCGTTTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTaaa                                             <  1:419144/141‑1 (MQ=255)
                                                                                        ttAGTCTTTTATGGTGTACCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAg                                            <  1:927811/141‑1 (MQ=255)
                                                                                        ttAGTCTTTTATCGTGTTTCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATTGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAg                                            <  1:1050284/141‑1 (MQ=255)
                                                                                        ttAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAg                                            <  2:477343/141‑1 (MQ=255)
                                                                                          aGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCt                                          <  2:330156/141‑1 (MQ=255)
                                                                                          aGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCt                                          <  2:839561/141‑1 (MQ=255)
                                                                                          aGTCTTTTATCGGGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCt                                          <  1:178726/141‑1 (MQ=255)
                                                                                               tttATCGTGTATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTaa                                     <  2:325166/141‑1 (MQ=255)
                                                                                                  ttCGTGTTTCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATcc                                  <  2:362068/140‑1 (MQ=255)
                                                                                                       gtATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCg                             <  2:242137/141‑1 (MQ=255)
                                                                                                       gtATCTACTTGTCGATATGCTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCg                             >  1:635567/1‑141 (MQ=255)
                                                                                                        tATCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGt                            <  2:418653/141‑1 (MQ=255)
                                                                                                         aTCTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTc                           <  1:222357/141‑1 (MQ=255)
                                                                                                           cTACTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCAt                         <  2:614537/141‑1 (MQ=255)
                                                                                                             aCTTGTCGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAg                       <  2:463096/141‑1 (MQ=255)
                                                                                                             aCTTGACGATATGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAg                       <  2:770992/141‑1 (MQ=255)
                                                                                                               ttgtaaatatGGTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCg                     <  2:1010060/135‑1 (MQ=255)
                                                                                                               ttGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAgag                     >  2:707547/1‑139 (MQ=255)
                                                                                                                      tatGTTTGAATAATTCAGGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAAt              >  1:183974/1‑141 (MQ=255)
                                                                                                                        tGTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAg            >  2:981840/1‑141 (MQ=255)
                                                                                                                         gTTTGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGc           <  1:646464/141‑1 (MQ=255)
                                                                                                                          tttGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGCTATCATATTCATTAATATTGAGACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGAGTTTAAATAGCt          >  2:362034/1‑141 (MQ=255)
                                                                                                                          tttGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGGCATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCt          >  2:362141/1‑141 (MQ=255)
                                                                                                                           ttGAATAATTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCCTATTCATTACTATTGATCCGTACAACTGGACATGCATTAAAGCTATTCATCCAGCCGTCATACCGTTTCACTAACtt         >  2:442602/1‑141 (MQ=255)
                                                                                                                                  atTCAAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGCCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCTTCATAGCGTTTAAATATCTTTTTCCAg  >  1:928011/1‑141 (MQ=255)
                                                                                                                                      |                                                                                                                                        
ATGTGAATTAATGTATTCGATAGAATTGATAAAATAATCGATAATATCTATGGAATATACATTGAATATAAACATATACGCTATGTCTTTAGTCTTTTATCGTGTATCTACTTGTCGATATGTTTGAATAATTCGAGCAATTTTGTTTATCATAGGATTTAAAGATTCGGGGTCCTTATGGATATCATATTCATTAATATTGATACGTACAACTGGACATGCATTAAAGCTATTAATCCAATCGTCATAGCGTTTAAATAGCTTTTTCCAG  >  CP000730/607750‑608020

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