BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,456,723	A→C	8.6%	D78A (GAT→GCT)	yfdF →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,456,723	0	A	C	8.6%	889.5 / 68.8	307	D78A (GAT→GCT)	yfdF	uncharacterized protein
Reads supporting (aligned to +/- strand): ref base A (148/132); new base C (11/16); total (159/148)
Fisher's exact test for biased strand distribution p-value = 3.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 328 total aligned reads displayed.

GTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTGAGAAGACTCATGATATCAATAACAAGG  >  CP009273/2456578‑2456872

gTACAAGAGTTCAAAAACCTTTTTAATGGTGAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTatattgctatt                                                                                                                                                   >  2:4459691/1‑150 (MQ=255)

gTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTatattgctatt                                                                                                                                                   >  2:1299643/1‑150 (MQ=255)

 tACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTatattgctattg                                                                                                                                                  <  2:81335/150‑1 (MQ=255)

     agagTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTc                                                                                                                                              >  2:1550617/1‑150 (MQ=255)

     agagTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTc                                                                                                                                              >  1:1901927/1‑150 (MQ=255)

     agagTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTc                                                                                                                                              >  1:3820794/1‑150 (MQ=255)

      tagTTCAAAAACCTTTTTAATGGTAAGGAAGTAATAATCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTcc                                                                                                                                             <  2:2691517/149‑1 (MQ=255)

            aaaaaCCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGAt                                                                                                                                       >  1:1684482/1‑150 (MQ=255)

             aaaaCCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATa                                                                                                                                      >  1:2238509/1‑150 (MQ=255)

              aaaCCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATAc                                                                                                                                     >  1:4052995/1‑150 (MQ=255)

              aaaCCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATAc                                                                                                                                     >  1:1800778/1‑150 (MQ=255)

                 ccTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGCTATTGACTCCAATGATACAGa                                                                                                                                  <  2:3157940/150‑1 (MQ=255)

                 ccTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGCTATTGACTCCAATGATACAGa                                                                                                                                  <  2:892624/150‑1 (MQ=255)

                 ccTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGa                                                                                                                                  >  2:4188581/1‑150 (MQ=255)

                  cattttAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGAc                                                                                                                                 <  2:690491/148‑1 (MQ=255)

                     tttAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCat                                                                                                                              <  2:84628/150‑1 (MQ=255)

                        aaTGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATAtt                                                                                                                           >  1:3154842/1‑150 (MQ=255)

                            gTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTa                                                                                                                       >  1:1413747/1‑150 (MQ=255)

                             tAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATAGTGACTCCAATGATACAGACCATATTACTAt                                                                                                                      >  2:3067882/1‑150 (MQ=255)

                              aaggaaggAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTAtt                                                                                                                     >  2:2783818/1‑150 (MQ=255)

                                 gaaggaaTAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGt                                                                                                                  <  2:3811743/150‑1 (MQ=255)

                                  aaggaaTAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGCTATTGACTCCAATGATACAGACCATATTACTATTGGTa                                                                                                                 >  2:1890603/1‑150 (MQ=255)

                                  aaggaaTAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATAAGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTa                                                                                                                 <  1:2638304/150‑1 (MQ=255)

                                      aaTAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGa                                                                                                             <  1:4531868/150‑1 (MQ=255)

                                        tAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACa                                                                                                           >  2:2691845/1‑150 (MQ=255)

                                        tAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACa                                                                                                           <  2:1735117/150‑1 (MQ=255)

                                         aaGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGCTATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACaa                                                                                                          <  1:1878885/150‑1 (MQ=255)

                                          aGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAAc                                                                                                         <  1:1625931/150‑1 (MQ=255)

                                           gCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAAcc                                                                                                        >  2:3434234/1‑150 (MQ=255)

                                              ccTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAg                                                                                                     <  2:3493119/150‑1 (MQ=255)

                                               cTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGa                                                                                                    >  1:318400/1‑150 (MQ=255)

                                                tGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGac                                                                                                   >  1:4412613/1‑149 (MQ=255)

                                                tGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGaa                                                                                                   <  1:795812/150‑1 (MQ=255)

                                                tGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGaa                                                                                                   >  1:2678703/1‑150 (MQ=255)

                                                   aTTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCa                                                                                                <  1:1660065/150‑1 (MQ=255)

                                                   aTTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCa                                                                                                <  2:635263/150‑1 (MQ=255)

                                                   aTTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCa                                                                                                >  2:3207803/1‑150 (MQ=255)

                                                     tAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAAt                                                                                              >  1:3432506/1‑150 (MQ=255)

                                                         cATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGaa                                                                                          <  1:4544678/150‑1 (MQ=255)

                                                            cTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTACTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTa                                                                                       >  2:4377123/1‑150 (MQ=255)

                                                              aCTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTAcc                                                                                     <  2:36210/150‑1 (MQ=255)

                                                              aCTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTAcc                                                                                     >  1:2959725/1‑150 (MQ=255)

                                                                   aGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGa                                                                                >  2:355824/1‑150 (MQ=255)

                                                                     cTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAAc                                                                              >  1:1583379/1‑150 (MQ=255)

                                                                       tataAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACta                                                                            <  1:277098/150‑1 (MQ=255)

                                                                            aaaaCGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAg                                                                       >  1:1228578/1‑150 (MQ=255)

                                                                            aaaaCGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAg                                                                       <  1:1252929/150‑1 (MQ=255)

                                                                                cGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACaa                                                                   >  1:1035366/1‑150 (MQ=255)

                                                                                  ccATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGCTATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAAACAATAGAATTACCTGCGAACTATAAAGACAAAg                                                                 >  2:1324451/1‑150 (MQ=255)

                                                                                   cATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACATGCGAACTATAAAGACAAAGa                                                                >  2:2640412/1‑150 (MQ=255)

                                                                                    aTACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGaa                                                               <  2:3442587/150‑1 (MQ=255)

                                                                                    aTACGAGTAAACGACGATCCTTAAAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGaa                                                               >  1:3752971/1‑150 (MQ=255)

                                                                                      aCGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACt                                                             <  2:2590665/150‑1 (MQ=255)

                                                                                      aCGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACt                                                             <  2:1540737/150‑1 (MQ=255)

                                                                                      aCGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACt                                                             >  2:2081262/1‑150 (MQ=255)

                                                                                       cGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTc                                                            <  1:1975768/150‑1 (MQ=255)

                                                                                        gAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCg                                                           <  1:1206478/150‑1 (MQ=255)

                                                                                          gTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTc                                                         <  2:1970397/150‑1 (MQ=255)

                                                                                          gTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTc                                                         <  2:317977/150‑1 (MQ=255)

                                                                                             aaCGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGt                                                      <  2:103876/150‑1 (MQ=255)

                                                                                                gacgaTCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACt                                                   <  2:236188/150‑1 (MQ=255)

                                                                                                   gaTCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATc                                                >  2:4069323/1‑150 (MQ=255)

                                                                                                    aTCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCa                                               >  1:4463413/1‑150 (MQ=255)

                                                                                                     tCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCAt                                              <  2:428078/150‑1 (MQ=255)

                                                                                                     tCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCAt                                              <  2:4153833/150‑1 (MQ=255)

                                                                                                      ccTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCAtt                                             <  2:480722/150‑1 (MQ=255)

                                                                                                         tataGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAAt                                          <  1:2581235/150‑1 (MQ=255)

                                                                                                          ataGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATg                                         >  1:3048981/1‑150 (MQ=255)

                                                                                                            aGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACGAAGAACCAATAGAATTACCTGCGAACTATCAAGACAAAGAACTCGTCCGTACTATCATTAATGac                                       >  2:2843311/1‑150 (MQ=255)

                                                                                                            aGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGac                                       >  2:2177757/1‑150 (MQ=255)

                                                                                                              aTTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGacaa                                     <  2:4475770/150‑1 (MQ=255)

                                                                                                              aTTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGacaa                                     <  2:1478476/150‑1 (MQ=255)

                                                                                                                  aatattaataATTCCTCAGTTACTTATATTGCTATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATa                                 <  1:762315/150‑1 (MQ=255)

                                                                                                                  aatattaataATTCCTCAGTTACTTATATTGCTATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATa                                 >  2:891578/1‑150 (MQ=255)

                                                                                                                  aatattaataATTCCTCAGTTACTTATATTGCTATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATa                                 >  2:806352/1‑150 (MQ=255)

                                                                                                                   atattaataATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAg                                <  1:3003281/150‑1 (MQ=255)

                                                                                                                   atattaataATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAg                                >  2:1427679/1‑150 (MQ=255)

                                                                                                                   atattaataATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAg                                <  2:3749778/150‑1 (MQ=255)

                                                                                                                   atattaataATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAg                                >  2:3942128/1‑150 (MQ=255)

                                                                                                                    tattaataATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGAGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGt                               >  2:2258534/1‑150 (MQ=255)

                                                                                                                     attaataATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGtt                              <  2:2451696/150‑1 (MQ=255)

                                                                                                                        aataatTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTgag                           >  1:4282926/1‑150 (MQ=255)

                                                                                                                        aataatTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTgag                           >  1:1480884/1‑150 (MQ=255)

                                                                                                                         ataatTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTgaga                          >  1:1312015/1‑150 (MQ=255)

                                                                                                                         ataatTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTgaga                          >  1:1313563/1‑150 (MQ=255)

                                                                                                                          taatTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTGagaa                         <  1:2376715/150‑1 (MQ=255)

                                                                                                                           aatTCCTCAGTTACTTATATTGCTATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACAAAGTTGagaag                        >  1:512466/1‑150 (MQ=255)

                                                                                                                           aatTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTGagaag                        >  1:3104532/1‑150 (MQ=255)

                                                                                                                             ttCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTGAGAAGAc                      <  2:2425980/150‑1 (MQ=255)

                                                                                                                             ttCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTGAGAAGAc                      <  2:123715/150‑1 (MQ=255)

                                                                                                                                cTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACCAAGAACTCGTCCGTACTATCATTAAGGACCACATAGTTGAGAAGACTCa                   >  2:3062580/1‑150 (MQ=255)

                                                                                                                                cTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTGAGAAGACTCa                   <  1:330693/150‑1 (MQ=255)

                                                                                                                                 tCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTGAGAAGACTCAt                  <  2:1950313/150‑1 (MQ=255)

                                                                                                                                       aCTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTGAGAAGACTCATGATATc            >  2:336085/1‑150 (MQ=255)

                                                                                                                                        cTTATATTGCTATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTGAGAAGACTCATGATATCa           <  1:3574460/150‑1 (MQ=255)

                                                                                                                                        cTTATATTGCTATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTGAGAAGACTCATGATATCa           <  1:4516235/150‑1 (MQ=255)

                                                                                                                                         ttATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTGAGAAGACTCATGATATCaa          >  1:3359429/1‑150 (MQ=255)

                                                                                                                                          tataTTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTGAGAAGACTCATGATATCAAt         <  2:1081283/150‑1 (MQ=255)

                                                                                                                                                gatattgaCTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTGAGAAGACTCATGATATCAATAACAAg   <  1:3540347/150‑1 (MQ=255)

                                                                                                                                                 atattgaCTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTGAGAAGACTCATGATATCAATAACAAgg  >  2:4188920/1‑150 (MQ=255)

GTACAAGAGTTCAAAAACCTTTTTAATGGTAAGGAAGGAATAAGCACCTGTATTAAACATCTACTTGAGCTTATAAAAAACGCCATACGAGTAAACGACGATCCTTATAGATTTAATATTAATAATTCCTCAGTTACTTATATTGATATTGACTCCAATGATACAGACCATATTACTATTGGTATCGACAACCAAGAACCAATAGAATTACCTGCGAACTATAAAGACAAAGAACTCGTCCGTACTATCATTAATGACAACATAGTTGAGAAGACTCATGATATCAATAACAAGG  >  CP009273/2456578‑2456872

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