Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 394,213 G→A 7.4% V140V (GTC→GTT)  USA300HOU_0361 ← hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730394,2130GA 7.4% 211.1 / 4.0 81V140V (GTC→GTT) USA300HOU_0361hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base G (47/27);  new base A (2/4);  total (49/32)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.69e-01

ACCAAGCGCCATAAACATTTGTCGTCATTGAAAAGATAGCGTAGATAGCTGTATATATAAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGCATCATTACAATACTGAAAATAACGACACCTACATCGATGGCT  >  CP000730/394074‑394351
                                                                                                                                            |                                                                                                                                          
aCCAAGCGCCATAAACATTTGTCGTCATTGAAAAGATAGCGTAGATAGCTGTATATATAAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCg                                                                                                                                            >  1:500652/1‑140 (MQ=255)
   aaGCGCCATAAACATTTGTCGTCATTGAAAAGATAGCGTAGATAGCTGTATATATAAGTGAAAATATCGTACCAATCAATGCGATGATACAGATACTA‑ATAGCTGTAGCCTTTTAAAGTGATTTGAATATTGGAGCGACt                                                                                                                                         >  1:582560/1‑140 (MQ=255)
    aGCGCCATAAACATTTGTCGTCATTGAAAAGATAGCGTAGATAGCTGTATATATAATTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGTAGCGACTg                                                                                                                                        >  1:89388/1‑140 (MQ=255)
    aGCGCCATAAACATTTGTCGTCATTGAAAAGATAGCGTAGATAGCTGTATATATAAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTg                                                                                                                                        >  2:545654/1‑140 (MQ=255)
    aGCGCCATAAACATTTGTCGTCATTGAAAAGATAGCGTAGATAGCTGTATATATAAGTGAAAATATCGTACCAATCAATGCGATGATACCAATACTA‑ATAGCTGTATCCTTTTAACGTAATTTGAATATTTGAGAGACTg                                                                                                                                        >  1:139442/1‑140 (MQ=255)
    aGCGCCATAAACATTTGTCGTCATTGAAAAGATAGCGTAGATAGCTGTATATATAAGTAAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTg                                                                                                                                        >  1:42558/1‑140 (MQ=255)
      cgcCATAAACATTTGTCGTCATTGAAAAGATAGCGTAGATAGCTGTATATATAAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGATCGACTGCg                                                                                                                                      >  2:663075/1‑140 (MQ=255)
      cgcCATAAACATTTGTCGTCATTGAAAAGATAGCGTAGATAGCTGTATATATAAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACAGCg                                                                                                                                      >  2:512066/1‑140 (MQ=255)
            aaaCATTTGTCGTCATTGAAAAGATAGCGTAGATAGCTGTATATATAAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATa                                                                                                                                <  1:74271/140‑1 (MQ=255)
                   tGTCGTCATTGAAAAGATAGCGTAGATAGCTGTATATATAAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATagcagcag                                                                                                                         <  1:288319/140‑1 (MQ=255)
                      cgtcATTGAAAAGGTAGCGTAGATAGCTGTATATATAAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGcac                                                                                                                      <  1:147407/140‑1 (MQ=255)
                      cgtcATTGAAAAGATAGCGTAGATAGCTGTATATATAAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGcac                                                                                                                      <  1:147604/140‑1 (MQ=255)
                      cgtcATTGAAAAGATAGCGTAGATAGCTGTATATATAAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGcac                                                                                                                      <  1:475227/140‑1 (MQ=255)
                       gtcATTGAAAAGATAGCGTAGATAGCTGTATATATAAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGcacc                                                                                                                     >  2:600337/1‑140 (MQ=255)
                                   aTAGCGTAGATAGCTGTGTATCTGCGTGTAATTATCGTATCTATCTATGCGATGTTACCGATAGTA‑ATGGCTGTATCCTTGTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCtgt                                                                                                         <  1:78935/140‑1 (MQ=255)
                                   aTAGCGTAGATAGCTGTATATATAAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCtgt                                                                                                         <  1:359427/140‑1 (MQ=255)
                                   aTAGCGTAGATAGCTGTATATATAAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCtgt                                                                                                         >  2:25667/1‑140 (MQ=255)
                                   aTAGCGTAGATAGCTGTATATATAAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGAACCACAAACGCCtgt                                                                                                         >  1:282722/1‑140 (MQ=255)
                                    tatcGTAGATAGCTGTATATATAAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCtgtg                                                                                                        <  2:120049/137‑1 (MQ=255)
                                      gCGTAGATAGCTGTATATATAAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGACCGTGcc                                                                                                      >  2:218708/1‑140 (MQ=255)
                                      gCGTAGATAGCTGTATATATAAGTGAAAATATCGTACCAATCAATGCGATGATACAGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACACCTGTGcc                                                                                                      >  2:450615/1‑140 (MQ=255)
                                      gCGTAGATAGCTGTATATATAAATGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTGGAATATTGGACCGACTGCGGCAATAGCAGCAGCACCACAACCGCGCGTGcc                                                                                                      >  2:80360/1‑140 (MQ=255)
                                       cGTAGATAGCTGTATATATAAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCg                                                                                                     >  1:290429/1‑140 (MQ=255)
                                       cGTAGATAGCTGTATATATAAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGACAATAGCAGCAGCACCACAGACGCCTGTGcag                                                                                                     >  1:242833/1‑138 (MQ=255)
                                        gTAGATAGCTGTATATATAAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGACGATAGCAGCAGCACCACAAACGCCTGTGCCGa                                                                                                    >  1:357991/1‑140 (MQ=255)
                                            aTAGCTGTAAATATAAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAAGAGCAGCAGCACCACAAACGCCTGTGCCGACAcc                                                                                                <  1:97782/140‑1 (MQ=255)
                                                   tatatataAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTAGCATTTTCACGTGATTTGAATATTGGAGCGAATGCGGCAATAGAAGCAGCACCACAAACGCCTGTGCCGACCCCTAGtaat                                                                                         >  2:51548/1‑140 (MQ=255)
                                                     tatataAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGCTTTGAATATTGGAGCGACTGCGGCAATAGCCGCAGCACCCCAAACGCCCGTGCCCACACCCAGgaataa                                                                                       >  2:649018/1‑140 (MQ=255)
                                                       tataAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCATGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACATCTAGTAATAATg                                                                                     >  1:52241/1‑140 (MQ=255)
                                                       tataAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATg                                                                                     <  2:571577/140‑1 (MQ=255)
                                                        ataAGTAAAAATATCGTACCAACCAATGCGATGCTACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGc                                                                                    <  2:396645/140‑1 (MQ=255)
                                                         taAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTAAAGTGATATAAATATTGGAGCGACTGCGGCAATAGCAGCAGAACCACAAACGCCTGTGCCGACACCTAGTAATAATGCg                                                                                   >  1:333357/1‑140 (MQ=255)
                                                          aaGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGa                                                                                  >  1:545033/1‑140 (MQ=255)
                                                           aGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGAt                                                                                 >  1:239102/1‑140 (MQ=255)
                                                                aaaTATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCACTAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGtttt                                                                            >  2:566694/1‑140 (MQ=255)
                                                                aaaTATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCATCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGtttt                                                                            >  1:125486/1‑140 (MQ=255)
                                                                  atatCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTg                                                                          >  1:631197/1‑140 (MQ=255)
                                                                  atatCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTg                                                                          >  2:446077/1‑140 (MQ=255)
                                                                   tatCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGACGACACCTAGTAAAAATGCGATGTTTTTGt                                                                         >  1:321645/1‑140 (MQ=255)
                                                                    atCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTACGGACACCTAGTAATAATGCGATGTTTTTGTc                                                                        >  2:115802/1‑140 (MQ=255)
                                                                      cGTACCAATAAATACGATGATACCGA‑ACAACACATCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGGTTTTGTCAc                                                                      >  2:113305/1‑140 (MQ=255)
                                                                        tACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCAACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGTCGACACCTAGTAATAATGCGATGTTTTTGTCACCa                                                                    >  1:269417/1‑140 (MQ=255)
                                                                        tACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCAACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCa                                                                    >  1:610142/1‑140 (MQ=255)
                                                                         aCCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCAt                                                                   <  1:6539/140‑1 (MQ=255)
                                                                          cgatttcgttctatgatACCGAGACTA‑ATAGCTGCATGCTTTTCACTTGATTCGATTATTGTAGCAACTGCGGCTATAGCAGCAGCACCATAAACGCCGGTGCCGACACCTAGTAATAATGCGCTGTTTTTGTCACCATg                                                                  <  2:184472/128‑1 (MQ=255)
                                                                          ccAATCAATGCGATGATACCGATTCTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCAACTGCGGCAATAGCAGCAGCACCACAACCGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATg                                                                  <  2:588406/140‑1 (MQ=255)
                                                                                 atGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCCCCATGCACCAGt                                                           >  2:134017/1‑140 (MQ=255)
                                                                                 atGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGt                                                           >  2:531958/1‑140 (MQ=255)
                                                                                   gCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGttt                                                         >  2:9981/1‑140 (MQ=255)
                                                                                     gatgatACCGATACTA‑TTATCTGTATCCTTTTTACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCAACACAAACGCCTGTGCCGACACCTAGTAATGATGCGATGTTTTTGTCACCATGCAACAGTttgt                                                       <  1:648947/140‑1 (MQ=255)
                                                                                     gatgatACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTttgt                                                       <  1:86923/140‑1 (MQ=255)
                                                                                     gatgatACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTttgt                                                       >  1:402607/1‑140 (MQ=255)
                                                                                     gatgatACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTttgt                                                       <  2:500652/140‑1 (MQ=255)
                                                                                        gatACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTTTGTTGa                                                    <  1:600337/140‑1 (MQ=255)
                                                                                           aCCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTTTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTTTGTTGACaa                                                 <  2:70272/140‑1 (MQ=255)
                                                                                            ccGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGATGCTCCACAAACGCCTGTTCGGACACCGAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTTTGTTGACaaa                                                <  1:156466/140‑1 (MQ=255)
                                                                                                  cTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGcat                                          >  2:375006/1‑140 (MQ=255)
                                                                                                  cTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGaat                                          >  2:375037/1‑140 (MQ=255)
                                                                                                  cTA‑ATAGCTGTATCCTTTCCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGcat                                          >  2:573330/1‑140 (MQ=255)
                                                                                                  cTA‑ATAGCTGTATCATTTTCACGTGAATTGAATATTGGAGCGACTGCGGCAATAGCCGCAGCAGCACAAACGCCTGTGCCGAGACATAGTAATAATTCGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGcat                                          >  2:166682/1‑140 (MQ=255)
                                                                                                    a‑aTAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGcatca                                        <  1:566559/140‑1 (MQ=255)
                                                                                                      aTAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGcatcat                                       >  1:614425/1‑140 (MQ=255)
                                                                                                       tAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATACTGCGATGTTTTTGTCACCATGCAACAGTTTGTTAACAAAGAGCAACAtt                                      >  2:585386/1‑140 (MQ=255)
                                                                                                       tAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATCTTTTTGTCACCATGCAACAGTTTGTTGAGCAAGAGCATCAtt                                      >  2:423723/1‑140 (MQ=255)
                                                                                                         gCTGTATCCGTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGCATCATTAc                                    <  2:343907/140‑1 (MQ=255)
                                                                                                         gCTGAATCCTTTACACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGCATCATTAc                                    >  1:470636/1‑140 (MQ=255)
                                                                                                                ccTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGAGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGCATCATTACAATACTg                             <  1:519479/140‑1 (MQ=255)
                                                                                                                    ttCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGCATCATTACAATACTGaaaa                         >  1:641799/1‑140 (MQ=255)
                                                                                                                        cGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCTTGTGCCGTTACCTTGTAATAATGCGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGCATCATTACAATACTGAAAATAAc                     <  1:226511/140‑1 (MQ=255)
                                                                                                                         gTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGCATCATTAAAATACTGAAAATAACg                    >  2:618019/1‑140 (MQ=255)
                                                                                                                           gATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGCATCATTACAATACTGAAAATAACGac                  >  2:161319/1‑140 (MQ=255)
                                                                                                                           gATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGCATCATTACAATACTGAAAATAACGac                  >  2:44383/1‑140 (MQ=255)
                                                                                                                                gAATATTGGAGCAACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGTCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGCATCATTACAATACTGAAAATAACGACACCTa             <  2:269417/140‑1 (MQ=255)
                                                                                                                                gAATACTGGAGCAACTGCGGCAATAGCAGCAGCACCACAACCGCCTGTTCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGCTTCATTACAATACTGAAAATAACGACACCTa             <  2:610142/140‑1 (MQ=255)
                                                                                                                                 aaTTTTGGAGCGACTGCGGCAATAGCAGCAGCACCACACACGCCTGTGCCGACACCTATTAATAATGCGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGCATCATTACAATACTGAAAATAACGACACCTAc            <  2:208880/140‑1 (MQ=255)
                                                                                                                                   tatTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGCATCATTACAATACTGAAAATAACGACACCTACAt          <  1:115802/140‑1 (MQ=255)
                                                                                                                                    atTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGCATCATTACAATACTGAAAATAACGACACCTACATc         <  2:42558/140‑1 (MQ=255)
                                                                                                                                     ttGGAGCCACTGCGGAAGTTGCGGCAGCACCAGAAACGCCTCGGCCGTCACCTAGTAATAATGCAATTTTTTTCTCACCATGCAACAGTTTGTTGACAAAGTGCTTCATTACAATACTGAAAATAACGACACCTACATCg        <  2:437559/140‑1 (MQ=255)
                                                                                                                                           cGATTGCGGCAATAGCAGCAGCACCACAAACGCCTGTCCCTACACCTAGTAATAACGCGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGCATCATTACAATACTGAAAATAACGACACCTACATCGATGGCt  <  2:125486/140‑1 (MQ=255)
                                                                                                                                           cGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGCATCATTACAATACTGAAAATAACGACACCTACATCGATGGCt  >  2:288484/1‑140 (MQ=255)
                                                                                                                                           cGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAAAAATGAGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGCATCATTACAATACTGAAAATAACGACACCTACATCGATGGCt  >  1:19401/1‑140 (MQ=255)
                                                                                                                                           cGACTGCGGCAATACCAGCAGCACCACAAACGCCTGCGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGCATCATTACAATACTGAAAATAACGACACCTACATCGATGGCt  <  2:43516/140‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                          
ACCAAGCGCCATAAACATTTGTCGTCATTGAAAAGATAGCGTAGATAGCTGTATATATAAGTGAAAATATCGTACCAATCAATGCGATGATACCGATACTA‑ATAGCTGTATCCTTTTCACGTGATTTGAATATTGGAGCGACTGCGGCAATAGCAGCAGCACCACAAACGCCTGTGCCGACACCTAGTAATAATGCGATGTTTTTGTCACCATGCAACAGTTTGTTGACAAAGAGCATCATTACAATACTGAAAATAACGACACCTACATCGATGGCT  >  CP000730/394074‑394351

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