BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,930	T→A	56.1%	intergenic (‑1013/‑155)	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,930	0	T	A	56.1%	33.9 / 202.9	151	intergenic (‑1013/‑155)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (38/27); new base A (50/32); total (91/60)
Fisher's exact test for biased strand distribution p-value = 8.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

Only 100 of 154 total aligned reads displayed.

TATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTA  >  CP000730/1929790‑1930065

tataTTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCATATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttttt                                                                                                                                         <  2:504767/141‑1 (MQ=255)

 atatTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttttttt                                                                                                                                        >  2:145239/1‑140 (MQ=255)

   attcacagatggtATAATAATAGAGTCGCATATCTCTAAGCAGTCAATT‑‑GAAGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTCTTATACttt                                                                                                                                     >  1:931029/8‑140 (MQ=255)

     tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                                                                                                   >  1:739977/6‑141 (MQ=255)

      tGCCGATGATATAATAGTAGAGTCGCATATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGGCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTAtt                                                                                                                                   >  2:687741/1‑141 (MQ=255)

      cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                                  >  1:724216/5‑141 (MQ=255)

      cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                                  >  1:27522/5‑141 (MQ=255)

      cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAACTT‑‑ACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                                  >  1:485173/5‑141 (MQ=255)

       acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                                                                                 <  2:751897/138‑1 (MQ=255)

       acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                                                                                 <  1:65185/138‑1 (MQ=255)

       acagatggtATAATAAGAGCGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGGGAGGAGGTGCA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATACCATCTGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                                                                                 <  2:1009607/138‑1 (MQ=255)

        cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                                >  2:1108754/3‑141 (MQ=255)

         agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                                                                               >  1:817079/2‑141 (MQ=255)

          gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                                                                              >  1:123878/1‑140 (MQ=255)

          gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                                                                              >  1:1074319/1‑140 (MQ=255)

           atgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTTCCATTGCGTTTTTTTCTTATTGGCTa                                                                                                                              <  2:531792/141‑1 (MQ=255)

            tgttataATAAGAGACGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAATAGAGGAGGTGCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                                            <  1:222035/138‑4 (MQ=255)

            tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGGTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                            >  1:334817/4‑141 (MQ=255)

            tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                            >  2:19058/4‑141 (MQ=255)

            tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                            >  2:115399/4‑141 (MQ=255)

            tggtataAGAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                                            <  1:965861/138‑4 (MQ=255)

            tgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTATACGAAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGATGTTTTCTTGtttggtca                                                                                                                             >  2:582386/1‑135 (MQ=255)

             ggaataagaATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                           <  1:33335/138‑1 (MQ=255)

             gatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCGGAGGCCGAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAg                                                                                                                            >  1:306137/1‑141 (MQ=255)

             gatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAg                                                                                                                            >  1:857951/1‑141 (MQ=255)

              gtataATAATAGAGTCGGCTATCTCTCAGGGGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCGCCAGTCATCAGTGGCTGTGCTGTTGCGTGTTATACTTATTGGTTTAGTa                                                                                                                          >  2:256610/2‑141 (MQ=255)

              gtataATAATAGAGTCGCCTATCTCTCAGGCGTCGATT‑‑GACGAAGAGAGGAGCTGCATTAC‑T‑T‑TGCTAATAATTTTCGTTCACATCATAGCACCAATCATCAGTGGCTGTGCTGTTGCGTATTATATGTATTGGCTTAGTa                                                                                                                          >  2:348131/2‑141 (MQ=255)

              gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                                                                                          >  1:132719/2‑141 (MQ=255)

              gtataATAATAGAGTCGCCTAGCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCCTACGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                                                                                          >  1:571971/2‑141 (MQ=255)

              gtataATAATAGAGCCGCCTATCTCTCAGGAGTCAA‑TTA‑ACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTCATCTTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATAATACTTTTTGCCTCAGTa                                                                                                                          >  2:861262/2‑141 (MQ=255)

               tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                         <  2:768041/141‑2 (MQ=255)

               tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                         >  1:869104/1‑140 (MQ=255)

               tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCAGAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                         >  1:915620/1‑140 (MQ=255)

               tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCAAATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATCCTTATTGGCTTAGtca                                                                                                                         >  2:518844/1‑139 (MQ=255)

               tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGACGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                         >  2:714008/1‑140 (MQ=255)

                  aaatatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACTAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                                                                                      <  2:816243/136‑1 (MQ=255)

                    tagtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGc                                                                                                                     <  1:1004251/141‑1 (MQ=255)

                    tagtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGc                                                                                                                     <  1:216377/141‑1 (MQ=255)

                     aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                     <  2:646883/137‑1 (MQ=255)

                      atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                                                                                                  >  2:1066255/2‑140 (MQ=255)

                      atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGCAAACGCaac                                                                                                                  >  2:851496/2‑140 (MQ=255)

                       tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                                                                                 <  1:571420/141‑3 (MQ=255)

                        agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                                >  2:595827/1‑138 (MQ=255)

                        agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                                <  2:298789/141‑4 (MQ=255)

                           gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                             >  1:855455/1‑135 (MQ=255)

                            cagCCTATCTCTCAATCCTCAATTTATACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTCCTTATTGGCTAAGTAGACGCAATACaaa                                                                                                             <  1:467949/139‑1 (MQ=255)

                              gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGAA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaaaaa                                                                                                           >  2:551636/1‑140 (MQ=255)

                              gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                          <  1:741885/141‑10 (MQ=255)

                               ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                                                                         >  2:748463/1‑131 (MQ=255)

                                cTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAg                                                                                                         <  1:603150/141‑1 (MQ=255)

                                cTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAg                                                                                                         <  1:285229/141‑1 (MQ=255)

                                cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                        >  2:870350/1‑130 (MQ=255)

                                cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                        >  2:318013/1‑130 (MQ=255)

                                cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataagtaa                                                                                                        >  2:742861/1‑130 (MQ=255)

                                 tATCTCTCAGGCTTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCGCCAGTCATTATTGGCTGTGCTGTGGCTTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                       <  2:1064077/141‑13 (MQ=255)

                                 tATCTCTCAGGCGTCAATT‑‑GACGAAGGGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                       <  2:510179/141‑13 (MQ=255)

                                 tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGCTGCGTATTATACTGATGGTCTTAGTAAACGCaacaaataatcaat                                                                                                       >  2:882173/1‑129 (MQ=255)

                                  aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTATACGCaacaaataattaata                                                                                                      <  1:607034/141‑14 (MQ=255)

                                  aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                      >  1:279443/1‑128 (MQ=255)

                                  aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                      >  2:474869/1‑128 (MQ=255)

                                  aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                      >  1:361686/1‑128 (MQ=255)

                                  aTCTCTCACGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATGATTTTCGTTCACATCATAGCACCAGTCATCAGTGGGTGTGCTGTTGCGTATTATACTTAGTGGCTTAGTAAACGCaacaaataattaata                                                                                                      >  2:815114/1‑128 (MQ=255)

                                   tctctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTCTTGGCTAAGTAGACGCAATACACAATAGGta                                                                                                      >  2:897083/1‑140 (MQ=255)

                                    ctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACACATAAttaatagg                                                                                                    >  2:340398/1‑134 (MQ=255)

                                    ctctcAGGCGTCAATATAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCCATACAAAATAGGTGa                                                                                                     >  1:820861/1‑141 (MQ=255)

                                             gTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGTTTATCTTCGTTCAAATCCTATCTTCCGTGATCAGTGTCTGTGCTATCGCGTTTTTTACTTATTGTCTTACTAACCGCAATACAAAATAGGTGACATATGGcc                                                                                            <  2:458261/141‑1 (MQ=255)

                                              tCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCAATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATCGGTGACATATAGCCg                                                                                           >  1:1002530/1‑141 (MQ=255)

                                                              gaggagGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAACAAAccc                                                                           >  2:1071356/1‑141 (MQ=255)

                                                               aggaggTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTATTATTGGTTTTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCCATAAAAATcccc                                                                          >  1:526248/1‑141 (MQ=255)

                                                                ggaggTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCt                                                                         <  2:857951/141‑1 (MQ=255)

                                                                ggaggTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCCGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATCGGTGACATATAGCCGCACCAATAAAAATCCCCt                                                                         <  2:729976/141‑1 (MQ=255)

                                                                 gaggTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTc                                                                        <  1:1101970/141‑1 (MQ=255)

                                                                  aggTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCa                                                                       >  2:903794/1‑141 (MQ=255)

                                                                   ggTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCAc                                                                      <  1:1035769/141‑1 (MQ=255)

                                                                      gtA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTTGCTAAGTAGACTCAATACAGAATAGGTGACATATAGACGAACCAATAAAAATCCCCTGAcgaa                                                                   >  1:194561/1‑138 (MQ=255)

                                                                          taA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGc                                                                >  2:632019/1‑141 (MQ=255)

                                                                          taA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGAAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGc                                                                >  1:636235/1‑141 (MQ=255)

                                                                           aa‑G‑G‑TAATGCTTGTTTGCGTCCACATCCTAGCACCAGTCATCAGTGGCTTTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACACATAGCCGCACCAATAACAATCCCCTCACTACCGCa                                                               <  2:213763/141‑1 (MQ=255)

                                                                            a‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTTTCTCAGTAGACTCTATACGAATTAGTGGCCATATATCCGCCCCAATAAAAATCCCCTCACTACCGCaa                                                              >  2:231315/1‑141 (MQ=255)

                                                                            a‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCaa                                                              <  1:951174/141‑1 (MQ=255)

                                                                            a‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCaa                                                              >  2:331773/1‑141 (MQ=255)

                                                                            a‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCaa                                                              >  1:861554/1‑141 (MQ=255)

                                                                                  tGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGGTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGAAACAATAAAAATCCCCTCACTACCGCAAATa                                                           <  2:118213/141‑1 (MQ=255)

                                                                                  tGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATa                                                           <  1:719582/141‑1 (MQ=255)

                                                                                   gATGCTGATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAg                                                          <  1:390148/141‑1 (MQ=255)

                                                                                        ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgg                                                     >  2:988190/1‑141 (MQ=255)

                                                                                         tactttCGTTCACATCATCGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTGTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGCGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                                    <  2:279408/138‑1 (MQ=255)

                                                                                           tcttCGGTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                  <  1:897083/139‑1 (MQ=255)

                                                                                               cGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTgg                                              >  1:444873/1‑141 (MQ=255)

                                                                                               cGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACAGATAGCCGCACCAATAAAAATCCCCTCCCTACCGCAAAAAGTGAGGGGATTgg                                              >  1:21586/1‑141 (MQ=255)

                                                                                                  tCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgtgt                                           >  1:622748/1‑141 (MQ=255)

                                                                                                  tCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgtgt                                           >  2:183488/1‑141 (MQ=255)

                                                                                                            gCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAAt                                 <  1:342894/141‑1 (MQ=255)

                                                                                                               ccAGTCATACATGGCTGTTCCATTTCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTGAAGACt                              >  2:1072782/1‑141 (MQ=255)

                                                                                                                     atcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAtttt                        >  1:651664/1‑141 (MQ=255)

                                                                                                                     atcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACACAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGCTTGGTGTATAAGTAAATACTTAtttt                        >  1:1042337/1‑141 (MQ=255)

                                                                                                                            gCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTc                 <  2:840814/141‑1 (MQ=255)

                                                                                                                               gtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACATATTTTCGTTGTCTTa              >  1:185970/1‑141 (MQ=255)

                                                                                                                                      tGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATAc       <  1:778228/141‑1 (MQ=255)

                                                                                                                                        cGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTg     <  1:988190/141‑1 (MQ=255)

TATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTA  >  CP000730/1929790‑1930065

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