Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,005,078 T→C 5.8% I508T (ATT→ACT)  USA300HOU_0958 → M03 family oligopeptidase F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,005,0780TC 5.8% 205.6 / 3.2 69I508T (ATT→ACT) USA300HOU_0958M03 family oligopeptidase F
Reads supporting (aligned to +/- strand):  ref base T (22/43);  new base C (2/2);  total (24/45)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

AATGTTCGCAGAATTTGAGCATAAAATTCATGCAATTGAAGAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCTTAAGTCATCAAATTTTAACAGAAGGTAAGCCAGCAGTAGATAGATATATTAAT  >  CP000730/1004941‑1005217
                                                                                                                                         |                                                                                                                                           
aaTGTTCGCAGAATTTGAGCATAAAATTCATGCAATTGAAGAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAg                                                                                                                                          <  1:133637/141‑1 (MQ=255)
  tGTTCGCAGAAGTTGAGCATAAAATTCATGCAATTGAAGAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTa                                                                                                                                        <  1:391044/141‑1 (MQ=255)
     tCGCAGAATTTGAGCATAAAATTCATGCAATTGAAGAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATAGTTCGGTGACTCAGTGTTATTGGATGAAGATATTTGTAAgg                                                                                                                                     <  1:510680/141‑1 (MQ=255)
     tCGCAGAATTTGAGCATAAAATTCATGCAATTGAAGAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAgg                                                                                                                                     <  1:383067/141‑1 (MQ=255)
     tCGCAGAATTTGAGCACAAAAGTCAAGCAATTGAAGAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAgg                                                                                                                                     <  1:179696/141‑1 (MQ=255)
       gCAGAATTTGAGCATAAAATTCATGCAATTGAAGAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTCAGGaa                                                                                                                                   >  2:399590/1‑141 (MQ=255)
            aTTTGAGCATAAAATTCATGCAATTGAAGAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTc                                                                                                                              >  1:561940/1‑141 (MQ=255)
               tGAGCATAAAATTCATGCAATTGAAGAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACg                                                                                                                           >  1:409778/1‑141 (MQ=255)
               tGAGCATAAAATTCATGCAATTGAAGAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACg                                                                                                                           <  2:8644/141‑1 (MQ=255)
               tGAGCATAAAATTCATGCAATTGAAGAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATACTAGTAAGGAATGGTCACg                                                                                                                           >  2:195985/1‑141 (MQ=255)
               tGAGCATAAAATTCATGCAATTGAAGAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATACTAGTAAGGAATGGTCACg                                                                                                                           >  2:637924/1‑141 (MQ=255)
                  gCATAAAATTCATGCAATTGAAGAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTAt                                                                                                                        <  1:672953/141‑1 (MQ=255)
                    aTAAAATTCATGCAATTGAAGAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTc                                                                                                                      <  2:304432/141‑1 (MQ=255)
                      aaaaGTCATGCAATTGAAGAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCca                                                                                                                    <  2:442064/141‑1 (MQ=255)
                       aaaCTCATGCAATTGAAGCAGCGGATGTACCATTAACGCCAACGAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCcac                                                                                                                   <  2:660534/141‑1 (MQ=255)
                          ttCCTGCAATGGAAGAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACt                                                                                                                <  1:287942/141‑1 (MQ=255)
                           tCAGGCAATTGAATAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACtt                                                                                                               <  2:276879/141‑1 (MQ=255)
                            cATGCAATTGAAGAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTc                                                                                                              <  2:188068/141‑1 (MQ=255)
                            cATGCAAGTGAAGAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTTATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTc                                                                                                              <  2:356998/141‑1 (MQ=255)
                             aTGCAATTGAAGAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCAAACTTCt                                                                                                             >  2:340543/1‑141 (MQ=255)
                              tGCAATTGAAGAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCta                                                                                                            <  2:617594/141‑1 (MQ=255)
                              tGCAATTGAAGAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCta                                                                                                            <  2:653933/141‑1 (MQ=255)
                              tGCAATTGAAGAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCta                                                                                                            >  2:673688/1‑141 (MQ=255)
                                  aTTGAAGAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGCAAATGCCTGGGGAAAGAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTTCACACTTCTATATg                                                                                                        <  1:364521/141‑1 (MQ=255)
                                     gaagaaGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAat                                                                                                     <  1:594383/141‑1 (MQ=255)
                                     gaagaaGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAat                                                                                                     <  2:358279/141‑1 (MQ=255)
                                      aagaagCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATACTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAatt                                                                                                    <  1:637924/141‑1 (MQ=255)
                                      aagaagCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATACTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAatt                                                                                                    <  1:195985/141‑1 (MQ=255)
                                         aagCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGGAAGGAATGGTCGCGTATTACACACTTCaatcaaaaatata                                                                                                 >  1:315321/1‑131 (MQ=255)
                                          agCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAAttatta                                                                                                <  1:75333/141‑1 (MQ=255)
                                            cAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATg                                                                                              >  1:59806/1‑141 (MQ=255)
                                            cAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATg                                                                                              >  2:414600/1‑141 (MQ=255)
                                                           cgtcaacgtgtatggatgaAGATTTTGCCAAATTATATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCaa                                                                               <  2:249453/130‑1 (MQ=255)
                                                                  tAGAATGAATGAGGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAGACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTa                                                                        >  1:268398/1‑141 (MQ=255)
                                                                  tAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAGCTGGTTa                                                                        >  1:110276/1‑141 (MQ=255)
                                                                  tAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTa                                                                        >  1:659966/1‑141 (MQ=255)
                                                                  tAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATAAATACGCCACAGGTTa                                                                        >  2:131738/1‑141 (MQ=255)
                                                                  tAGAATGAATGAAGAAGATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGCATTATTATGTATATCAATACACAACTGGTTa                                                                        >  1:110239/1‑141 (MQ=255)
                                                                       tgaatgaaGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTCCAGTg                                                                   >  1:667220/1‑141 (MQ=255)
                                                                        gaatgaaGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGc                                                                  <  1:590339/141‑1 (MQ=255)
                                                                          atgaaGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTGCATTGCAg                                                                >  1:613772/1‑141 (MQ=255)
                                                                          atgaaGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAg                                                                <  1:170903/141‑1 (MQ=255)
                                                                           tgaaGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGc                                                               <  2:409778/141‑1 (MQ=255)
                                                                             aagaaTATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTc                                                             <  2:513437/141‑1 (MQ=255)
                                                                             aagaaTATGCCAAATTAAATAAATAATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTc                                                             >  1:506941/1‑141 (MQ=255)
                                                                              agaaTATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACCGTGCAGCCCa                                                            >  2:28198/1‑141 (MQ=255)
                                                                                    tGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATCTTCGCAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCtt                                                      >  1:334088/1‑141 (MQ=255)
                                                                                          aTTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCTTAAGtca                                                <  2:571142/141‑1 (MQ=255)
                                                                                                    atattCTTCGGTGATTCTGTAGAAACTGATGAAGATATGAGTAAGGTATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGAAACTTCTTACCGGGCAGTTCAAAGCTTAAGTCATCAAATTTTa                                      <  2:252684/140‑1 (MQ=255)
                                                                                                          tgctgctaTTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCTTAAGTCATCAAATTTTAACAGaa                                <  2:307710/134‑1 (MQ=255)
                                                                                                             ggTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATTTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCTTAAGTCATCAAATTTTAACAGAAGGt                             <  2:498812/141‑1 (MQ=255)
                                                                                                               tGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCTCGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCTTAAGTCATCAAATTTTAACAGAAGGTaa                           <  1:470333/141‑1 (MQ=255)
                                                                                                                   tCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCGCGTATTCCACACTTCTGTATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCTTAAGCCATCAACTTTTAACACAAGGTAAGAca                       >  2:400731/1‑141 (MQ=255)
                                                                                                                   tCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCTTAAGTCATCAAATTTTAACAGAAGGTAAGCca                       <  1:673688/141‑1 (MQ=255)
                                                                                                                      gTAGAAACTGATGAAGATATTAGTAAGTAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCTTAAGTCATCAAATTTTAACAGAAGGTAAGCcagca                    <  1:642912/141‑1 (MQ=255)
                                                                                                                      gTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCTTAAGTCATCAAATTTTAACAGAAGGTAAGCcagca                    <  1:593655/141‑1 (MQ=255)
                                                                                                                      gTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCTTAAGTCATCAAATTTTAACAGAAGGTAAGCcagca                    <  1:391553/141‑1 (MQ=255)
                                                                                                                      gTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCTTAAGTCATCAAATTTTAACAGAAGATAAGCcagca                    <  1:503286/141‑1 (MQ=255)
                                                                                                                      gTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCTTAAGTCAACCAATAATAACAGAAGGTAAGCcagca                    >  2:518055/1‑141 (MQ=255)
                                                                                                                       tAGAAACTGATGAAGATATTAGTCAGGAATGCGCACGGATTCCACACCTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCTTAAGTCATCAAATTTTAACAGACGGTAAGCcagcag                   <  2:315265/141‑1 (MQ=255)
                                                                                                                       tAGAAACTGATGAAGATATTAGTAAGGTCAAGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCTTAAGTCATCAAATTTTAACAGAAGGTAAGCcagcag                   <  2:315321/141‑1 (MQ=255)
                                                                                                                       gaGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCTTAAGTCATCAAATTTTAACAGAAGGTAAGCcagcag                   <  2:281501/140‑1 (MQ=255)
                                                                                                                       gaGAAACCGATGAAGATATTAGTAAGGAATGGGCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCTTAAGTCATCAAATTTTAACAGAAGGTAAGCcagcag                   <  2:422004/140‑1 (MQ=255)
                                                                                                                        aGAAACTGATGAAGATATTAGTAAGGGATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCTTAAGTCATCAAATTTTAACAGAAGGTAAGCCAGCAGt                  <  2:66564/141‑1 (MQ=255)
                                                                                                                           aaCTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCTTAAGTCATCAAATTTTATCAGAAGGTAAGCCAGCAGtagg               <  2:506941/141‑2 (MQ=255)
                                                                                                                              tgatgaAGATATTAGTAAGGAATGGTCACGTATTCCACATTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCTTAAGTCATCTAAGTTTAACAGAAGGTAAGCCAGCAGtagatag            <  2:29792/141‑1 (MQ=255)
                                                                                                                              tgatgaAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCTTAAGTCATCAAATTTTAACAGAAGGTAAGCCAGCagaagaaaa            >  2:605748/1‑138 (MQ=255)
                                                                                                                                 tgaAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCTTAAGTCATCAAATTTTAACAGAAGGTAAGCCAGCAGTAgatagata         >  1:397994/1‑141 (MQ=255)
                                                                                                                                  taaGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCTTAAGTCATCAAATTTTAACAGAAGGTAAGCCAGCAGTAGATAGatat        <  2:20350/140‑1 (MQ=255)
                                                                                                                                     aatatTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATTTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCTTAAGTCATCAAATTTTAACAGAAGGTCAGCCAGCTGTAGATAGATATAtt     <  2:669594/140‑1 (MQ=255)
                                                                                                                                        atTAGTAAGTAATGGTCACGTATTCCCCACTTCGCTATGTATTACTATGTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGTTTAAGTCATCAAATTTTAACAGAAGGTAAGCCAGCAGTAGATAGATATATTAAt  <  2:537908/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                           
AATGTTCGCAGAATTTGAGCATAAAATTCATGCAATTGAAGAAGCAGGTGAACCATTAACGCCAACTAGAATGAATGAAGAATATGCCAAATTAAATAAATTATACTTCGGTGATTCTGTAGAAACTGATGAAGATATTAGTAAGGAATGGTCACGTATTCCACACTTCTATATGAATTATTATGTATATCAATACGCAACTGGTTACAGTGCAGCTCAAAGCTTAAGTCATCAAATTTTAACAGAAGGTAAGCCAGCAGTAGATAGATATATTAAT  >  CP000730/1004941‑1005217

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