Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 1,419,527 G→T P7T (CCT→ACT)  USA300HOU_1317 ← hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,419,5270GT98.6% 222.4 / ‑2.9 71P7T (CCT→ACT) USA300HOU_1317hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (38/31);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

TTTGGAATAAAAAATGCAAGCACACCTAATATAACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGGCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAATATTCTATCTACTAATAGAAA  >  CP000730/1419394‑1419656
                                                                                                                                     |                                                                                                                                 
tttGGAATAAAAAATGCAAGCACACCTAATATAACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCg                                                                                                                            <  1:557829/141‑1 (MQ=255)
tttGGAATAAAAAATGCAAGCACACCTAATATAACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCg                                                                                                                            >  2:418063/1‑141 (MQ=255)
tttGGAATAAAAAATGCAAGCACACCTAATATAACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCg                                                                                                                            <  1:773370/141‑1 (MQ=255)
 ttGGAATAAAAAATGCAAGCACACCTAATATAACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGa                                                                                                                           >  2:258919/1‑141 (MQ=255)
  tGGAATAAAAAATGCAAGCACACCTAATATAACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGaa                                                                                                                          >  2:353404/1‑141 (MQ=255)
  tGGAATAAAAAATGCAAGCACACCTAATATAACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGaa                                                                                                                          >  1:561763/1‑141 (MQ=255)
  cggAATAAAAAATGCAAGCACACCTAATATAACAAATGTAATTACTGCAGATACAAACGATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCGATAAATAATTAGTCATATCGaa                                                                                                                          >  2:283946/2‑141 (MQ=255)
   ggAATAAAAAATGCAAGCACACCTAATATAACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGaat                                                                                                                         >  2:28735/1‑141 (MQ=255)
    gAATAAAAAATGCAAGCACACCTAATATAACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGaata                                                                                                                        <  2:91518/141‑1 (MQ=255)
    gAATAAAAAATGCAAGCACACCTAATATAACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGaata                                                                                                                        <  1:683846/141‑1 (MQ=255)
    gAATAAAAAATGCAAGCACACCTAATATAACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGaata                                                                                                                        >  2:807335/1‑141 (MQ=255)
      aTAAAAAATGCAAGCACACCTAATATAACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGaataat                                                                                                                      >  2:332946/1‑141 (MQ=255)
      aTAAAAAATGCAAGCACACCTAATATAACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGaataat                                                                                                                      >  1:289432/1‑141 (MQ=255)
          aaaaTGCAAGCACACCTAATATAACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCa                                                                                                                  <  1:286709/141‑1 (MQ=255)
                cAAGCACACCTAATATAACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATattca                                                                                                            <  2:413539/141‑1 (MQ=255)
                 aaGCACACCTAATATAACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATattcat                                                                                                           >  1:461847/1‑141 (MQ=255)
                  aGCACACCTAATATAACAAATGCAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTTATAAAATCATTATTAATGTTGTGATCCATAATTAATTATTCATATCGAATAATGTCATattcatt                                                                                                          >  2:28278/1‑141 (MQ=255)
                   gCACACCTAATATAACAAATGTAATCACTGGAGATACAAACCAGTTATATAAAACTAAGCAACAGAATATGGTTAATAAAATCATTATTAATGTGGTGATCCATAAATAATTAGTAATATCGAATAATGTCATattcatta                                                                                                         >  2:151581/1‑140 (MQ=255)
                           aaTATAACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTa                                                                                                 >  2:39213/1‑141 (MQ=255)
                               taACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTc                                                                                             >  2:663305/1‑141 (MQ=255)
                               taACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTc                                                                                             <  1:609038/141‑1 (MQ=255)
                                aaCAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCa                                                                                            >  1:843867/1‑141 (MQ=255)
                                  cAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCAtt                                                                                          >  2:348817/1‑141 (MQ=255)
                                  cAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCAtt                                                                                          >  1:844040/1‑141 (MQ=255)
                                    aaTGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTAc                                                                                        >  1:283061/1‑141 (MQ=255)
                                        tAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTc                                                                                    >  1:137272/1‑141 (MQ=255)
                                            tACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTg                                                                                >  1:172557/1‑141 (MQ=255)
                                              cTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGta                                                                              <  1:767425/141‑1 (MQ=255)
                                                gCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGtata                                                                            <  2:617966/141‑1 (MQ=255)
                                                 cAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATAc                                                                           <  2:775/141‑1 (MQ=255)
                                                      aCAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACAtttt                                                                      >  1:494152/1‑141 (MQ=255)
                                                        aaaCCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTttat                                                                    <  1:397707/141‑1 (MQ=255)
                                                          aCCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTttatta                                                                  >  2:396869/1‑141 (MQ=255)
                                                             atttatttAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATaa                                                               >  2:815059/1‑141 (MQ=255)
                                                               ttatttaAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAAt                                                             <  1:749144/141‑1 (MQ=255)
                                                                 atttaAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAAttt                                                           >  2:146636/1‑141 (MQ=255)
                                                                  tttaAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAAtttt                                                          <  2:503438/141‑1 (MQ=255)
                                                                     aaaaCTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGGCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATATATaa                                                       <  2:311456/141‑1 (MQ=255)
                                                                        aCTAAGCAACAGAATATTGTCAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTaaaaa                                                    <  1:497737/141‑1 (MQ=255)
                                                                         cTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAAc                                                   <  1:247061/141‑1 (MQ=255)
                                                                          tAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACt                                                  <  1:258919/141‑1 (MQ=255)
                                                                            aGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTa                                                >  1:187030/1‑141 (MQ=255)
                                                                              cAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTaaa                                              >  1:527223/1‑141 (MQ=255)
                                                                               aaCAGACTATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTaacc                                             <  2:608678/141‑3 (MQ=255)
                                                                               aaCAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATAAGAAATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTaaaa                                             >  1:235256/1‑140 (MQ=255)
                                                                                aCAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACa                                            >  1:322444/1‑141 (MQ=255)
                                                                                aCAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACa                                            >  2:236780/1‑141 (MQ=255)
                                                                                aCAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACa                                            <  2:41967/141‑1 (MQ=255)
                                                                                      tatTGTTAATAAAATCATTATTAATGTTGTGACCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGca                                      <  2:691256/141‑1 (MQ=255)
                                                                                        ttGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGcagt                                    <  2:691032/141‑1 (MQ=255)
                                                                                        ttGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGcagt                                    <  2:494152/141‑1 (MQ=255)
                                                                                            taataaAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGcagtcagt                                <  1:627707/141‑1 (MQ=255)
                                                                                              ataaAATCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGttt                              <  1:348817/141‑1 (MQ=255)
                                                                                                  aaTCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAg                          <  2:238672/141‑1 (MQ=255)
                                                                                                  aaTCATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAg                          <  2:328615/141‑1 (MQ=255)
                                                                                                     cATTATTAATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCaa                       <  2:325188/141‑1 (MQ=255)
                                                                                                        tattaATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAatat                    <  2:361264/141‑1 (MQ=255)
                                                                                                         attaATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAATAtt                   <  2:391268/141‑1 (MQ=255)
                                                                                                          ttaATGTTGTTATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAATATtc                  <  2:857431/141‑1 (MQ=255)
                                                                                                           taATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAATATtct                 >  2:308959/1‑141 (MQ=255)
                                                                                                           taATGTTGTGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAATATtct                 >  1:767174/1‑141 (MQ=255)
                                                                                                                ttgtGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTaaaaa                                                    >  1:551424/1‑101 (MQ=255)
                                                                                                                ttgtGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTATATTTCATTACTTTACTTGTATAAATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTAAAGCAATATtctatcaa            >  2:109543/1‑139 (MQ=255)
                                                                                                                ttgtGATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTACTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTaaaaa                                                    <  2:551424/101‑1 (MQ=255)
                                                                                                                   tgATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAATATTCTATctacta         >  2:332801/1‑141 (MQ=255)
                                                                                                                   tgATCCATAAATAATTAGTCATATCGAATAATGTCATATTCATGCTCCTTTTATTTGAGTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTATACAATAGCCGTCTGTTTCAAGCAATATTCTATctacga         >  2:301019/1‑139 (MQ=255)
                                                                                                                      tCCATAAATAATTAGTCATATCGAATAATTTCATATTCATTCTCCTTTTATATAATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAATATtctaaaaaaaaaaa      <  2:561763/141‑11 (MQ=255)
                                                                                                                      tCCATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAATATTCTATCTACTAATa      >  1:749171/1‑141 (MQ=255)
                                                                                                                        cATAAATAATTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAATATTCTATCTACTAATAGa    >  2:611050/1‑141 (MQ=255)
                                                                                                                         ataaataaTTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAATATTCTATCTACTAATAGaa   >  2:240189/1‑141 (MQ=255)
                                                                                                                          taaataaTTAGTCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAATATTCTATCTACTAATAGaaa  >  2:485285/1‑141 (MQ=255)
                                                                                                                                     |                                                                                                                                 
TTTGGAATAAAAAATGCAAGCACACCTAATATAACAAATGTAATTACTGCAGATACAAACCATTTATTTAAAACTAAGCAACAGAATATTGTTAATAAAATCATTATTAATGTTGTGATCCATAAATAATTAGGCATATCGAATAATGTCATATTCATTCTCCTTTTATTTCATTACTTTCCTTGTATACATTTTATTATAAATTTTTAAAAACTTAAACAATAGCAGTCAGTTTCAAGCAATATTCTATCTACTAATAGAAA  >  CP000730/1419394‑1419656

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