Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,929,842 Δ1 bp 100% intergenic (‑925/‑243) USA300HOU_1794 ← / → USA300HOU_1796 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,929,8420T.57.7% 33.1 / 138.3 71intergenic (‑925/‑243)USA300HOU_1794/USA300HOU_1796hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (16/13);  new base . (22/19);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

TAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGC  >  CP000730/1929708‑1929977
                                                                                                                                      |                                                                                                                                            
tAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGc                                                                                                                                        >  2:233838/1‑141 (MQ=255)
   aaTTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCaga                                                                                                                                     <  2:55197/141‑1 (MQ=255)
     tttGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCagaga                                                                                                                                   <  1:297053/141‑1 (MQ=255)
     tttGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCagaga                                                                                                                                   <  2:372714/141‑1 (MQ=255)
        gTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAAGCGTCAATTTAGACGCAGGgaggg                                                                                                                                >  1:44754/1‑140 (MQ=255)
         tCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAgaggag                                                                                                                               >  2:289397/1‑141 (MQ=255)
                           aTTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCtt                                                                                                        <  2:201988/141‑1 (MQ=255)
                              cAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTAtt                                                                                                     >  2:93331/1‑141 (MQ=255)
                                aTCTATCAGCTAGAATCACATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACCCAGAGAGGAGGTGTATAAGGGGATGCTTAtttt                                                                                                   >  1:177470/1‑141 (MQ=255)
                                     tCAGCTAGAATCATATGAATTCTGATTATCAAGTGACTGCGCTTGTATATTTGCCGCTGATATAAAAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGAAGAGAGGAGGTGTATAAGGTGAGGCTTATTTTAGTTc                                                                                              >  2:85934/1‑141 (MQ=255)
                                     tCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGGCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTc                                                                                              >  2:288682/1‑141 (MQ=255)
                                            gAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTTAcatcat                                                                                       >  2:323653/1‑141 (MQ=255)
                                                          aTCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGtcatca                                                                         >  2:51878/1‑141 (MQ=255)
                                                           tcggaaggagggggcTGTGCTTGGATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTGAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAg                                                                        <  2:382244/127‑1 (MQ=255)
                                                           tCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAg                                                                        <  2:330569/141‑1 (MQ=255)
                                                              atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                    <  2:61360/127‑1 (MQ=255)
                                                               tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTAAAATCATAGACCCAGTCATCAGTGGCt                                                                   >  1:365510/14‑141 (MQ=255)
                                                                cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                  <  1:261897/129‑1 (MQ=255)
                                                                  gaagtagatatgCCTGTATATTCACAGATGCTACCATAACAGCGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                <  1:9700/131‑1 (MQ=255)
                                                                    agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                                              >  1:300757/9‑140 (MQ=255)
                                                                     gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                             <  1:341053/134‑3 (MQ=255)
                                                                     gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCGATCTCTCAGGCGTCAATTAACGAAGAGAGGAGCTTCATTACTTTGCTAATCATTTTCGTTCCCCTCATCGCACCAGTCATCAGTGGCTGTGctg                                                             >  2:76933/8‑139 (MQ=255)
                                                                     gTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGcc                                                              >  2:211986/1‑141 (MQ=255)
                                                                      tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                            <  1:264032/135‑4 (MQ=255)
                                                                       agatatgCTTGTATATTCACAGATGGTATAATAATAGAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTAGTTCACATAATAGCACCAGTCATCAGTGCATGTGctgtt                                                           >  2:248424/6‑137 (MQ=255)
                                                                        gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                          <  1:307878/137‑6 (MQ=255)
                                                                        gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                          <  1:124587/137‑6 (MQ=255)
                                                                          tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgcg                                                        <  2:73250/139‑8 (MQ=255)
                                                                          tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                                                        >  2:339333/3‑141 (MQ=255)
                                                                           atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                       >  1:23854/2‑141 (MQ=255)
                                                                                tGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat                                                  <  1:303477/141‑6 (MQ=255)
                                                                                  tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCGTCAGTGGCTGTGCTGTTGCGtattatac                                                <  2:67438/141‑8 (MQ=255)
                                                                                    tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACtt                                              >  2:79705/9‑141 (MQ=255)
                                                                                     tttcggcGATGATATATTAGTAGAGGCGCCTATCTGTCAGGCGTCAACTTAGACGCAGAGAGGATGTGTATAAGGTGATGCTTATTTTCGTTCACATCGTAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCtt                                              <  1:85934/135‑1 (MQ=255)
                                                                                     atTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCTTTTTTTTCtt                                              <  1:288682/141‑1 (MQ=255)
                                                                                       ttGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTAt                                            <  2:18332/141‑1 (MQ=255)
                                                                                       tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                           >  1:37368/6‑141 (MQ=255)
                                                                                          cagacggtataATTATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCAATTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                        <  1:87477/134‑1 (MQ=255)
                                                                                           cGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGCGAGGAGGTGTATAAGGCGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc                                        <  1:260833/141‑1 (MQ=255)
                                                                                             agcatatAATAGGAGAGTCGCCTATCTCTCAGGCGTCAATTTAGCCGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTa                                      <  2:215536/138‑1 (MQ=255)
                                                                                              tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATCCTTATTGGCTTAg                                    >  1:220004/4‑141 (MQ=255)
                                                                                              tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                    <  1:141444/138‑4 (MQ=255)
                                                                                              tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTACTCATTTTCGTTCACGTCATAGCACCTGTCATCAGTGGCTGTGCTGCTGCGTATTATACTTATTGGCTTAg                                    <  2:275815/138‑1 (MQ=255)
                                                                                              tgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCACCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTaa                                     >  2:238635/1‑141 (MQ=255)
                                                                                               ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAgg                                   >  2:291720/3‑140 (MQ=255)
                                                                                               ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                   >  2:319625/3‑141 (MQ=255)
                                                                                                  ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                >  1:313259/1‑139 (MQ=255)
                                                                                                   taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAAGATCGGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                               <  1:30533/141‑4 (MQ=255)
                                                                                                   gaataatAGGGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGCGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTCTTGGCTTAGTaaac                               <  2:272455/140‑4 (MQ=255)
                                                                                                     ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTATTAAACGc                             >  2:159873/5‑141 (MQ=255)
                                                                                                     ataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACg                              >  1:335014/1‑141 (MQ=255)
                                                                                                     ataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCCCTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACg                              >  1:290362/1‑141 (MQ=255)
                                                                                                      taatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                            <  1:297687/138‑1 (MQ=255)
                                                                                                      caatatagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                            <  2:196202/135‑1 (MQ=255)
                                                                                                        atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GAGGAAGAGAGGAGGTGCATTACTGTGCTAATCATTTTCGTTCGCATCATAGCACCAGTCATCAGTCGCTGTGCTGTGGCGTCTTATACTTATTGGCTTAGTCAACGCaac                          >  2:301096/2‑140 (MQ=255)
                                                                                                             gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaa                      <  1:24325/141‑1 (MQ=255)
                                                                                                             gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacgaata                     >  1:12308/1‑135 (MQ=255)
                                                                                                             gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                     >  2:286085/1‑135 (MQ=255)
                                                                                                              tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                    >  2:314126/1‑134 (MQ=255)
                                                                                                                gcctttccgtcAGGCTTCAATTTAAGCGCAGAGGGGAGGTGTATAAGGTGATGCTTATGTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTGTTTCTTATTGGCTAAGTAGACGCAATACAAAAt                   <  2:77131/132‑1 (MQ=255)
                                                                                                                gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAATCATCAGTGGCTTTGCTGTTGCGTATCATACTTATTGGCTTAGTAAACGCaacaaatgatt                  >  2:102184/1‑132 (MQ=255)
                                                                                                                gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                  >  1:21717/1‑132 (MQ=255)
                                                                                                                gCCTATCTCGCAGGCGTCAAGT‑‑GGCGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTGCGTTCACATCATAGCACCAGTCATCCGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                  <  1:76933/141‑10 (MQ=255)
                                                                                                                  cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                <  1:312343/141‑12 (MQ=255)
                                                                                                                  cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTCAATGGCTTAGTAAACGCaacaaataattaa                >  1:187558/1‑130 (MQ=255)
                                                                                                                   tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTAATTGGCTTAGTAAACGCaacaaataattaaa               >  2:105831/1‑129 (MQ=255)
                                                                                                                   tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGCAGGTGCATTACTTGGCTAATCATTTTCACTCACAACATAGCACAAGTCATGAGTGGCTGTGCTGATGCGTATTATACTTATTGGCTTAGTAAACGCaacaactaattaat               >  2:149117/1‑129 (MQ=255)
                                                                                                                     tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTACTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag             >  2:295300/1‑127 (MQ=255)
                                                                                                                         tcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCCTTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATaacatacct          >  2:278132/1‑133 (MQ=255)
                                                                                                                         tcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACat          <  1:93331/141‑1 (MQ=255)
                                                                                                                                 cAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGc  >  2:352045/1‑141 (MQ=255)
                                                                                                                                      |                                                                                                                                            
TAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTATAAGGTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGC  >  CP000730/1929708‑1929977

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