BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,941	A→T	50.5%	intergenic (‑1024/‑144)	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,941	0	A	T	50.5%	6.6 / 196.0	137	intergenic (‑1024/‑144)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (32/36); new base T (44/25); total (76/61)
Fisher's exact test for biased strand distribution p-value = 5.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

Only 100 of 139 total aligned reads displayed.

ATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTG  >  CP000730/1929801‑1930078

 tgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTCCGTTTTTTTCTTATTGGCCaa                                                                                                                                                 >  1:169448/1‑141 (MQ=255)

  ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                                               >  1:612586/3‑141 (MQ=255)

  ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                                               >  1:204431/3‑141 (MQ=255)

  ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                                               <  2:930193/139‑1 (MQ=255)

  ggtacaataatAGCGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTAGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                                               <  2:566423/136‑1 (MQ=255)

  gatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAg                                                                                                                                                <  1:649515/141‑1 (MQ=255)

   gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                                                                                                              >  1:139950/2‑141 (MQ=255)

     ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                                                            >  2:677812/1‑139 (MQ=255)

     ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                                                            >  1:4748/1‑139 (MQ=255)

      taataagagcggcGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAc                                                                                                                                           <  1:42065/129‑1 (MQ=255)

      taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                                                                                                                           <  1:772575/141‑4 (MQ=255)

        ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTCC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCGGTCAGCAGCGGCTTTGCTGTTGCGTATTATTCTTATTGGCTTAGTAAACGc                                                                                                                                         >  2:198069/5‑141 (MQ=255)

        ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                         <  1:184848/141‑1 (MQ=255)

        ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                         <  1:38992/141‑1 (MQ=255)

        ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATAATACGGATTGGCTTAGTAAACGc                                                                                                                                         <  1:226877/141‑1 (MQ=255)

        ataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACg                                                                                                                                          >  2:70133/1‑141 (MQ=255)

           gtagAGACGACTATCTCTCAGGCGCCAAGTTAGACGCAGAGGGGATGTG‑TA‑TAA‑G‑GGGATTCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCa‑‑a                                                                                                                                     <  2:414113/141‑1 (MQ=255)

           atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCTGTTGCTTTGCTGTTGCGTATTATATTTATTGGCTTATTAAACGCa‑‑ac                                                                                                                                    >  2:147231/2‑140 (MQ=255)

           atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑ac                                                                                                                                    >  2:633990/2‑140 (MQ=255)

            tagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑GTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTACTTAGACGCA‑‑At                                                                                                                                    >  2:259973/1‑141 (MQ=255)

            tagAGTCGCCTATCTCTCAGGCGTCAATGTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑At                                                                                                                                    <  1:346058/141‑1 (MQ=255)

             agagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATa                                                                                                                                   <  1:913395/141‑1 (MQ=255)

             agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaa                                                                                                                                  >  1:864551/1‑138 (MQ=255)

             agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaa                                                                                                                                  >  2:682593/1‑138 (MQ=255)

              gagTCGCCTATCTTTCATGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaa                                                                                                                                 <  2:378383/141‑5 (MQ=255)

              gagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATAc                                                                                                                                  <  2:928912/141‑1 (MQ=255)

              gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCAGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAGGCa‑‑acaaa                                                                                                                                 >  1:68318/1‑137 (MQ=255)

              gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑agaat                                                                                                                                 >  1:442004/1‑137 (MQ=255)

               agTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACa                                                                                                                                 >  1:707462/1‑141 (MQ=255)

                gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaata                                                                                                                               <  1:716526/141‑7 (MQ=255)

                gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACAAGACATCAGTTGCTGTGCTGTTGCGTATTATACTTGTTGGCTTAGTAAACGCa‑‑acaaata                                                                                                                               <  1:689121/141‑7 (MQ=255)

                  ggCCTATCTCTCCGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTC‑C‑TT‑‑TGCTAATCATTTCGTTCACATCATAGCACCACTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGTCTTAGTAAACGCa‑‑acaaataatt                                                                                                                            <  2:337322/140‑10 (MQ=255)

                  cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataat                                                                                                                             >  2:272426/1‑133 (MQ=255)

                  cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataat                                                                                                                             >  2:139467/1‑133 (MQ=255)

                   gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                            <  2:718271/141‑10 (MQ=255)

                   gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                            <  2:134626/141‑10 (MQ=255)

                   gCCTATCTCTCAGGCGGCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                            <  1:49383/141‑10 (MQ=255)

                   gCCGATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TTA‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AAt                                                                                                                           <  2:494823/141‑1 (MQ=255)

                    ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGCAGACGCA‑‑GTACaa‑‑aaga                                                                                                                          >  2:891759/1‑139 (MQ=255)

                    ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatta                                                                                                                           >  2:945098/1‑131 (MQ=255)

                     cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaa                                                                                                                          <  2:864551/141‑12 (MQ=255)

                     cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaa                                                                                                                          >  1:769679/1‑130 (MQ=255)

                      tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCAA‑TAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGTCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaaa                                                                                                                         >  1:526362/1‑129 (MQ=255)

                       ttctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTCCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                        <  2:695177/140‑14 (MQ=255)

                       aTCTCTCAGGCGTCACTTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGt                                                                                                                       >  1:937040/1‑141 (MQ=255)

                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGTTTAGTAAACGCa‑‑acaaataattgata                                                                                                                        >  1:714068/1‑128 (MQ=255)

                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaatg                                                                                                                        >  1:699730/1‑128 (MQ=255)

                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                        >  1:167979/1‑128 (MQ=255)

                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                        >  1:632144/1‑128 (MQ=255)

                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                        <  2:800637/141‑14 (MQ=255)

                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                        >  1:299714/1‑128 (MQ=255)

                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGACTTAGTAAACGCa‑‑acaaataattaata                                                                                                                        >  1:283630/1‑128 (MQ=255)

                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAATGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                        >  1:542668/1‑128 (MQ=255)

                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTGCACATCATAGCACCAGTCATGAGTGGCTGTGCTGTTGCGTATTATACTGATTGGCTTAGTAGACGCa‑‑acaaataattaaaa                                                                                                                        >  1:95698/1‑128 (MQ=255)

                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattcata                                                                                                                        >  2:327059/1‑128 (MQ=255)

                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGATAATCAT‑‑TATCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGCTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                        >  2:649384/1‑128 (MQ=255)

                       aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTTTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                        >  1:839580/1‑128 (MQ=255)

                        tctctcATGGGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaatag                                                                                                                       <  2:365915/141‑15 (MQ=255)

                        tctctcAGGCGTCAATTTAGACGCAGAGAGGCGGCG‑TA‑CAA‑G‑GTGATGCTTAT‑‑TTTAGTTCACAACATAGAACCAGTCATCAGTGGATGTGCCATTGCGTTTTTTGCTTATTGGCTAAGTTGGCGCA‑‑ATCCAA‑‑ACTAGGTg                                                                                                                      >  2:453324/1‑141 (MQ=255)

                        tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaatag                                                                                                                       >  1:553255/1‑127 (MQ=255)

                             cAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AGACAA‑‑AATAGGTGAC‑ata                                                                                                                >  1:736943/1‑141 (MQ=255)

                             cAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCAGAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ata                                                                                                                <  2:428250/141‑1 (MQ=255)

                              atgCGTAAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑GTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGGGCCATTGCGGTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑atat                                                                                                               <  1:550462/139‑1 (MQ=255)

                               ggCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑GTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑Atat‑a                                                                                                             <  2:398467/141‑1 (MQ=255)

                                gCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTTGATGCA‑‑GTACAA‑‑AATAGGTGAC‑ATAT‑Ag                                                                                                            >  1:679796/1‑141 (MQ=255)

                                 cGTTAATTTAGACGCAGGGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGc                                                                                                           <  1:309534/141‑1 (MQ=255)

                                  gTCAATTTTGACGCAGGGAGGAAGTT‑GA‑TAT‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGGGCTATTTCGTTTTTTTCTTATTGGCTAAGTATACGCA‑‑ATACAA‑‑TATAGGTGAC‑ATAT‑AGcc                                                                                                          <  2:90529/141‑1 (MQ=255)

                                     aaTTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCa                                                                                                      <  1:101620/141‑1 (MQ=255)

                                       tttAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GGGATGCTTAT‑‑GTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGAA‑‑ATACAC‑‑CATAGGAGAC‑ATAT‑AGCC‑GCAcc                                                                                                    <  2:693638/141‑1 (MQ=255)

                                         tAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑CTACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCCCCaa                                                                                                  >  1:377853/1‑141 (MQ=255)

                                           gACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATa                                                                                                <  1:883763/141‑1 (MQ=255)

                                           gACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATa                                                                                                <  1:781378/141‑1 (MQ=255)

                                             cGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATaaa                                                                                              <  1:725120/141‑1 (MQ=255)

                                                 gagagGAGGTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATc                                                                                          >  1:441603/1‑141 (MQ=255)

                                                     ggaggTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCt                                                                                      >  1:348371/1‑141 (MQ=255)

                                                         gtg‑tA‑TAA‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCAca                                                                                  >  2:185922/1‑140 (MQ=255)

                                                                  a‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCaa                                                                           >  1:757352/1‑141 (MQ=255)

                                                                     ttttaTGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTCGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCaaa                                                                          >  2:222905/5‑141 (MQ=255)

                                                                          tGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTg                                                                     <  2:500526/141‑1 (MQ=255)

                                                                            cTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAACAGTGAg                                                                   <  2:520782/141‑1 (MQ=255)

                                                                              tAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                                                 >  1:599473/1‑141 (MQ=255)

                                                                              tAT‑‑TTTCGGTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTGTTTTTCTTATTGCCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                                                 <  2:645611/141‑1 (MQ=255)

                                                                                        ttCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgtg                                                         >  1:623820/1‑141 (MQ=255)

                                                                                         taacaTCAAAGCACCAGCCATCAGTGGCTGTGCCAATGCGTTTTTTTCTTATTGGCTAAGCAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgtgt                                                        <  2:265796/139‑1 (MQ=255)

                                                                                                atAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTTTTTTCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTa                                                 >  2:78495/1‑141 (MQ=255)

                                                                                                   gCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGTGACTTTATAAGCCAGCAACGA‑AAAAATCCCCACACTACTGC‑AATAGTGTGGGGATTGGTGcattaa                                                   >  2:690474/1‑138 (MQ=255)

                                                                                                     cccAGTCATCAGTGGCTGGGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATAc                                            <  1:495684/140‑1 (MQ=255)

                                                                                                           catcaGTGGCTGTGCCATTGCGTTTTTTTCTTGTTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAttt                                      <  2:22524/141‑1 (MQ=255)

                                                                                                               aGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGt                                  >  1:11169/1‑141 (MQ=255)

                                                                                                               aGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGt                                  >  2:573681/1‑141 (MQ=255)

                                                                                                                      gtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTa                           <  1:677448/141‑1 (MQ=255)

                                                                                                                      gtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTa                           >  2:286405/1‑141 (MQ=255)

                                                                                                                      gtgCCATTGCGTTTTTTTCTTATTGGCTAAGGAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTa                           <  1:829688/141‑1 (MQ=255)

                                                                                                                                  ttttttttTATTGGATAAGTAGACGCA‑‑ATACAT‑‑AATAGGTGAC‑ACAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTa               <  2:377853/141‑1 (MQ=255)

                                                                                                                                     tttCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAAtt            <  1:729231/141‑1 (MQ=255)

                                                                                                                                     tttATTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAAtt            <  2:747109/141‑1 (MQ=255)

                                                                                                                                       tCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAAtttt          <  1:733059/141‑1 (MQ=255)

                                                                                                                                        cTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAAttttt         >  1:464843/1‑141 (MQ=255)

                                                                                                                                          taattgCTAAGTAAACGCA‑‑ATACAA‑‑AATAGGGGAC‑ATAT‑AGCC‑GCGCTGAGGAAAATCGCCATACTATCGCATATTTTCAGGGGATTGATGTATTATGAAATACTTATTTTCCTTGTCTTAATTATACTGCTAATATTTCt       <  2:394821/136‑1 (MQ=255)

                                                                                                                                               gCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTg  <  2:937040/141‑1 (MQ=255)

ATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTG  >  CP000730/1929801‑1930078

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