Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 250,944 C→A 8.8% Q109K (CAA→AAA)  USA300HOU_0224 → possible dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730250,9440CA 8.8% 178.9 / 2.3 65Q109K (CAA→AAA) USA300HOU_0224possible dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (35/23);  new base A (2/5);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTGTTAAAAGATGACACGATTGATGTTATCCATGTTTGTACGCCAAATGACTCGCATTGTGAAATTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGGTTATCAAAATCGTTTCCGAGCAGATAGTCAATTTTTACATCAAGCAGCGCAACGTGGCGACTTAGGAGACATTTACTTCGGAAAGGC  >  CP000730/250804‑251080
                                                                                                                                            |                                                                                                                                        
cATTGTTAAAAGATGACACGATTGATGTTATCCATGTTTGTACGCCAAATGACTCGCATTGTGAAATTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTc                                                                                                                                          <  1:595549/141‑1 (MQ=255)
   tGTTAAACGATGACACGGTTGATGTTATCCATGTTTGTACGACAAATGACTCGCATTGTGAAATGACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTCTCACACACAGGCGGCTAGACAGACAGCAGAAGCGCaaa                                                                                                                                       >  2:485048/1‑141 (MQ=255)
   tGTTAAAAGATGACAAGCGTTTTGTTATCCATGTTTGTACGCCAAATGACTCGCATTGTGAAATTACTGTAGCAGGGTTGCAAGCTGGTAAACATGTGATATGTGAAAAACCAATGGCTAAAACGACAACAGAAGCTCaaa                                                                                                                                       >  2:493947/1‑141 (MQ=255)
    gTTAAAAGATGACACGATTGATGTTATCCATGTTTGTACGCCAAATGACTCGCATTGTGAAATTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCaaaa                                                                                                                                      >  1:139883/1‑141 (MQ=255)
    gTTAAAAGATGACACGATTGATGTTATCCATGTTTGTACGCCAAATGACTCGCATTGTGAAATTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGAGAAAAACCAATGGCTAAAACGACAGCAGAAGCTAaaaa                                                                                                                                      >  2:136638/1‑141 (MQ=255)
             tGACACGATTGATGTTATCCATGTTTGTACGCCAAATGACTCGCATTGTGACATTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAAGGGCTAAAACGACGGCAGAAGCTCAAAAAATgataga                                                                                                                             >  1:403466/1‑141 (MQ=255)
             tGACACGATTGATGTTATCCATGTTTGTACGCCAAATGACTCGCATTGTGAAATTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATgataga                                                                                                                             >  1:484028/1‑141 (MQ=255)
             tGACACGATTGATGTTATCCATGTTTGTACGCCAAATGACTCGCATTGTGAAATTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAAAgataga                                                                                                                             >  2:226110/1‑141 (MQ=255)
             tGACACGATTGATGTTATCCATGTTTGTACACCAATTGACTCGCATGGTGAAATTACTGTAGCAGGGTTGGATGCTAGTAAACATGTGATGTGTGGAAGACCAATGGCTAAAAGGACAGCAGAATCTCAAAAAATaataga                                                                                                                             >  2:185415/1‑141 (MQ=255)
                cacGTTTGATGTTATCCAGGATTGTACGATAAGTGACTAGCGTTGTGAAATTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGAAAAAAAAAAAAAAATGATAGATAc                                                                                                                          <  2:341022/141‑1 (MQ=255)
                 acGATTGATGTTATCCATGTTTGTACGCCAAATGACTCGCATTGTGAAATTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACa                                                                                                                         >  2:30994/1‑141 (MQ=255)
                  cGATTGATGTTATCCATGTTTGTACGCCAAATGACTCGCATTGTGAAATTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAg                                                                                                                        <  1:264012/141‑1 (MQ=255)
                   gATTGATGTTATCCATGTTTGTACGCCAAATGACTCGCATTGTGAAATTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGc                                                                                                                       <  1:47646/141‑1 (MQ=255)
                      tGAAGTTATCCAGGTTTGTACGCCAAATGACTCGCATTGTGAAATTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTaa                                                                                                                    <  2:116531/141‑1 (MQ=255)
                       gATGTTATCCATGTTTGTACGCCAAATGACTCGCATTGTGAAATTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTaaa                                                                                                                   <  1:505931/141‑1 (MQ=255)
                       gATGTTATCCATGTTTGCACCCCAAATGACCCGCATTGTGAAATTACTATAAGAGGTTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTaaa                                                                                                                   <  1:131016/141‑1 (MQ=255)
                            tATCCATGTTTGTAAGCCAAATGACTCGCTTTGTGAAATTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAAc                                                                                                              >  1:687592/1‑141 (MQ=255)
                             aTCCATGTTTGTACGCCAAATGACTCGCATTGTGAAATTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATCCAGCTAAATCAACa                                                                                                             >  1:661485/1‑141 (MQ=255)
                             aTCCATGTTTGTACGCCAAATGACTCGCATTGTGAAATTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATCCAACTAAAacaaca                                                                                                             >  1:100984/1‑141 (MQ=255)
                               ccATGTTTGGACGCCAAATGACTCGCATTGTGAAATTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAgg                                                                                                           <  2:560039/141‑1 (MQ=255)
                               acATGTTTAGACGCAAAAAGACGAGAAATATGAAATTACTGAAGAAGGCTTGCATGCACGTAAACATGTGATGTGAGAAAAATAAATGGCTAAAACGAAAAAAAAAAAAAAAAAAATGATAGATACAGCTAAATCAACAgg                                                                                                           <  2:428217/140‑1 (MQ=255)
                               acATGTTCATACGCCAAATGAATAGAATAATTAAATTACGGAAGCAGGATTGCATGCCGGTAAACATGTGATGTGTGAAACTCCAATGTCTTAACCGAAAAAAAAAAAAAAAAAAATGATAGATACAGCTAAATCAACAgg                                                                                                           <  2:44050/140‑1 (MQ=255)
                                  tggttGTACCCCTAATGACTCGCAGTGTGAAATTACAGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGTCTAAAACGACAGCAGAATCTCAAAAAATGATAGATACAGCTAAATCAACAGGTaa                                                                                                        <  2:398719/138‑1 (MQ=255)
                                  tGTTTGTACGCCAAATGACTCGCATTGTGAAATTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTaa                                                                                                        <  1:77722/141‑1 (MQ=255)
                                  aatttGTACGACAAATGACTAGAAAATTAAAAATACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGAAAAAAAAAAAAAAATGATAGATACAGCTAAATCAACAGGTaa                                                                                                        <  2:485046/139‑1 (MQ=255)
                                    gttGTACGCCAAATGACTCGCATTGTGAAATTACTGTAGCAGGGTTGCATGCTCGTAAACATAAAAAAAATAAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTaaaa                                                                                                      <  2:381883/140‑1 (MQ=255)
                                       gTACGCCAAATGACTCGCATTGTGAAATTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAAt                                                                                                   <  1:563418/141‑1 (MQ=255)
                                         gcGCCAAATGGCTCGCAGTGTGAAATTACGGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACAACTGGCTAAAAAGACAGTAGTAGCTCAAAAAATTATAGATACAGCTAAATCAACAGGTAAAAAATTa                                                                                                 <  1:583115/140‑1 (MQ=255)
                                         aCGCCAAATGACTCGCATTGTGAAATTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTa                                                                                                 <  1:572769/141‑1 (MQ=255)
                                         aCGCCAAATGACTCGCATTGTGAAATTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTa                                                                                                 <  1:253042/141‑1 (MQ=255)
                                         aCGCCAAAGGACTCGCATTGTGAAATTACTGTAGCAGGGTTGGATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTa                                                                                                 <  1:247737/141‑1 (MQ=255)
                                          cGCCAAATGACTCGCATTGTGAAATTACTGTAACAGGGTTGAATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAAAGACAGAAGAAGCTCAAAAAAGGATAGATAAAGCACAATAAAAAAGTAGAAAACTaa                                                                                                >  2:265601/1‑141 (MQ=255)
                                                aTGACTCGCATTGTGAAATTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCACCAGGTAAAAAATTAACAATAg                                                                                          >  1:55212/1‑141 (MQ=255)
                                                  ggccctcatgtggAAATTACGGTAGCAGGGTTGCATGCTGGTAAACATGTCATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGGt                                                                                        <  2:474012/129‑1 (MQ=255)
                                                    cTCGCATTGTGAAATAACCGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATAAACAATAGGTTa                                                                                      >  2:251874/1‑141 (MQ=255)
                                                         tttGTGAAATTACTGCAGCAGAGTTGTATGCTGGAAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAACCTCAAAACATGATAGTTACAGCTAAATCAACAGGTAAAATATTAACAATAGGTTATCaaa                                                                                 <  1:423339/140‑1 (MQ=255)
                                                           tatgAAATTCCTGTAGCAGGGTTGAATGCTGGTAAACAAGTGAGGGGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGGTTATCaaaac                                                                               >  2:62721/3‑140 (MQ=255)
                                                              gacaTTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGCTGTGTGAAAAACAAATGGCTACAACGACAGCAGAAGCTAAAAAAATGATAGATAAAGCTAAATCAACAGGTAAAAAATGAAACATAGGTTATCAAAAtaga                                                                            >  2:586435/4‑138 (MQ=255)
                                                                taTTACTGAAGCAGGGTTGAATGCTGGTACAATTGGGTTACGTGCAATACCAATGACTCGAGCGAGAATATAAGCGCGGATAACGATAGGGACGGCAAAGTGAAGAGGTAAAAAATTAACAATAGGTTATCAAAATCGttt                                                                          <  2:576773/140‑1 (MQ=255)
                                                                aaTTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGGTTATCAAAATCGttt                                                                          <  1:141760/141‑1 (MQ=255)
                                                                   tACTGTAGCAGGGTTGCGCGCGTATTATGATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATATAGCTAAATCAACAGGTAAAAAATTAACAATGGGTTATCAAAATCGTTTCCg                                                                       <  2:173575/141‑1 (MQ=255)
                                                                      tGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCACCAGGTAAAAAATTAACAATAGGTTATCAAAATCGTTTCCGGGc                                                                    >  1:153969/1‑141 (MQ=255)
                                                                        tAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCGAAATCAACAGGTAAAAAATTAACAATAGGTTATAAAAATCGTTTCCGAGAAg                                                                  >  2:100344/1‑141 (MQ=255)
                                                                           cAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGGTTATCAAAATCGTTTCCGAGCAGATa                                                               >  1:605940/1‑141 (MQ=255)
                                                                               gTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGGTTATCAAAATCGTTTCGCAGCAGATAGTCa                                                           >  1:53820/1‑141 (MQ=255)
                                                                                 tgaatgcctctAAACATGTGATGTGAAAAAAACCAAAGGCTAAAACGACAGAAAAAAAAAAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGGTTATCAAAATCGTTTCCGAGCAGATAGTCAAt                                                         <  2:694108/131‑1 (MQ=255)
                                                                                 tgaaTGCTGGTAAACATGTGATGTGTGACAAACCAATGGCTAAAACGACAGCAGAAGATCAAACAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGGTTATCAAAAGCGTTTCCGAGCAGATAGTCAAt                                                         <  2:277399/138‑1 (MQ=255)
                                                                                  gCATGCTGGTAAAAATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGAGACAGCTAAATCAACAGGGAAAAAATTAACAATAGGTTATCAAAATCGTTTCAGAGCCGATAGTCAAtt                                                        >  2:425367/1‑141 (MQ=255)
                                                                                   cATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGGTTAACAAAATCGTTTCCGGGCAGAAAGGCAAttt                                                       >  1:486154/1‑141 (MQ=255)
                                                                                             aaCATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGGTTATCAAAATCGTTTCCGAGCAGATAGTCAATTTTTACATCAag                                             >  2:657793/1‑141 (MQ=255)
                                                                                             aaCATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGGTTATCAAAATCGTTTCCGAGCAGATAGTCAATTTTTACATCAag                                             <  2:656893/141‑1 (MQ=255)
                                                                                                   tgatgtgTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGGTTATCAAAATCGTTTCCGAGCAGATAGTCAATTTTTACATCAAGCAgcgc                                       >  2:333524/1‑141 (MQ=255)
                                                                                                               aacgaaTGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGGTTATCAAAAACGTTTCCGAGCAGATAGTCAATTTTTACATCAAGCAGCGCAACGTGGCGACt                           >  1:499471/5‑141 (MQ=255)
                                                                                                                  cAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGGTTATCAAAATCGTTTCCGAGCAGATAGTCAATTTTTACATCAAGCAGCGCAACGTGGCGACTTAg                        >  1:627135/1‑141 (MQ=255)
                                                                                                                      ggCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGGTTATCAAAATCGTTTCCGAGCAGATAGTCAATTTTTACATCAAGCAGCGCAACGTGGCGACTTAGgaga                    >  1:592090/1‑141 (MQ=255)
                                                                                                                              cGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGGTTATCAAAATCGTTTCCGAGCAGATAGTCAATTTTTACATCAAGCAGCGCAACGTGGCGACTTAGGAGCCATTTACt            >  1:70898/1‑141 (MQ=255)
                                                                                                                               gACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGGTTATCAAAATCGTTTCCGAGCAGATAGTCAATTTTTACATCAAGCAGCGCAACGTGGCGACTTAGGAGACATTTACtt           >  1:97318/1‑141 (MQ=255)
                                                                                                                               gACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGGTTATCAAAATCGTTTCCGAGCAGATAGTCAATTTTTACATCAAGCAGCGCAACGTGGCGACTTAGGAGACATTTACtt           >  1:481271/1‑141 (MQ=255)
                                                                                                                                 catcacaaGATCAAAAAATAAGAAATACAGCTAAATCAACCGGTAACAACTTAACCATAGGTTATCACAATCGTTTCCGAGCAGATAGTCAATTTTTACATCAAGCAGCGCAACGTTGCGACTTAGGAGACATTTACTTCg         <  2:261193/135‑1 (MQ=255)
                                                                                                                                 cagcagAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGGTTATCAAAATCGTTTCCGAGCAGATAGTCAATTTTTACATCAAGCAGCGCAACGTGGCGACTTAGGAGACATTTACTTCg         >  1:419894/1‑141 (MQ=255)
                                                                                                                                    cagAAGCTCAAAAAATGATAGATACAGCTAACTCAACAGGTAAAAAATTAACAATAGGTTATCAAAATCGTTTCCGAGCAGATAGTCAATTTTTACATCAAGCAGCGCAACGTGGCGACTTAGGAGAAATTTACTTTGGaa      >  2:616567/1‑141 (MQ=255)
                                                                                                                                    cagAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGGTTATCAAAATCGTTTCCGAGCAGATAGTCAATTTTTACATCAAGCAGCGCAACGTGGCGACTTAGGAGACATTTACTTCGGaa      >  1:43766/1‑141 (MQ=255)
                                                                                                                                     agAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGGTTATCAAAATCGTTTCCGAGCAGATAGTCAATTTTTACATCAAGCAGCGCAACGTGGCGACTTAGGAGACATTTCCTTCGGaaa     >  1:638175/1‑141 (MQ=255)
                                                                                                                                     agAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGGTTATCAAAATCGTTTCCGAGCAGATAGTCAATTTTTACATCAAGCAGCGCAACGTGGCGACTTAGGAGACATTTACTTCGGaaa     >  1:439523/1‑141 (MQ=255)
                                                                                                                                     agAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGCTTATCAAAATCGTTTCCGAGCAGATAGTCAATTTTTACATCAAGCAGCGCAACGTGGCGACTTAGGAGACATTTACTTCGGaaa     <  2:403466/141‑1 (MQ=255)
                                                                                                                                        aGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGGTTATCAAAATCGTTTCCGAGCAGATAGTCAATTTTTACATTAATCAGCGCAACGTGGCGACTTAGGAGACATTTACTTCGGAAAGGc  >  2:138014/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                        
CATTGTTAAAAGATGACACGATTGATGTTATCCATGTTTGTACGCCAAATGACTCGCATTGTGAAATTACTGTAGCAGGGTTGCATGCTGGTAAACATGTGATGTGTGAAAAACCAATGGCTAAAACGACAGCAGAAGCTCAAAAAATGATAGATACAGCTAAATCAACAGGTAAAAAATTAACAATAGGTTATCAAAATCGTTTCCGAGCAGATAGTCAATTTTTACATCAAGCAGCGCAACGTGGCGACTTAGGAGACATTTACTTCGGAAAGGC  >  CP000730/250804‑251080

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: ‑