BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,929,848	0	C	A	75.9%	782.8 / inf	500	intergenic (‑931/‑237)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (58/62); new base A (192/185); total (252/248)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-02
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 503 total aligned reads displayed.

TAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAA  >  CP000730/1929708‑1929988

                        aCGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTCTCTATCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATg                                                                                                                      <  2:1101758/141‑1 (MQ=255)

                        aCGATTCAATCTATCAGCTAGAATAATATAACGTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATg                                                                                                                      <  1:991895/141‑1 (MQ=255)

                         cGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGGG‑TA‑TAA‑G‑GTGATGc                                                                                                                     >  1:895352/1‑141 (MQ=255)

                                     tCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATCATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTc                                                                                                         >  1:163959/1‑141 (MQ=255)

                                                cATATAACTTATCAGACGCACGTGGCGCGTCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCa                                                                                              >  1:583086/1‑141 (MQ=255)

                                                 atatAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCAc                                                                                             >  2:905817/1‑141 (MQ=255)

                                                              atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                               <  2:308531/127‑1 (MQ=255)

                                                               tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCACCAGTGGCt                                                                              >  2:846613/14‑141 (MQ=255)

                                                                 ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                            <  1:910644/130‑1 (MQ=255)

                                                                  gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                           <  1:1083004/131‑1 (MQ=255)

                                                                  gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                           <  1:109766/131‑1 (MQ=255)

                                                                  gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                           <  1:329649/131‑1 (MQ=255)

                                                                  gaagcagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                           <  1:439366/131‑1 (MQ=255)

                                                                   aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                          <  1:339545/132‑1 (MQ=255)

                                                                    agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                                                         <  1:677263/133‑2 (MQ=255)

                                                                    aagagatatgCTTGTATATTCACAGATCGTATAATAATAGAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCGCCAGTCATCAGTGGCTGTGct                                                                         <  2:37459/133‑2 (MQ=255)

                                                                     gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTgcgcgg                                                                        <  2:142043/134‑3 (MQ=255)

                                                                     gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                        >  1:419062/8‑139 (MQ=255)

                                                                     gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                        <  2:950646/134‑3 (MQ=255)

                                                                     gtagatatgCTTGTATATTCACAGATGGAATAATAATAGAGTCGCCTAACTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCATTTTCGTTCACATCACAGCACCAGTCATCAGTGGCTGTGctg                                                                        >  1:728846/8‑139 (MQ=255)

                                                                      tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                       <  2:561478/135‑4 (MQ=255)

                                                                       agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgtt                                                                      >  2:200616/6‑137 (MQ=255)

                                                                       agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgtt                                                                      <  1:653811/136‑5 (MQ=255)

                                                                       agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgtt                                                                      >  2:150655/6‑137 (MQ=255)

                                                                        gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                     >  1:505757/5‑136 (MQ=255)

                                                                        gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                     <  1:316134/137‑6 (MQ=255)

                                                                        gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                     >  1:691440/5‑136 (MQ=255)

                                                                        gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAATGGCTGGGctgttg                                                                     >  2:1089659/5‑136 (MQ=255)

                                                                        gataagCTTGTATATTCACAGATGGTATAATAATAGGGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATACGTGGCTGTGctgttg                                                                     >  2:1009697/6‑136 (MQ=255)

                                                                         atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                                                                    >  1:196253/4‑135 (MQ=255)

                                                                         atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCATCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                                                                    >  1:801572/4‑135 (MQ=255)

                                                                          tgtgCTCGTATACTTGCCGATGTTATAATAGAAGAGTCGCCTATCTCTCAGTCGTCAATTTAGACGCAGAGAGGAGGTG‑GA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGc                                                                    <  2:245263/141‑1 (MQ=255)

                                                                          tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                                                                   >  1:1206117/3‑141 (MQ=255)

                                                                          tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCGATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                                                                   >  1:824611/3‑141 (MQ=255)

                                                                              cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                               <  1:388032/141‑4 (MQ=255)

                                                                                tGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttt                                                              <  2:585974/141‑1 (MQ=255)

                                                                                 gTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttttt                                                             <  2:831638/141‑1 (MQ=255)

                                                                                  tataTTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttttt                                                            <  2:180408/141‑1 (MQ=255)

                                                                                  tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                                                           <  1:451677/141‑8 (MQ=255)

                                                                                  tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                                                           <  2:710299/141‑8 (MQ=255)

                                                                                  tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                                                           <  2:26203/141‑8 (MQ=255)

                                                                                  tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAAAT‑‑GAATTGGTGTGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                                                           <  1:275746/141‑8 (MQ=255)

                                                                                   atatTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTGCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTc                                                           <  2:349701/141‑1 (MQ=255)

                                                                                    tatTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCt                                                          >  1:856711/1‑141 (MQ=255)

                                                                                     attcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatactta                                                        <  2:1061553/134‑8 (MQ=255)

                                                                                      ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                       >  1:1101130/7‑141 (MQ=255)

                                                                                      ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                       >  1:745376/7‑141 (MQ=255)

                                                                                       tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                      <  1:534030/136‑1 (MQ=255)

                                                                                        cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                     >  1:130881/5‑141 (MQ=255)

                                                                                         acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                    <  1:1000429/138‑1 (MQ=255)

                                                                                         acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                    <  1:361485/138‑1 (MQ=255)

                                                                                          cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                   >  1:1026034/3‑141 (MQ=255)

                                                                                           cGATGATATTATAGTCGAGTCGCCTATCTCTCAGTCTTCAATTTAGGCGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc                                                   <  2:1175439/141‑1 (MQ=255)

                                                                                            ggtggtataATAATAGAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGGGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                 <  2:129821/136‑2 (MQ=255)

                                                                                            gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAACACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTAGCtt                                                 >  2:923182/1‑140 (MQ=255)

                                                                                              tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGTAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTTCTGTTGCGTATTATACTTATTGGCTTAg                                               >  1:233175/4‑141 (MQ=255)

                                                                                              tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                               <  2:963276/138‑4 (MQ=255)

                                                                                              tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                               <  2:413987/138‑4 (MQ=255)

                                                                                              tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                               <  1:408428/138‑4 (MQ=255)

                                                                                              tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGAGACT‑TAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                               >  1:1068638/4‑141 (MQ=255)

                                                                                               ggtataATAATAGAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                              <  2:644721/139‑1 (MQ=255)

                                                                                                gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                             >  2:1128469/2‑141 (MQ=255)

                                                                                                gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                             <  1:1168156/140‑1 (MQ=255)

                                                                                                gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTCCTGTTCGGTATTATACTTATTGGCttggta                                             >  2:114599/2‑136 (MQ=255)

                                                                                                  ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATCGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATCATACTTATTGGCTTAGTaaa                                           <  1:781360/141‑3 (MQ=255)

                                                                                                  ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                           <  1:1125940/141‑3 (MQ=255)

                                                                                                  ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                           <  2:621543/141‑3 (MQ=255)

                                                                                                   taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGTTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                          <  2:691277/141‑4 (MQ=255)

                                                                                                   taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                          <  2:822873/141‑4 (MQ=255)

                                                                                                   taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                          <  2:484733/141‑4 (MQ=255)

                                                                                                       aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                      <  1:954760/139‑1 (MQ=255)

                                                                                                       aatagAGTCGCCTATCTCTCAGGCGTCAA‑‑TCGACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                      <  2:886531/139‑1 (MQ=255)

                                                                                                         tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                    <  2:448196/141‑3 (MQ=255)

                                                                                                         taaagcctccTTTCTCTCAGGCGTCAATT‑‑GATGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCTTTTACATAATAGCACCAGTCATCAGTGGCTGGGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                    <  1:30344/133‑3 (MQ=255)

                                                                                                            agTCGTCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                 <  1:637794/141‑6 (MQ=255)

                                                                                                                 ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                            >  1:21150/1‑131 (MQ=255)

                                                                                                                 ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                            >  1:1149035/1‑131 (MQ=255)

                                                                                                                 ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTAACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                            >  2:28146/1‑131 (MQ=255)

                                                                                                                 ccTATCTCTAAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                            >  2:774800/1‑131 (MQ=255)

                                                                                                                  cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAACCGCaacaaataattaa                           <  1:713053/141‑12 (MQ=255)

                                                                                                                   tATCTCTCAGGCGTGAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                          >  1:664766/1‑129 (MQ=255)

                                                                                                                   tATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAgg                           >  2:1200657/1‑141 (MQ=255)

                                                                                                                   tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                          >  1:484965/1‑129 (MQ=255)

                                                                                                                   tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                          <  1:968682/141‑13 (MQ=255)

                                                                                                                    aTCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAggg                          >  2:655623/1‑140 (MQ=255)

                                                                                                                    aTCTCTCAGGCGTCAATT‑‑GACGACGAGAGGAGGTG‑CACTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTGAGTAAACGCaacaaagaagtaata                         >  2:542516/1‑128 (MQ=255)

                                                                                                                    aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTCC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAGACGCAACACataactaata                         >  2:518825/1‑135 (MQ=255)

                                                                                                                    aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                         >  2:325084/1‑128 (MQ=255)

                                                                                                                    aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                         >  2:982572/1‑128 (MQ=255)

                                                                                                                    aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                         >  1:1056533/1‑128 (MQ=255)

                                                                                                                    aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGCTAAGCATTTTGGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaagaattaata                         >  1:619380/1‑128 (MQ=255)

                                                                                                                    aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTCTTATACTTATTGGCTTAGTAAACGCaacaaataattaata                         >  2:480161/1‑128 (MQ=255)

                                                                                                                    aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                         >  1:93053/1‑128 (MQ=255)

                                                                                                                    aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaatacttaata                         >  2:744070/1‑128 (MQ=255)

                                                                                                                    aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                         >  1:117993/1‑128 (MQ=255)

                                                                                                                          cAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACata                    >  2:793050/1‑141 (MQ=255)

                                                                                                                           aGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACatat                   >  1:257150/1‑141 (MQ=255)

                                                                                                                                  aaTTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCa            <  2:338711/141‑1 (MQ=255)

                                                                                                                                     ttAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCa         >  1:17623/1‑141 (MQ=255)

                                                                                                                                            cAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaaa  <  2:1031649/141‑1 (MQ=255)

TAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAA  >  CP000730/1929708‑1929988

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