BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,946	G→A	19.3%	intergenic (‑1029/‑139)	USA300HOU_1794 ← / → USA300HOU_1796	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,929,946	0	G	A	19.3%	321.6 / 64.3	145	intergenic (‑1029/‑139)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (68/49); new base A (13/15); total (81/64)
Fisher's exact test for biased strand distribution p-value = 2.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

Only 100 of 145 total aligned reads displayed.

AATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAATA  >  CP000730/1929808‑1930085

 ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                                <  2:939624/137‑1 (MQ=255)

 ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                                <  2:1369602/137‑1 (MQ=255)

 ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                                <  2:1101199/137‑1 (MQ=255)

  tagtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGTCTAAGTAGACGc                                                                                                                                                >  2:650653/1‑141 (MQ=255)

     tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                                                                                                            >  2:907281/1‑139 (MQ=255)

      agagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑a                                                                                                                                          >  1:1503160/1‑141 (MQ=255)

      agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                                                           >  1:483866/1‑138 (MQ=255)

      agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                                                           <  2:1170798/141‑4 (MQ=255)

      agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCA‑TTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                                                           >  1:493091/1‑138 (MQ=255)

       gagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑Ac                                                                                                                                         <  1:1643263/141‑1 (MQ=255)

        agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                                                         >  1:659176/1‑136 (MQ=255)

        agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                                                         <  1:333207/141‑6 (MQ=255)

         gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                                                        <  1:1116210/141‑7 (MQ=255)

         gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                                                        >  2:1063340/1‑135 (MQ=255)

         gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGATGGGG‑TAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                                                        >  1:798539/1‑135 (MQ=255)

          tCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACaa‑‑a                                                                                                                                    <  1:1464183/141‑1 (MQ=255)

           cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCA‑TTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                                                                                                                                      >  1:476466/1‑133 (MQ=255)

            gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGCGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AAt                                                                                                                                  >  2:1102600/1‑141 (MQ=255)

            gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AAt                                                                                                                                  >  2:1452066/1‑141 (MQ=255)

            gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                                                     <  2:377700/141‑10 (MQ=255)

            gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                                                     <  1:107660/141‑10 (MQ=255)

            gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                                                     <  1:1499478/141‑10 (MQ=255)

              cTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AAAAg                                                                                                                                >  1:828419/1‑141 (MQ=255)

              cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                                                   >  1:1433287/1‑130 (MQ=255)

              cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCA‑TTACTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGGAAACGCaacaaataattaa                                                                                                                                   >  1:389714/1‑130 (MQ=255)

              cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCA‑TTA‑CTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                                                   <  2:434342/141‑12 (MQ=255)

                aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑TT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGGAAACGCaacaaataattaata                                                                                                                                 <  2:389714/141‑14 (MQ=255)

                       aGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑atat                                                                                                                      <  1:1058421/141‑1 (MQ=255)

                         gCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑Ag                                                                                                                   >  2:398026/1‑141 (MQ=255)

                            tCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAC‑‑AATAGGTAAC‑ATAT‑ACCC‑g                                                                                                               >  2:593632/1‑141 (MQ=255)

                             cAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑Gc                                                                                                              >  1:1014603/1‑141 (MQ=255)

                             cAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑Gc                                                                                                              >  1:234803/1‑141 (MQ=255)

                              aaTTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGCCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCa                                                                                                             <  2:1303011/141‑1 (MQ=255)

                                tttAGACGCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCAcc                                                                                                           <  2:82037/141‑1 (MQ=255)

                                 ttAGACGCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCa                                                                                                          >  2:785120/1‑141 (MQ=255)

                                 ttAGACGCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCa                                                                                                          <  1:616746/141‑1 (MQ=255)

                                  tAGACGCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCaa                                                                                                         >  2:307028/1‑141 (MQ=255)

                                   aGACGCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAAt                                                                                                        >  2:1034142/1‑141 (MQ=255)

                                   aGACGCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAAt                                                                                                        <  2:184202/141‑1 (MQ=255)

                                   aGACGCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAAt                                                                                                        <  2:1503160/141‑1 (MQ=255)

                                    gACGCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATa                                                                                                       <  2:219729/141‑1 (MQ=255)

                                      cGCAGAGAGGAGGT‑GTA‑TAA‑GT‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATaaa                                                                                                     <  2:148404/141‑1 (MQ=255)

                                       gCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATaaaa                                                                                                    >  1:1089195/1‑141 (MQ=255)

                                       gCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATaaaa                                                                                                    >  1:71364/1‑141 (MQ=255)

                                       gCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATaaaa                                                                                                    <  1:1004352/141‑1 (MQ=255)

                                                       a‑taA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCg                                                                                     <  2:778769/141‑1 (MQ=255)

                                                             gg‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTatttgcgg                                                                                 <  1:433283/141‑8 (MQ=255)

                                                             gg‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCaaa                                                                                 <  2:273451/141‑1 (MQ=255)

                                                             gg‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCaaa                                                                                 <  1:1002215/141‑1 (MQ=255)

                                                              g‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTatttgcggt                                                                                <  2:310055/141‑9 (MQ=255)

                                                                  aTGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCCCTACCGCAAATAGt                                                                             >  2:1552005/1‑141 (MQ=255)

                                                                     cTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACaa‑‑aa                                                                                                                                   <  1:1587018/87‑1 (MQ=255)

                                                                     cTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACaa‑‑aa                                                                                                                                   >  2:1587018/1‑87 (MQ=255)

                                                                     cTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAg                                                                          <  1:24189/141‑1 (MQ=255)

                                                                       tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                                                        <  1:633162/141‑1 (MQ=255)

                                                                        aTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgggg                                                                       >  2:276834/1‑141 (MQ=255)

                                                                         ttttCGTTCGCATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                                      <  2:893844/141‑1 (MQ=255)

                                                                         ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGTCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTATCGCAAATAGTGAGGGGa                                                                      <  2:1382912/141‑1 (MQ=255)

                                                                         ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                                      >  1:1374596/1‑141 (MQ=255)

                                                                         ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                                      <  2:1357559/141‑1 (MQ=255)

                                                                         ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                                      <  1:1315681/141‑1 (MQ=255)

                                                                           ttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGAtt                                                                    >  1:898666/1‑141 (MQ=255)

                                                                           ttCGTTCACATCATAGCACCAGTAATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGAtt                                                                    <  1:180279/141‑1 (MQ=255)

                                                                            tCGTTCACATCATAGCCCCATCCATCCGTGGTTTTGCCATGTCGGTGGTTTTTTATTTTCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTg                                                                   <  1:1562916/141‑1 (MQ=255)

                                                                               ttCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTTGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgtg                                                                >  2:843509/1‑141 (MQ=255)

                                                                               ttCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgtg                                                                >  1:76360/1‑141 (MQ=255)

                                                                                   catcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGtata                                                            <  1:307028/141‑1 (MQ=255)

                                                                                    atcatAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGA‑AAAAATCCCCACACTACTGC‑AATAGTGCGGGGATTGgtg                                                                >  2:1641040/1‑141 (MQ=255)

                                                                                     tcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAg                                                          <  1:1484012/141‑1 (MQ=255)

                                                                                     tcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAg                                                          >  2:743996/1‑141 (MQ=255)

                                                                                        tAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTaa                                                       >  2:551125/1‑141 (MQ=255)

                                                                                        tAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTaa                                                       >  1:1625092/1‑141 (MQ=255)

                                                                                        tAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTaa                                                       >  1:670369/1‑141 (MQ=255)

                                                                                           cACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAt                                                                                                                                             <  2:1492543/59‑1 (MQ=255)

                                                                                           cACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAt                                                                                                                                             >  1:1492543/1‑59 (MQ=255)

                                                                                           cACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATa                                                    >  2:1654052/1‑141 (MQ=255)

                                                                                            aCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATATTGAGGGGATTGGTGTATAAGTAAATAc                                                   <  1:386347/141‑1 (MQ=255)

                                                                                            aCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATAc                                                   <  2:627953/141‑1 (MQ=255)

                                                                                            aCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATAc                                                   <  1:781169/141‑1 (MQ=255)

                                                                                              cAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTTTTTTTTTTTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACtt                                                 >  2:1303069/1‑141 (MQ=255)

                                                                                                   atcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAtttt                                            <  1:759162/141‑1 (MQ=255)

                                                                                                    tccgTGGCTGTGCCATTGCGTTTTTTTCTTATTTGTTTAGTAGACGCGAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATTGTGAGGGAGATGGTGTAAAAGTATGTACGTAtgtgt                                           >  2:745494/4‑137 (MQ=255)

                                                                                                       gTGGCTGTGCCATTGCGTTTTTTTTTTTTTTTTTTAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCCCCAATAACAATCCCCTCACTACCGCAAATAGTTAGGGGATTGGTGTATACTTAAATTCTTTTTTTCGtt                                        >  2:845157/1‑141 (MQ=255)

                                                                                                        tGGCTGTGCCATTGCGTTTTTTTTTTTTTTGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTTTTTTCGTTg                                       >  2:503287/1‑141 (MQ=255)

                                                                                                          gCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGCGGATTGTTGTATAAGTAAATACTTATTTTCGTTGTc                                     <  1:793150/141‑1 (MQ=255)

                                                                                                             gtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTa                                  >  2:769257/1‑141 (MQ=255)

                                                                                                              tgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGc                                                                                                                  <  2:1185039/61‑1 (MQ=255)

                                                                                                              tgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGc                                                                                                                  >  1:1185039/1‑61 (MQ=255)

                                                                                                               gCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAAt                                >  2:822892/1‑141 (MQ=255)

                                                                                                                        tttttttCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCt                       <  2:350340/141‑1 (MQ=255)

                                                                                                                              tttttttttttAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTTAC‑ATAT‑GGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGTGTTTTGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAAtttt                 >  2:1321120/11‑141 (MQ=255)

                                                                                                                                ttATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTc               >  1:624558/1‑141 (MQ=255)

                                                                                                                                ttATTGGCAAAGTAGACGCAAT‑‑ACCA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATACCCTTCCTACCGCAAATAGGGAGGGGGTCGGTGTATAAGTGAATACTTATTTTCGTTGTGTTCACTCTACTGCttcttcgta               >  2:93127/1‑133 (MQ=255)

                                                                                                                                  tttGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCtt             <  2:1127000/140‑1 (MQ=255)

                                                                                                                                  aTTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCtt             >  1:1488464/1‑141 (MQ=255)

                                                                                                                                      gCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTg         >  2:615529/1‑141 (MQ=255)

                                                                                                                                         aaGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAAttaaaatttttataa      <  2:562697/141‑14 (MQ=255)

                                                                                                                                         aaGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTaa      <  2:427389/141‑1 (MQ=255)

                                                                                                                                           gTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTaaaa    >  1:1384142/1‑141 (MQ=255)

                                                                                                                                             aGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAata  <  1:398026/141‑1 (MQ=255)

AATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑GG‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAATA  >  CP000730/1929808‑1930085

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