BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,946	G→A	39.5%	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	39.5%	83.8 / 91.1	108	intergenic (‑1029/‑139)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (33/31); new base A (27/15); total (62/46)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 110 total aligned reads displayed.

ATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAATAT  >  CP000730/1929806‑1930086

ataataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTAGTGGCTAGGTAg                                                                                                                                                      >  2:185483/1‑141 (MQ=255)

 taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAc                                                                                                                                                    >  2:292369/1‑141 (MQ=255)

 taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAc                                                                                                                                                    >  1:613689/1‑141 (MQ=255)

  aataGTTGGCTCGCCTGTCTCCTAGTCGTCAATTTAGACGCAGAGAGTAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGAc                                                                                                                                                    <  2:743024/141‑1 (MQ=255)

   ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                                  <  1:109307/137‑1 (MQ=255)

    tagtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCCGTGGCTGTGCCATTTCGTTTTTTTCTTATTGACTAGTTAGACGc                                                                                                                                                  >  2:416952/1‑141 (MQ=255)

    aaataggggcGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                                                                                                 <  2:5565/132‑1 (MQ=255)

      atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑ac                                                                                                                                             >  1:689720/2‑140 (MQ=255)

       tagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑At                                                                                                                                             <  2:234925/141‑1 (MQ=255)

       aagagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑aca                                                                                                                                            <  2:584958/140‑3 (MQ=255)

        agcgTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaa                                                                                                                                           <  2:721540/138‑4 (MQ=255)

        agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaa                                                                                                                                           >  1:331254/1‑138 (MQ=255)

        agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaa                                                                                                                                           <  1:70717/141‑4 (MQ=255)

        agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaa                                                                                                                                           >  1:265802/1‑138 (MQ=255)

        agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaa                                                                                                                                           >  1:504381/1‑138 (MQ=255)

         gattCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaa                                                                                                                                          <  2:62385/138‑5 (MQ=255)

         gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaa                                                                                                                                          >  1:775920/1‑137 (MQ=255)

         gagTCGCCTATCTCTCAGGCGTAAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaa                                                                                                                                          <  2:46574/141‑5 (MQ=255)

          agTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACa                                                                                                                                          >  2:88928/1‑141 (MQ=255)

           gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaata                                                                                                                                        >  1:400650/1‑135 (MQ=255)

           gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaata                                                                                                                                        <  1:527599/141‑7 (MQ=255)

            gcggcttactctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACaa‑‑a                                                                                                                                      <  2:785247/133‑1 (MQ=255)

             cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataat                                                                                                                                      >  1:314514/1‑133 (MQ=255)

             cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataat                                                                                                                                      >  1:768522/1‑133 (MQ=255)

              gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                                     <  1:342295/141‑10 (MQ=255)

              gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                                     >  2:185601/1‑132 (MQ=255)

               gatATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatta                                                                                                                                    >  2:7948/3‑131 (MQ=255)

               ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATa                                                                                                                                   <  1:55901/141‑1 (MQ=255)

               ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatta                                                                                                                                    >  2:65493/1‑131 (MQ=255)

                cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaccaaatgattaa                                                                                                                                     >  1:173919/1‑129 (MQ=255)

                cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaa                                                                                                                                   >  2:150269/1‑130 (MQ=255)

                cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaa                                                                                                                                   <  1:71526/141‑12 (MQ=255)

                 tATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTTTTTTTTGTTAAGTAGACGCA‑‑ATACAG‑‑AATAgg                                                                                                                                 >  2:507565/1‑141 (MQ=255)

                 tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaat                                                                                                                                  >  2:710614/1‑129 (MQ=255)

                 tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaat                                                                                                                                  >  2:732980/1‑129 (MQ=255)

                 tATCTCTCAGGCGTAAATGTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAgg                                                                                                                                 <  2:614323/141‑1 (MQ=255)

                  aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGATGTGCA‑TTA‑C‑TTTGCTAATCCTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                                 <  2:422526/141‑14 (MQ=255)

                  aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                                 >  2:108745/1‑128 (MQ=255)

                  aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                                 >  2:176229/1‑128 (MQ=255)

                  aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                                 >  2:421488/1‑128 (MQ=255)

                  aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                                 >  1:15817/1‑128 (MQ=255)

                  aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTACTGGCTTAGTGAACGCa‑‑acaaataattaata                                                                                                                                 >  2:500437/1‑128 (MQ=255)

                  aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                                 <  2:422755/141‑14 (MQ=255)

                  aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaatacttaata                                                                                                                                 >  1:539569/1‑128 (MQ=255)

                  aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                                 >  2:217934/1‑128 (MQ=255)

                   tttctcAGGCATCAATTTAGACGAAGAGAGGATGTCTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTg                                                                                                                               <  1:46378/139‑1 (MQ=255)

                   tctctcAGGCGTCAATTT‑‑ACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaatag                                                                                                                                <  2:510007/141‑15 (MQ=255)

                   tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaatag                                                                                                                                <  2:610091/141‑15 (MQ=255)

                   tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑aaaaataattaatag                                                                                                                                >  2:607126/1‑127 (MQ=255)

                        cAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑a‑ta                                                                                                                        <  1:426748/141‑1 (MQ=255)

                        cAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑a‑ta                                                                                                                        <  1:671500/141‑1 (MQ=255)

                         aagCGTCACTTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTATTGCGTTTTTTTTTTTTTTTTTAAGTAGACGCA‑‑CGACAA‑‑AATAGGTGAA‑a‑tat                                                                                                                       >  2:582805/3‑141 (MQ=255)

                         aGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCGATTGCGTTTTTTTCTTATTGGCTCAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑a‑tat                                                                                                                       >  2:582931/1‑141 (MQ=255)

                         aGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTAACTGGAGATGTGAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAC‑‑AATAGGTGAC‑a‑tat                                                                                                                       <  1:413639/141‑1 (MQ=255)

                           gCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCGCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTACACGCA‑‑GTACAA‑‑AATAGGTGAC‑A‑TAT‑Ag                                                                                                                    >  1:255882/1‑141 (MQ=255)

                              tCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGCC‑A‑TAT‑AGccc                                                                                                                 >  2:424931/1‑140 (MQ=255)

                              tCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑g                                                                                                                >  1:543142/1‑141 (MQ=255)

                              tCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑g                                                                                                                >  2:203661/1‑141 (MQ=255)

                                     aGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAAt                                                                                                         >  2:419533/1‑141 (MQ=255)

                                     aGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTACTTTCATTCATATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAAt                                                                                                         <  1:534310/141‑1 (MQ=255)

                                      gACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCA‑GCACCAATa                                                                                                        >  2:233968/1‑141 (MQ=255)

                                             agagGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCCCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATcc                                                                                                 >  1:25230/1‑141 (MQ=255)

                                               aggaggTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATcccc                                                                                               >  1:80486/1‑141 (MQ=255)

                                                     tgtA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAGTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAA‑AGCC‑GCACCAATAAAAATCCCCTCACTa                                                                                         >  1:24220/1‑141 (MQ=255)

                                                        a‑taA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCg                                                                                      <  1:497200/141‑1 (MQ=255)

                                                          taA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGc                                                                                     <  2:651966/141‑1 (MQ=255)

                                                           aa‑G‑G‑TGATGCTTATTTTCGTTCACATCGTAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCa                                                                                    <  2:511883/141‑1 (MQ=255)

                                                           aa‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATACCCTCACTACCGCa                                                                                    >  2:20756/1‑141 (MQ=255)

                                                                g‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACCA‑‑AATCGATGAC‑C‑TAT‑AGCA‑GCACACAGCAACATCCCCTCACTGCCGccata                                                                                 >  2:541548/1‑137 (MQ=255)

                                                                       gttATTTTCGTTCACATCGTAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAg                                                                           <  2:73456/140‑1 (MQ=255)

                                                                       cTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTTTTTTTTTTTAAGTAGACGCA‑‑ATACCA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCGCCAATAAAAATCCCCTCACTACCGCAAATAGGGAg                                                                           >  2:661488/1‑141 (MQ=255)

                                                                       cTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTTTTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCCCACTCCCGCAAATAGTGAg                                                                           >  2:423313/1‑141 (MQ=255)

                                                                       cTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAg                                                                           <  2:120961/141‑1 (MQ=255)

                                                                          atattCGGTCACGTCGTAGGGGCAGTCATCAGGGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAG‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgggg                                                                        <  1:270743/138‑1 (MQ=255)

                                                                          aTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCCAATAGGGAgggg                                                                        >  2:122738/1‑141 (MQ=255)

                                                                          aTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AACAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgggg                                                                        <  1:605459/141‑1 (MQ=255)

                                                                                 ttCACATCATAGCACCAGTCATCAGTGGGTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAC‑‑AATAGGTGAC‑A‑TAT‑AGAC‑CCACCAATAAAAATCCCCTCACTACCGCAAATAGGGAGGGGGTTGgtg                                                                 >  1:21666/1‑141 (MQ=255)

                                                                                     catcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTATTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGtata                                                             >  2:422593/1‑141 (MQ=255)

                                                                                        catAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATACCCTCACTACCGCAAATCGTGAGGGGGTTGGTGTATAAGt                                                          >  1:161611/1‑141 (MQ=255)

                                                                                         atAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTa                                                         >  1:439671/1‑141 (MQ=255)

                                                                                          tAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTaa                                                        >  1:628677/1‑141 (MQ=255)

                                                                                           aGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGAC‑TTTATAAGCCAGCATCGA‑AAAAATCCCCACACTACTGC‑AATAGTGTGGGGATTGGTGcatta                                                            >  2:642208/1‑139 (MQ=255)

                                                                                            gCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTT‑TATAAGCCAGCATCGA‑AAAAATCCCCACACTACTGC‑AATAGTGTGGGGATTGgtgcattaa                                                           >  2:737388/1‑135 (MQ=255)

                                                                                              actaGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCGCCAATAAAAATTCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATAc                                                    <  1:227747/138‑1 (MQ=255)

                                                                                              aCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACCA‑‑CTTAGGGGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATAc                                                    <  1:574666/141‑1 (MQ=255)

                                                                                               ccAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCATATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACt                                                   <  1:216308/141‑1 (MQ=255)

                                                                                                   tcatcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAtt                                               <  1:736978/141‑1 (MQ=255)

                                                                                                    catcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATCCAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAttt                                              >  1:171201/1‑141 (MQ=255)

                                                                                                    catcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAttt                                              <  1:349128/141‑1 (MQ=255)

                                                                                                    catcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAttt                                              <  2:671676/141‑1 (MQ=255)

                                                                                                    catcaGTGGCCGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTCCCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAttt                                              <  1:486741/141‑1 (MQ=255)

                                                                                                     atcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAtttt                                             >  2:762766/1‑141 (MQ=255)

                                                                                                        aGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGt                                          <  1:21626/141‑1 (MQ=255)

                                                                                                            gCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTc                                      <  2:80486/141‑1 (MQ=255)

                                                                                                               gtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTa                                   >  1:277087/1‑141 (MQ=255)

                                                                                                                 gCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TGT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAAt                                 >  1:196011/1‑141 (MQ=255)

                                                                                                                      tGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATAc                            <  2:525578/141‑1 (MQ=255)

                                                                                                                                 cTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAAttttt                 >  2:583529/1‑141 (MQ=255)

                                                                                                                                   tATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCt               <  1:403482/141‑1 (MQ=255)

                                                                                                                                                gACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAatat  <  1:790969/141‑1 (MQ=255)

ATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑A‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAATAT  >  CP000730/1929806‑1930086

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