Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 824,661 Δ1 bp 4.6% coding (231/1071 nt) USA300HOU_0773 ← diguanylate cyclase GGDEF domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730824,6600C. 4.6% 285.0 / 9.3 87coding (232/1071 nt)USA300HOU_0773diguanylate cyclase GGDEF domain protein
Reads supporting (aligned to +/- strand):  ref base C (45/38);  new base . (2/2);  total (47/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.13e-01

CATTTTGTTTTAAAAATGGTCCAATAGCACTTGTTACTGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACAGATAATAATAATGATACGATTGTCATTAAAACAGTGACATAAGCCTTAGAA  >  CP000730/824525‑824796
                                                                                                                                       |                                                                                                                                        
catTTTGTTTTAAAAATGGTCCAATAGCACTTGTTACTGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCatat                                                                                                                                     >  2:948567/1‑141 (MQ=255)
 atTTTGTTTTAAAAATGGTCCAATAGCACTTGTTACTGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATAtt                                                                                                                                    >  1:283239/1‑141 (MQ=255)
   tttGTTTTAAAAATGGTCCAATAGCACTTGTTACTGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCGTATTTg                                                                                                                                  >  1:541247/1‑141 (MQ=255)
   tttGTTTTAAAAATGGTCCAATAGCACTTGTTACTGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTg                                                                                                                                  >  1:976299/1‑141 (MQ=255)
   tttGTTTTAAAAATGGTCCAATAGCACTTGTTACTGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTg                                                                                                                                  >  2:713851/1‑141 (MQ=255)
    ttGTTTTAAAAATGGTCCAATAGCACTTGTTACTGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGt                                                                                                                                 >  1:904708/1‑141 (MQ=255)
    ttGTTTAAAAAATGGTCCAATAGCACTTGTTACTGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTCGCAGATAATGTATAGACCATATTTGt                                                                                                                                 <  2:684300/141‑1 (MQ=255)
           aaaaaTGGGCCAATACCACTTGTTAGTGCAGCAAAAACAATCAACGTAAAGCCTGATATAATCGAGTTGTTAAATACTACTACTTTAACAATTGCTACAAGTACTGTGTCAGATAAGGTTTAGACCGTATTTGGAAATTTa                                                                                                                          <  2:278177/141‑1 (MQ=255)
             agaTGGTCCAATAGCACTTGTTACTGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTAcc                                                                                                                        <  2:195097/139‑1 (MQ=255)
             aaaTGGTCCAATAGCACTTGTTACTGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTAcc                                                                                                                        <  1:501829/141‑1 (MQ=255)
             aaaTGGTCCAATAGCAATTGTTACTGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACCATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTAcc                                                                                                                        <  2:236610/141‑1 (MQ=255)
               aTGGTCCAATAGCACTTGTTACTGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTa                                                                                                                      <  2:253665/141‑1 (MQ=255)
                tGGTCCAATAGCACTTGTTACTGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTaa                                                                                                                     >  1:512668/1‑141 (MQ=255)
                 ggTCCAATAGCACTTGTTACTGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTaaa                                                                                                                    <  1:199792/141‑1 (MQ=255)
                         aGCACTTGTTACTGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATaa                                                                                                            <  2:139766/141‑1 (MQ=255)
                           cACTTGTTACTGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAg                                                                                                          >  2:980973/1‑141 (MQ=255)
                            aCTTGTTACTGCAGCACTAACAATCAACGTAACACCGTACATAATCGAGTCGTTAAAAACTCCAATTCCAACAATTGCTATAAATACTGTGGCCGATCATGTCTAGACCATATCTGTAAATATCACTaaaagaagtaaaaa                                                                                                         >  2:235077/1‑131 (MQ=255)
                            aCTTGTTACTGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAgg                                                                                                         >  2:330957/1‑141 (MQ=255)
                              ttGTTACTGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGaa                                                                                                       <  2:828367/141‑1 (MQ=255)
                               tGTTACTGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAAc                                                                                                      >  2:341420/1‑141 (MQ=255)
                               tGTTACTGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTCTAAATCTACCTAAAAACCATAAAGGAAc                                                                                                      >  1:77992/1‑141 (MQ=255)
                                gTTACTGCAGCAATTACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACg                                                                                                     <  1:756673/141‑1 (MQ=255)
                                gTTACTGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACg                                                                                                     <  1:756551/141‑1 (MQ=255)
                                   aCTGCAGCAATAACCATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCCTATTTGTAAATCTACCTAAAAACAATAAAGAAACGAAt                                                                                                  >  2:125115/1‑141 (MQ=255)
                                   aCTGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGaac                                                                                                  >  2:225868/1‑140 (MQ=255)
                                   aCTGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATATACCTAAAAACAATAAAGGAACGaag                                                                                                  >  2:751229/1‑140 (MQ=255)
                                   aCTGCAGCAATAACAATAAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAAt                                                                                                  >  2:925333/1‑141 (MQ=255)
                                     tGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGtt                                                                                               >  1:285382/1‑141 (MQ=255)
                                     tGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGtt                                                                                               >  1:220782/1‑141 (MQ=255)
                                         gcaATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACCATAAAGGAACGAATGTTAAg                                                                                            <  1:668774/141‑1 (MQ=255)
                                         gcaATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAg                                                                                            <  1:887642/141‑1 (MQ=255)
                                         acaATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAg                                                                                            <  1:243736/140‑1 (MQ=255)
                                              aaCAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAAt                                                                                       <  2:520588/141‑1 (MQ=255)
                                              aaCAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAAt                                                                                       <  2:350425/141‑1 (MQ=255)
                                               aCAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAAtt                                                                                      >  2:845141/1‑141 (MQ=255)
                                               aCAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCAAAAAACAATAAAGGAACGAATGTTAAGTGAAtt                                                                                      >  2:4852/1‑141 (MQ=255)
                                                   taaaCGTAACACCTTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAaat                                                                                  <  2:632836/139‑1 (MQ=255)
                                                   tCAACTTAACACCGTACATAAACAAGTTGTTAAATACTACAATGTCAACAATTGCTACAATTACGGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAaat                                                                                  <  2:214157/141‑1 (MQ=255)
                                                   tCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAaat                                                                                  <  2:300047/141‑1 (MQ=255)
                                                   tCAAAGCAACTCCGTACATAAACGAGTGGTTAACTACTACAATTTCAACAAGTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAaat                                                                                  <  2:249786/141‑1 (MQ=255)
                                                   aaaaCGCAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCACCTTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAaat                                                                                  <  2:551103/139‑1 (MQ=255)
                                                     aaCGTAACCCCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAaataa                                                                                <  1:166764/141‑1 (MQ=255)
                                                     aaCGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTCCAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGACTGATAAGTGAATTAaataa                                                                                <  1:617617/141‑1 (MQ=255)
                                                     aaCGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTACATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAaataa                                                                                <  1:767243/141‑1 (MQ=255)
                                                     aaCGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAaataa                                                                                <  1:664920/141‑1 (MQ=255)
                                                       cGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACACTTGCTACAATTACTGTGGCAGATAATGTATAGACCATAGTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAaataatt                                                                              >  1:1047992/1‑140 (MQ=255)
                                                         tAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCtt                                                                            <  2:424954/141‑1 (MQ=255)
                                                          aaCACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTc                                                                           <  1:980941/141‑1 (MQ=255)
                                                           acacCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATATGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCa                                                                          <  1:676435/141‑1 (MQ=255)
                                                             acCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACg                                                                        <  2:377307/141‑1 (MQ=255)
                                                                   cATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAgg                                                                  >  1:278184/1‑141 (MQ=255)
                                                                   cATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAgg                                                                  >  1:74355/1‑141 (MQ=255)
                                                                   cATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAgg                                                                  >  1:818797/1‑141 (MQ=255)
                                                                      aaTCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGAt                                                               <  1:350426/141‑1 (MQ=255)
                                                                            gttgttAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTa                                                         >  1:443500/1‑141 (MQ=255)
                                                                            gttgttAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTa                                                         >  2:383359/1‑141 (MQ=255)
                                                                            gttgttAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGAAGGGGTa                                                         >  2:592311/1‑141 (MQ=255)
                                                                               gttAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGAc                                                      <  1:905634/141‑1 (MQ=255)
                                                                               gttAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAAAATTAACGGTAAGGGATAGGGTAGaa                                                      >  2:581681/1‑140 (MQ=255)
                                                                                ttAAATACTACCATTTCAACAATTGCTACAATTACGGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACa                                                     <  2:366319/141‑1 (MQ=255)
                                                                                ttAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACa                                                     >  2:788398/1‑141 (MQ=255)
                                                                                ttAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGGATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACa                                                     >  1:439852/1‑141 (MQ=255)
                                                                                ttAAATACTACAATTTCAACAATTGCTACAAGTACTGTGGCAGATAATGTATAGACCTTATTTGTAAATCTACCTAAAAACAATAAAGAAACGAATGTAAAGTGAATGAAATAATCTTGACGATAAggagtagtgaagaaa                                                     >  2:176551/1‑128 (MQ=255)
                                                                                  aaaTACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACAGa                                                   >  2:299659/1‑141 (MQ=255)
                                                                                  aaaTACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACAGa                                                   >  2:750719/1‑141 (MQ=255)
                                                                                  aaaTACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGAGAAGGGATAGGGTAGACAGa                                                   >  1:287436/1‑141 (MQ=255)
                                                                                  aaaTACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATGAAATAATCTTCACGATAAGGGATAGGGTAGACAGa                                                   >  1:1065320/1‑141 (MQ=255)
                                                                                     tactgcAATGTCAACAATTGCTACAAGTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACAGataa                                                <  2:428380/141‑1 (MQ=255)
                                                                                     tactacAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACAGataa                                                <  1:98623/141‑1 (MQ=255)
                                                                                     tactacAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACAGataa                                                <  1:1075566/141‑1 (MQ=255)
                                                                                      actacAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACAGataat                                               >  1:663679/1‑141 (MQ=255)
                                                                                        tacAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACAGataataa                                             >  2:426058/1‑141 (MQ=255)
                                                                                                caacaaTTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACAGATAATAATAATGATa                                     >  1:605566/1‑141 (MQ=255)
                                                                                                  acaaTTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGCATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACAGATACTAAGAATGATACg                                   >  1:76243/1‑141 (MQ=255)
                                                                                                     aTTGCTACAATTACTGTGGCAGATAATGTATAGACATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACAGATAATAATAATGATACGATTg                               <  2:220782/141‑1 (MQ=255)
                                                                                                     aTTGCTACAATTACTGTGGCAGATAATGTATAGACATATTTGTAAAGCTACCTAAAAACAATAAAGGAACGAATGTTAAGGGAATTAAATAATCTTCACGATAAGGGAGAGGGTAGACAGATATTAATAATGATACGATTg                               <  2:285382/141‑1 (MQ=255)
                                                                                                       tGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACAGATAATAATAATGATACGATTGt                              >  2:697805/1‑141 (MQ=255)
                                                                                                       tGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAAACAGATAATAATACTGATACGATTGt                              >  2:938772/1‑141 (MQ=255)
                                                                                                              aTTACTGTGGCAGATAATGTATAGACCATATTAGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACAGATAATAATAATGATACGATTGTAATTaaa                       >  1:867705/1‑141 (MQ=255)
                                                                                                               ttACTGTGGCAGATAATGTATAGACCATATTTGTAAATCAACCGAAAACCAATTAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACAGATAATAATAATGATACGATTGTCATTaaaa                      >  1:530701/1‑141 (MQ=255)
                                                                                                                   tgtgGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACAGATAATAATAATGATACGATTGTCATTAAAACAgg                  >  2:278050/1‑140 (MQ=255)
                                                                                                                           aTAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACAGATAATAATACTGATACGATTGTCATTAAAACAGTGACATAAg          >  1:673076/1‑141 (MQ=255)
                                                                                                                             aaTGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACAGATAATAATAATGATACGATTGTCATTAAAACAGTGACATAAGcc        <  1:980973/141‑1 (MQ=255)
                                                                                                                              aTGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACAGATAATAATAATGATACGATTGTCATTAAAACAGTGACATAAGCCt       >  2:841074/1‑141 (MQ=255)
                                                                                                                              aTGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACAGATAATAATAATGATACGATTGTCATTAAAACAGTGACATAAGCCt       <  1:299659/141‑1 (MQ=255)
                                                                                                                              aTGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACAGATAATAATAATGATACGATTGTCATTAAAACAGTGACATAAGCCt       <  2:530701/141‑1 (MQ=255)
                                                                                                                                   taGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACAGATAATAATAATGATACGATTGTCATTAAAACAGTGACATAAGCCTTAGaa  >  1:667348/1‑141 (MQ=255)
                                                                                                                                       |                                                                                                                                        
CATTTTGTTTTAAAAATGGTCCAATAGCACTTGTTACTGCAGCAATAACAATCAACGTAACACCGTACATAATCGAGTTGTTAAATACTACAATTTCAACAATTGCTACAATTACTGTGGCAGATAATGTATAGACCATATTTGTAAATCTACCTAAAAACAATAAAGGAACGAATGTTAAGTGAATTAAATAATCTTCACGATAAGGGATAGGGTAGACAGATAATAATAATGATACGATTGTCATTAAAACAGTGACATAAGCCTTAGAA  >  CP000730/824525‑824796

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