Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 290,156 A→T 5.5% intergenic (+47/‑180) USA300HOU_0257 → / → gatC1 PTS family galactitol (gat) porter component IIB/PTS family galactitol (gat) porter component IIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730290,1560AT 5.5% 232.8 / 5.4 73intergenic (+47/‑180)USA300HOU_0257/gatC1PTS family galactitol (gat) porter component IIB/PTS family galactitol (gat) porter component IIC
Reads supporting (aligned to +/- strand):  ref base A (36/33);  new base T (2/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

AAATTCCCAGTTATCAACGCACGTAATTTCTTAACTGGTATTGGTATTGAAGAAACAAAACAACAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGAAATGACTGTAACGTCAAATGTAATCACAATCGGAAA  >  CP000730/290020‑290296
                                                                                                                                        |                                                                                                                                            
aaaTTCCCAGTTATCAACGCACGTAATTTCTTAACTGGTATTGGTATTGAAGAAACAAAACAACAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGc                                                                                                                                          <  2:634198/141‑1 (MQ=255)
  aTTCCCAGTTATCAACGCACGTAATTTCTTAACTGGTATTGGTATTGAAGAAACAAAACAACAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCat                                                                                                                                        >  2:481066/1‑141 (MQ=255)
     cccAGTTATCAACGCACGTAATTTCTTAACTGGTATTGGTATTGAAGAAACAAAACAACAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCAtata                                                                                                                                     <  1:485441/141‑1 (MQ=255)
      ccAGTTATCAACGCACGTAATTTCTTAACTGGTATTGGTATTGAAGAAACAAAACAACAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATaa                                                                                                                                    >  2:622851/1‑141 (MQ=255)
       cAGTTATCAACGCACGTAATTTCTTAACTGGTATTGGTATTGAAGAAACAAAACAACAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATaaa                                                                                                                                   >  1:301556/1‑141 (MQ=255)
       cAGTTATCAACGCACGTAATTTCTTAACTGGTATTGGTATTGAAGAAACAAAACAACAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATaaa                                                                                                                                   >  2:623738/1‑141 (MQ=255)
       cAGTTATCAACGCACGTAATTTCTTAACTGGTATTGGTATTGAAGAAACAAAACAACAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATaaa                                                                                                                                   <  1:10653/141‑1 (MQ=255)
         gTTATCAACGCACGTAATTTCTTAACTGGTATTGGTATTGAAGAAACAAAACAACAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATCCAAATTGCATATAAATg                                                                                                                                 >  2:37758/1‑141 (MQ=255)
           tATCAACGCACGTAATTTCTTAACTGGTATTGGTATTGAAGAAACAAAACAACAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAATTTGCAtatattaata                                                                                                                               <  2:103152/141‑7 (MQ=255)
           tATCAACGCACGTAATTTCTTAACTGGTATTGGTATTGAAGAAACAAAACAACAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAATTTGCAtatattaata                                                                                                                               <  2:103038/141‑7 (MQ=255)
               aaCGCACGTAATTTCTTAACTGGTATTGGTATTGAAGAAACAAAACAACAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATc                                                                                                                           >  1:98584/1‑141 (MQ=255)
                       tgaTTTCTTAACTGGTATTGGTATTGAAGAAACAAAACAACAAATCTTAAGAGAGTTACAAAAATAACTGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACaaa                                                                                                                   <  2:609189/139‑1 (MQ=255)
                        aacttCTTAACTGGTATTGGTATTGAAGAAACAAAACAACAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAAt                                                                                                                  <  1:503280/138‑1 (MQ=255)
                          tttCTTAACTGGTATTGGTATTGAAGAAACAAAACAACAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGa                                                                                                                <  1:439729/141‑1 (MQ=255)
                            tCTTAACTGGTATTGGTATTGAAGAAACAAAACAACAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGAtt                                                                                                              <  1:79229/141‑1 (MQ=255)
                                 aCTGGTATTGGTATTGAAGAAACAAAACAACAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTtga                                                                                                         >  1:369464/1‑141 (MQ=255)
                                      tattggtattgAAGAAACAAAACAACAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACa                                                                                                    <  1:414860/141‑1 (MQ=255)
                                       attggtattgAAGAAACAAAACAACAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACAt                                                                                                   >  1:76055/1‑141 (MQ=255)
                                              ttgAAGAAACAAAACAACAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACg                                                                                            >  1:578471/1‑141 (MQ=255)
                                              ttgAAGAAACAAAACAACAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACg                                                                                            >  1:409636/1‑141 (MQ=255)
                                                    aaacaaaacaACAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGa                                                                                      >  2:28978/1‑141 (MQ=255)
                                                    aaacaaaacaACAAATCTTAAAAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGa                                                                                      <  2:357656/141‑1 (MQ=255)
                                                       caaaacaaCAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATaa                                                                                                                                    >  1:397749/1‑92 (MQ=255)
                                                       caaaacaaCAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATaa                                                                                                                                    <  2:397749/92‑1 (MQ=255)
                                                       caaaacaaCAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATTACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTg                                                                                   <  1:589485/141‑1 (MQ=255)
                                                       caaaacaaCAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTg                                                                                   >  1:580771/1‑141 (MQ=255)
                                                        aaaacaaCAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGt                                                                                  <  2:692518/141‑1 (MQ=255)
                                                          aacaacAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTaa                                                                                >  2:438959/1‑141 (MQ=255)
                                                             aacaaATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTc                                                                             <  1:146824/141‑1 (MQ=255)
                                                              acaaATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCa                                                                            <  2:707809/141‑1 (MQ=255)
                                                              acaaATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCa                                                                            <  2:45854/141‑1 (MQ=255)
                                                               caaATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTg                                                                                       <  2:209884/129‑1 (MQ=255)
                                                               caaATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTg                                                                                       >  1:209884/1‑129 (MQ=255)
                                                               caaATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCaa                                                                           >  1:352781/1‑141 (MQ=255)
                                                                 aaTCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCaaaa                                                                         >  1:504846/1‑141 (MQ=255)
                                                                       aaCAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCAAGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTg                                                                   >  1:404361/1‑141 (MQ=255)
                                                                        aCAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGc                                                                  >  1:384619/1‑141 (MQ=255)
                                                                        aCAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGc                                                                  >  1:114682/1‑141 (MQ=255)
                                                                               tACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGaa                                                           >  1:489924/1‑141 (MQ=255)
                                                                                aCAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTtga                                                                                                         <  2:74046/94‑1 (MQ=255)
                                                                                aCAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTtga                                                                                                         >  1:74046/1‑94 (MQ=255)
                                                                                      aTAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCa                                                    <  2:359639/141‑1 (MQ=255)
                                                                                         aCGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAgg                                                 >  1:495705/1‑141 (MQ=255)
                                                                                          cGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAggg                                                >  2:506805/1‑141 (MQ=255)
                                                                                           ggAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGt                                               >  1:380971/1‑141 (MQ=255)
                                                                                           ggAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGt                                               <  2:253412/141‑1 (MQ=255)
                                                                                             aaTTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTg                                             >  2:157850/1‑141 (MQ=255)
                                                                                                tGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATATATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCg                                          <  1:180808/141‑1 (MQ=255)
                                                                                                 ggtttGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGATAAATGATTGTTGATGACATAGAAACGCAGTGTCTGTAATTCAAAAGGACTGCAGTAGAATTTAGCAAGGGTTGTCGt                                         <  2:473575/137‑1 (MQ=255)
                                                                                                    atGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGa                                      <  1:497253/141‑1 (MQ=255)
                                                                                                    atGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGa                                      <  2:532104/141‑1 (MQ=255)
                                                                                                    atGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGa                                      <  1:3238/141‑1 (MQ=255)
                                                                                                             gggTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGAAATGACTGt                             <  2:547059/141‑1 (MQ=255)
                                                                                                             gggTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGAAATGACTGt                             >  2:265494/1‑141 (MQ=255)
                                                                                                             gggTCGAGAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCTAGGGATGTCGTTGAAAGGACTGt                             >  2:314115/1‑141 (MQ=255)
                                                                                                               gTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAAAGAAGAGACTAAAAATAAAAAGGACTACAGTGGAATTTAGAAAGGGGTGTCGTTGAAATGACTGTaa                           >  1:418706/1‑141 (MQ=255)
                                                                                                                 cGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGAAATGACTGTAACg                         <  2:214579/141‑1 (MQ=255)
                                                                                                                 cGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGAAATGACTGTAACg                         <  2:339256/141‑1 (MQ=255)
                                                                                                                   aTAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGAAATGACTGTAACGTc                       <  1:79402/141‑1 (MQ=255)
                                                                                                                    tAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGAAATGACTGTAACGTCa                      >  1:31985/1‑141 (MQ=255)
                                                                                                                     aaCATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGAAATGACTGTAACGTCaa                     >  2:15044/1‑141 (MQ=255)
                                                                                                                      aCATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGAAATGACTGTAACGTCaaa                    >  1:610568/1‑141 (MQ=255)
                                                                                                                      aCATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGAAATGACTGTAACGTCaaa                    >  1:363236/1‑141 (MQ=255)
                                                                                                                         tatTAAGATATGCAAATTGCATATAAATGACTATCACGATTAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAATGGACTGCAGTCGAAATTAGCAAGGGTTGTCGTCGCAAGGACTGAAACGTCAAATGt                 <  2:638776/141‑1 (MQ=255)
                                                                                                                           ttAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGAAATGACTGTAACGTCAAATGTaa               >  1:343482/1‑141 (MQ=255)
                                                                                                                           ttAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGAAATGACTGTAACGTCAAATGTaa               >  1:416576/1‑141 (MQ=255)
                                                                                                                            tAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGAAATGACTGTAACGTCAAATGTAAt              <  2:511861/141‑1 (MQ=255)
                                                                                                                             aaTATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGAAATGACTGTAACGTCAAATGTAATc             <  1:394315/141‑1 (MQ=255)
                                                                                                                             aaTATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGAAATGACTGTAACGTCAAATGTAATc             <  2:465807/141‑1 (MQ=255)
                                                                                                                                 tatGCAATTTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGAAATGACTGTAACGTCAAATGTAATCACaa         >  2:609376/1‑141 (MQ=255)
                                                                                                                                   tGCAATTTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGAAATGACTGTAACGTCAAATGTAATCACAATc       >  1:390959/1‑141 (MQ=255)
                                                                                                                                     cAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTc                                                              <  1:294642/84‑1 (MQ=255)
                                                                                                                                     cAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTc                                                              >  2:294642/1‑84 (MQ=255)
                                                                                                                                      aaaTTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGAAATGACTGTAACGTCAAATGTAATCACAATCGGa    <  2:426045/141‑1 (MQ=255)
                                                                                                                                      aaaTTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGAAATGACTGTAACGTCAAATGTAATCACAATCGGa    >  1:269860/1‑141 (MQ=255)
                                                                                                                                        aTTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGAAATGACTGTAACGTCAAATGTAATCACAATCGGaaa  >  2:64559/1‑141 (MQ=255)
                                                                                                                                        aTTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGAAATGACTGTAACGTCAAATGTAATCACAATCGGaaa  <  1:481066/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                            
AAATTCCCAGTTATCAACGCACGTAATTTCTTAACTGGTATTGGTATTGAAGAAACAAAACAACAAATCTTAACAGAGTTACAAAAATAACGGAATTGATATGTAACGTGGGTCGATAACATATTAATATATGCAAATTGCATATAAATGACTATCACGACAAATGATTGTTGATGACATAGAAACGCAGTGACTGTAATTCAAAAGGACTGCAGTCGAATTTAGCAAGGGTTGTCGTTGAAATGACTGTAACGTCAAATGTAATCACAATCGGAAA  >  CP000730/290020‑290296

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