BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	2,455,698	2 bp→AT	coding (1216‑1217/1260 nt)	USA300HOU_2315 →	hypothetical membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,455,698	0	C	A	89.3%	389.0 / 41.9	134	coding (1216/1260 nt)	USA300HOU_2315	hypothetical membrane protein
Reads supporting (aligned to +/- strand): ref base C (0/2); major base A (45/72); minor base . (6/8); total (52/82)
Fisher's exact test for biased strand distribution p-value = 7.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	CP000730	2,455,699	0	A	T	97.1%	397.8 / ‑2.7	139	H406L (CAT→CTT)	USA300HOU_2315	hypothetical membrane protein
Reads supporting (aligned to +/- strand): ref base A (1/3); new base T (51/82); total (52/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

AGCAGAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCCTCTTGAGTG  >  CP000730/2455559‑2455841

||

agcagAAGCAGTTAAACAACCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGaaa                                                                                                                                        <  2:1092155/148‑1 (MQ=255)

    gAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAAc                                                                                                                                    >  1:574101/1‑148 (MQ=255)

     aaTCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACt                                                                                                                                   <  1:300845/148‑1 (MQ=255)

           taaaaCAAGTACCACAAAGGGACGAATTAAAAAAAGATACAAATGAAGATGAAATAGAGGAATCATTAAAAGAACCTGTGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAAAAAAAAAAAAAACTgaagaa                                                                                                                             <  2:742453/146‑1 (MQ=255)

            tAAACAAGCACCACAAAGAGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAACGAACCTGCGACTTATAAAGAAGACAGACGTCACACAGTAGTAATTGATGCAAAAAAA‑TCATCGACAAAACTGAAGAAc                                                                                                                            >  2:551088/1‑148 (MQ=255)

                cAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCACTCATATAAATAAACTGCGGTGCGCAAATAAGACAGACGTTCATCAGTTGTAAGTTATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATc                                                                                                                        >  1:1031045/1‑148 (MQ=255)

                cAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAA‑TTATCGAAAAAACTGAAGAAAAATc                                                                                                                        >  1:597604/1‑148 (MQ=255)

                cAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTgaagaaacaac                                                                                                                        >  1:711429/1‑143 (MQ=255)

                cAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATc                                                                                                                        >  1:1354204/1‑148 (MQ=255)

                cAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATc                                                                                                                        >  1:1109203/1‑148 (MQ=255)

                  aGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAA‑TTATCGAAAAAACTGAAGAACAATCtt                                                                                                                      <  2:1193660/148‑1 (MQ=255)

                  aGCACCACAAAGTGACCATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAA‑TTATCGAAAAAACTGAAGAACAATCtt                                                                                                                      <  2:143317/148‑1 (MQ=255)

                    caccacAAAGTGACGATTTAACAAACGATTCAAATGACGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATCATCGAAAAAACTGAAGAACAATCTTCa                                                                                                                    <  1:916886/148‑1 (MQ=255)

                    caccacAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCa                                                                                                                    >  1:1057781/1‑148 (MQ=255)

                     actacaAAGTGACGATTTAACAAGCGATTAAAAAGAAGAAGAATCAGGGAAACCATTAAAATAACCTTCGACTTATAAAGAAGACAGACGTTCAATAGTTGTAATTGATGCAAAAAAA‑TTATCGACAAAACTTAAGAACAATCTTCAg                                                                                                                   <  2:992316/145‑1 (MQ=255)

                     accacAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAAGCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAg                                                                                                                   <  2:776136/148‑1 (MQ=255)

                      ccacAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGa                                                                                                                  <  2:212281/148‑1 (MQ=255)

                      ccacAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGa                                                                                                                  <  2:717151/148‑1 (MQ=255)

                       cacaAGGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCatat                                                                                                                 <  2:94089/148‑1 (MQ=255)

                       cacaAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGAt                                                                                                                 <  2:952090/148‑1 (MQ=255)

                        acaAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATa                                                                                                                <  2:1228290/148‑1 (MQ=255)

                         caAAGTGACGATTTAACAAACGATTCAAATCAAGAAGAAATAGGGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaa                                                                                                               <  2:768872/148‑1 (MQ=255)

                           aaGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCAGTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaaaa                                                                                                             <  2:57385/148‑1 (MQ=255)

                           aaGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATAATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATaaaa                                                                                                             >  2:882457/1‑148 (MQ=255)

                               gACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCAGCAGTTGTAATTGATGCAAAAAAA‑TTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACaa                                                                                                         >  1:737653/1‑148 (MQ=255)

                                aCGATTTAACAGACGATTCAAATGAAGATGAAAAAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTGCATCAGTTGTAAAAAAAAAA‑AAAAAATTATCGAAAAAACGGAAGAACAATCTTCAGATAAAAACaaa                                                                                                        <  2:788551/148‑1 (MQ=255)

                                 cGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaat                                                                                                       >  1:801241/1‑148 (MQ=255)

                                 cGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaat                                                                                                       >  1:256921/1‑148 (MQ=255)

                                   aTTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaataa                                                                                                     <  2:1169366/148‑1 (MQ=255)

                                    tttAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaataat                                                                                                    <  2:382573/148‑1 (MQ=255)

                                    tttAACAAACGATTCAAAAGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAACGAAGACAGACGTTCATCATTTGTAATGGATGAAAAAAAAC‑CAACGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaataat                                                                                                    <  1:737578/148‑1 (MQ=255)

                                    tttAACAAAAGATTCAAATCAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAaataat                                                                                                    <  2:1099431/148‑1 (MQ=255)

                                     ttAACCAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAataata                                                                                                   <  1:1072474/148‑1 (MQ=255)

                                              gATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAgcgc                                                                                          >  2:1049915/1‑148 (MQ=255)

                                                 tCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTg                                                                                       <  2:835086/148‑1 (MQ=255)

                                                  cAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGc                                                                                      <  1:263358/148‑1 (MQ=255)

                                                      tGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACCTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGt                                                                                  <  1:1190908/148‑1 (MQ=255)

                                                      tGAAGATGAAAGAGAGCAATCATTAACAGAACCTGCAACATATAAAGAAGACAGACGTTTATCAGTTGAAAAAAAAAAA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGt                                                                                  <  2:1398682/148‑1 (MQ=255)

                                                        aaGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTAc                                                                                >  2:481991/1‑148 (MQ=255)

                                                            ttaaaGAGAGCAATCAGTAAAACAACCTGCGACCTATAAACAAGACAGACGTTCACCAGTTGCAATTGATGCAAAAAAA‑AAAAAAAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACctct                                                                            <  2:347472/146‑1 (MQ=255)

                                                            tGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAA‑TTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACctct                                                                            <  1:1146270/148‑1 (MQ=255)

                                                             gAAATAGAGCAAACATTAAAAGGACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTTAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCtt                                                                           <  2:912988/148‑1 (MQ=255)

                                                             gAAATAGAGCAAACATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCtt                                                                           <  2:1090833/148‑1 (MQ=255)

                                                              aaaTAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTa                                                                          >  1:868400/1‑148 (MQ=255)

                                                               aagagagCAACCATTAAAAGGACCTGGGACTTATAAAGAAGACGGACGGTGGGGGGTGGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaa                                                                         <  1:960175/145‑1 (MQ=255)

                                                               aaTAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAACAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaa                                                                         <  1:937599/148‑1 (MQ=255)

                                                               aaTAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTAGCTAGTACCTCTTaa                                                                         >  1:111143/1‑148 (MQ=255)

                                                                  agagCACTCATTAAAAGCACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaaaaa                                                                      <  1:316006/148‑1 (MQ=255)

                                                                   gagCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTAAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAg                                                                     >  1:44455/1‑148 (MQ=255)

                                                                      cAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAA‑TTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATACAAAAAAAGCGCTTCGTAGTCACTCTTAAAAAgaga                                                                  >  2:348165/1‑146 (MQ=255)

                                                                             tAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAg                                                           >  2:1260200/1‑148 (MQ=255)

                                                                             tAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGAAAATGCTAg                                                           >  2:1347094/1‑148 (MQ=255)

                                                                              aaaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAA‑TTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCCCTTAAAAAGATGATGCTAGc                                                          >  2:407455/1‑148 (MQ=255)

                                                                               aaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAAATGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATAATGCTAgaa                                                         >  2:660503/1‑146 (MQ=255)

                                                                                 aGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAg                                                       >  2:797933/1‑148 (MQ=255)

                                                                                  gAACCCGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgc                                                      <  1:4906/148‑1 (MQ=255)

                                                                                   aaCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAGTAACAAAGCGATCGATAGTACCTCTTAAAAAGATTATGCTAGCAAgcg                                                     <  1:1037281/148‑1 (MQ=255)

                                                                                    cccTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATACCAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgcgc                                                    <  1:198476/147‑1 (MQ=255)

                                                                                    aCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgcgc                                                    <  2:334587/148‑1 (MQ=255)

                                                                                    aCCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAgcgc                                                    <  2:980972/148‑1 (MQ=255)

                                                                                     ccTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCt                                                   <  2:617328/148‑1 (MQ=255)

                                                                                      cTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAA‑TTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCtt                                                  <  1:512127/148‑1 (MQ=255)

                                                                                      cTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCtt                                                  <  1:21793/148‑1 (MQ=255)

                                                                                        gCGACTTATAAAGAAGAAAGACGTGCATCAGGTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCGTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCtttt                                                <  2:337103/148‑1 (MQ=255)

                                                                                         caaCTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTc                                               <  2:1062018/146‑1 (MQ=255)

                                                                                         cGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTc                                               <  1:969764/148‑1 (MQ=255)

                                                                                         cGAATTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTc                                               <  1:421168/148‑1 (MQ=255)

                                                                                          gACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTct                                              <  1:1347512/148‑1 (MQ=255)

                                                                                          gACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATGATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTct                                              <  1:1316950/148‑1 (MQ=255)

                                                                                           aCTTAAAAACAAGACAGACGTTCAGCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACATTTAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTcta                                             <  1:41646/148‑1 (MQ=255)

                                                                                            cTTATACAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctat                                            <  2:903934/148‑1 (MQ=255)

                                                                                            cTTATAAAGAAGACCGACGTTTATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATCACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctat                                            <  1:46801/148‑1 (MQ=255)

                                                                                              tataAAGAAGACAGATGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatac                                          <  1:1233383/148‑1 (MQ=255)

                                                                                               ataAAGAAGACAGACCTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatact                                         <  2:604551/148‑1 (MQ=255)

                                                                                                taAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatacta                                        >  1:811401/1‑148 (MQ=255)

                                                                                                taAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatacta                                        >  2:950846/1‑148 (MQ=255)

                                                                                                gaaaGAAGACAGACGTTCATCAGTGGTAATTGATTCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTctatacta                                        <  2:204818/147‑1 (MQ=255)

                                                                                                   agaagaCAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACtagat                                     >  1:97035/1‑148 (MQ=255)

                                                                                                      agacagacGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATAtatt                                  >  1:124600/1‑148 (MQ=255)

                                                                                                         cagacGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAGCTGAAGAACAATATTCAGATAAAAACAAATAATAACCAAGCGCCTGCTAGCACCTCTTAACAAGAAGATCCTACCCAGCGCTTTTCCATCCTACACattttt                               >  2:1264515/1‑148 (MQ=255)

                                                                                                          agacGTTCATCAGTTGTAATTGATGCA‑AAAAAATGATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTAttt                              >  2:873506/1‑148 (MQ=255)

                                                                                                           gacTTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGCTAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTAtttt                             <  1:949920/148‑1 (MQ=255)

                                                                                                              gTTCATCAGTTGTAATTGATGCAAAAAAA‑TTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCGTAAAAAGATGATGCTAGCGAGCGCTTTTCTATACTATATATTATTTTTCt                          >  1:719346/1‑148 (MQ=255)

                                                                                                               gtcatcaGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTa                         <  2:1313297/147‑1 (MQ=255)

                                                                                                                  atcaGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTaaaa                      <  2:256921/148‑1 (MQ=255)

                                                                                                                  atcaGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTaaaa                      <  2:801241/148‑1 (MQ=255)

                                                                                                                   tcaGTTGTAATTGATGCA‑AAAAAATTATCGAANAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAAt                     <  1:159242/148‑1 (MQ=255)

                                                                                                                   aaaGTTGTAATTAATGCAAAAAAAC‑TATCGAAAAAACTAAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATTATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAAt                     <  1:999665/146‑1 (MQ=255)

                                                                                                                     aGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAAttt                   <  1:539900/148‑1 (MQ=255)

                                                                                                                       ttGTAATTGATG‑AAAAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTa                 <  1:548658/148‑1 (MQ=255)

                                                                                                                        tGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTaa                >  1:435714/1‑148 (MQ=255)

                                                                                                                          tAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCACGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAaca              >  2:1380389/1‑148 (MQ=255)

                                                                                                                          tAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATAATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTCAAATTGTAcca              >  1:268498/1‑148 (MQ=255)

                                                                                                                             ttGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACAccc           >  1:1295268/1‑148 (MQ=255)

                                                                                                                             ttGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACAccc           >  1:477577/1‑148 (MQ=255)

                                                                                                                                 atgcaaaaaaaTTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCctct        <  1:1152721/143‑1 (MQ=255)

                                                                                                                                 atgcaaaaaaaTTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCctct        <  1:453538/143‑1 (MQ=255)

                                                                                                                                      aaaaaaTTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCCTCTTGAGt   <  1:1273214/148‑1 (MQ=255)

                                                                                                                                      aaaaaaTTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCCTCTTGAGt   <  1:809975/148‑1 (MQ=255)

                                                                                                                                      aaaaaaTTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCCTCTTGAGt   <  1:1147573/148‑1 (MQ=255)

||

AGCAGAATCAGTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTTAACACCCTCTTGAGTG  >  CP000730/2455559‑2455841

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