BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,662,198	0	C	A	51.3%	20.8 / inf	224	A548E (GCA→GAA)	ynfF	S‑ and N‑oxide reductase, A subunit, periplasmic
Reads supporting (aligned to +/- strand): ref base C (60/49); new base A (61/54); total (121/103)
Fisher's exact test for biased strand distribution p-value = 7.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 224 total aligned reads displayed.

GTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTG  >  NC_000913/1661970‑1662434

gTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGa                                                                                                                                                                                                                                     >  1:1063554/1‑238 (MQ=255)

       aTTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAAc                                                                                                                                                                                                                             <  1:1641432/239‑1 (MQ=255)

        ttGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAAc                                                                                                                                                                                                                             >  2:1641432/1‑238 (MQ=255)

         tGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAAccc                                                                                                                                                                                                                           <  2:2216737/239‑1 (MQ=255)

               atcaaCATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCt                                                                                                                                                                                                                       >  2:1540072/1‑237 (MQ=255)

               atcaaCATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCt                                                                                                                                                                                                                       >  1:1548081/1‑237 (MQ=255)

               atcaaCATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGCGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCt                                                                                                                                                                                                                       >  1:824696/1‑237 (MQ=255)

                     aTGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGAt                                                                                                                                                                                                                >  2:405326/1‑238 (MQ=255)

                     aTGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGAt                                                                                                                                                                                                                >  1:1336104/1‑238 (MQ=255)

                        gCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACt                                                                                                                                                                                                                    >  1:2286370/1‑231 (MQ=255)

                        gCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACt                                                                                                                                                                                                                    <  2:2286370/231‑1 (MQ=255)

                        gCGACATCAATCACACCCATGAGGCCCTTCAGGATGACAGCAAGAGCGAGATGATTGTAGGCATTGACCACTTCATGACGGTCTCGGCTAAGAATTGCGATATACTGTTGCCCCACCTGATGCCTACAGAGCAAGAAGACATTATCTTTCATGAATCTGCAGGAAATATGGGCTATGTGATCATCGCCCAACCCGCAACCTCAGCAAAATTTGAACTCAACCCCATCTACTGGATGCtt                                                                                                                                                                                                            >  2:907793/1‑238 (MQ=255)

                               caatcaCACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAg                                                                                                                                                                                                     <  1:1057437/239‑1 (MQ=255)

                                 atcaCACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTc                                                                                                                                                                                                   >  2:1820368/1‑239 (MQ=255)

                                       cccATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCaa                                                                                                                                                                                              >  2:84345/1‑238 (MQ=255)

                                           tGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTgg                                                                                                                                                                                                                  <  2:449281/214‑1 (MQ=255)

                                           tGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTgg                                                                                                                                                                                                                  >  1:449281/1‑214 (MQ=255)

                                                tGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATg                                                                                                                                                                                                               <  1:1551530/212‑1 (MQ=255)

                                                  cTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCaa                                                                                                                                                                                                                                            <  1:1046974/181‑1 (MQ=255)

                                                  cTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCaaa                                                                                                                                                                                                                              >  2:475553/1‑195 (MQ=255)

                                                  cTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCaaa                                                                                                                                                                                                                              <  1:475553/195‑1 (MQ=255)

                                                                 aaGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTAt                                                                                                                                                                        <  1:2359075/234‑1 (MQ=255)

                                                                  aGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCaaa                                                                                                                                                                    >  2:1560243/1‑237 (MQ=255)

                                                                                  tgGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCa                                                                                                                                                  <  1:1263341/239‑1 (MQ=255)

                                                                                        tGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGc                                                                                                                                            <  2:316316/239‑1 (MQ=255)

                                                                                         gACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAAccc                                                                                                                                                                                                                           <  2:1985365/159‑1 (MQ=255)

                                                                                         gACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAAccc                                                                                                                                                                                                                           >  1:213913/1‑159 (MQ=255)

                                                                                             aCTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGaa                                                                                                                                       <  1:398106/239‑1 (MQ=255)

                                                                                                 cgtGACGGCATCGGCTATGTACTGCGATATCCTGATGCCCTACCTAATTCCAACAGAGCAAGACGACCTTATCACTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGGAACCTCAGCAAAATTTGAACGCAAACCTATCTACTGGATGCTGGGTGAAGTCGCCAAACGCTTAGGTCCAGACGTTTATCAAACCTTTACTGAAGGTCGCAATCAGCATGAATGGa                                                                                                                                   <  2:948175/237‑1 (MQ=255)

                                                                                                         ccTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCaaa                                                                                                                                                                                                                              <  2:1602284/140‑1 (MQ=255)

                                                                                                         ccTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCa                                                                                                                                                                                                                                >  2:1556572/1‑138 (MQ=255)

                                                                                                            cGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCt                                                                                                                                                                                         <  1:1080343/174‑1 (MQ=255)

                                                                                                            cGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATAtct                                                                                                                         >  1:361263/1‑238 (MQ=255)

                                                                                                                       tGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCa                                                                                                                                                                                                                                <  1:716069/124‑1 (MQ=255)

                                                                                                                       tGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCa                                                                                                                                                                                                                                >  2:716069/1‑124 (MQ=255)

                                                                                                                       tGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAAttt                                                                                                                                                                                                                                       >  1:530292/1‑117 (MQ=255)

                                                                                                                       tGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAAttt                                                                                                                                                                                                                                       <  2:530292/117‑1 (MQ=255)

                                                                                                                            tatCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAAccc                                                                                                                                                                                                                           <  2:236529/124‑1 (MQ=255)

                                                                                                                            tatCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCaa                                                                                                                                                                                              <  2:2278866/153‑1 (MQ=255)

                                                                                                                            tatCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAAc                                                                                                                                                                   <  2:1466549/180‑1 (MQ=255)

                                                                                                                              tCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGa                                                                                                       >  1:1829025/1‑238 (MQ=255)

                                                                                                                              tCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAAtt                                                                                                                                                                                                                                        <  2:281241/109‑1 (MQ=255)

                                                                                                                                cTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACg                                                                                                    <  2:105507/239‑1 (MQ=255)

                                                                                                                                 tGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGt                                                                                                   <  1:611246/239‑1 (MQ=255)

                                                                                                                                   ttGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGAc                                                                                                                                                                                   >  1:1588649/1‑157 (MQ=255)

                                                                                                                                   ttGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAAttt                                                                                                                                                                                                                                       <  1:1154887/105‑1 (MQ=255)

                                                                                                                                   ttGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCa                                                                                                                                                                                 <  2:380039/159‑1 (MQ=255)

                                                                                                                                   ttGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCa                                                                                                                                                                                 >  1:380039/1‑159 (MQ=255)

                                                                                                                                   ttGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTa                                                                                                  >  1:1489210/1‑238 (MQ=255)

                                                                                                                                   ttGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCCCCATGCGAAAACGAAGGAACGTa                                                                                                  >  2:867823/1‑238 (MQ=255)

                                                                                                                                      cccGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAAcc                                                                                                                                                                                                                            <  2:755937/113‑1 (MQ=255)

                                                                                                                                        cGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAAc                                                                                                >  2:2031377/1‑235 (MQ=255)

                                                                                                                                          aCCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGcc                                                                                                                                                                                                <  1:1715597/137‑1 (MQ=255)

                                                                                                                                                tGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGcc                                                                                     <  2:1336104/238‑1 (MQ=255)

                                                                                                                                                     aCAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGaa                                                                                                                                                                                                                                    >  1:997384/1‑90 (MQ=255)

                                                                                                                                                     aCAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACt                                                                                >  1:1724095/1‑238 (MQ=255)

                                                                                                                                                     aCAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACt                                                                                >  1:2171806/1‑238 (MQ=255)

                                                                                                                                                         agCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGAcc                                                                                                                                                                                  >  2:1058016/1‑136 (MQ=255)

                                                                                                                                                         agCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGt                                                                                                                                                <  1:1463687/170‑1 (MQ=255)

                                                                                                                                                           cAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTg                                                                                                                                                                                                        <  2:1442256/112‑1 (MQ=255)

                                                                                                                                                                agaCCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACg                                                                             <  1:1364515/230‑1 (MQ=255)

                                                                                                                                                                   ccTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGaaa                                                                  >  1:373565/1‑238 (MQ=255)

                                                                                                                                                                      tATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCtt                                                                                                                                                                                        >  2:1959384/1‑117 (MQ=255)

                                                                                                                                                                      tATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCtt                                                                                                                                                                                        <  1:1959384/117‑1 (MQ=255)

                                                                                                                                                                       aTCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCa                                                                                                                                                                                                                                >  2:70591/1‑76 (MQ=255)

                                                                                                                                                                       aTCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCa                                                                                                                                                                                                                                <  1:70591/76‑1 (MQ=255)

                                                                                                                                                                           ctcATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCggg                                                         <  1:262746/239‑1 (MQ=255)

                                                                                                                                                                            tcATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACg                                                                                                                                                                             <  1:567410/122‑1 (MQ=255)

                                                                                                                                                                            tcATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACg                                                                                                                                                                             >  2:567410/1‑122 (MQ=255)

                                                                                                                                                                             cATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCaa                                                                                                                                                                     >  2:2177702/1‑129 (MQ=255)

                                                                                                                                                                                   tCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCAt                                                                                                                                                                                                                         <  2:1825959/71‑1 (MQ=255)

                                                                                                                                                                                            aaTATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGa                                                                                                                                                                                                       <  2:895520/80‑1 (MQ=255)

                                                                                                                                                                                            aaTATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCggg                                                         >  1:300478/1‑222 (MQ=255)

                                                                                                                                                                                            aaTATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTAAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCggg                                                         <  2:300478/222‑1 (MQ=255)

                                                                                                                                                                                                 gggCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAAt                                                                                                                                      <  1:1438845/140‑1 (MQ=255)

                                                                                                                                                                                                 gggCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCgg                                   <  2:1136542/239‑1 (MQ=255)

                                                                                                                                                                                                 gggCTATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCgg                                   >  1:1908966/1‑239 (MQ=255)

                                                                                                                                                                                                     tATGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACg                                                                                                    <  1:2040745/170‑1 (MQ=255)

                                                                                                                                                                                                      aTGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCa                                                                                                                                                                                                                                >  1:2196506/1‑45 (MQ=255)

                                                                                                                                                                                                      aTGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGa                                                                                                                                   <  2:1762885/138‑1 (MQ=255)

                                                                                                                                                                                                      aTGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGa                                                                                                                                   >  1:1762885/1‑138 (MQ=255)

                                                                                                                                                                                                      aTGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAg                                                                                                         >  2:2116396/1‑164 (MQ=255)

                                                                                                                                                                                                      aTGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAg                                                                                                         <  1:2116396/164‑1 (MQ=255)

                                                                                                                                                                                                      aTGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAAc                                                                                                <  1:893575/173‑1 (MQ=255)

                                                                                                                                                                                                      aTGTGATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCa                                                                                                                                                                                                                                >  1:1535948/1‑45 (MQ=255)

                                                                                                                                                                                                          gATCCTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACt                          <  2:1514522/239‑1 (MQ=255)

                                                                                                                                                                                                            tCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTgg                                                                                                                                                                                                                  >  1:797839/1‑53 (MQ=255)

                                                                                                                                                                                                             ccTCGCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATg                                                                                                                                     >  2:1374802/1‑129 (MQ=255)

                                                                                                                                                                                                              cTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATg                                                                                                                                                                                                               >  1:852275/1‑54 (MQ=255)

                                                                                                                                                                                                               tCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTc                     >  1:752233/1‑239 (MQ=255)

                                                                                                                                                                                                                 gCCCAACCCGCAACCTCAGCAAAATTTGAACGCaaa                                                                                                                                                                                                                              >  2:1717465/1‑36 (MQ=255)

                                                                                                                                                                                                                 gCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACCGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGCAGGTCGCCGTCAGCATGAATGGATCAAATATCTCAATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGCAGCACACTAAGTc                     >  2:1236658/1‑237 (MQ=255)

                                                                                                                                                                                                                 gCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAAcc                                                                                                                                                                                                                            >  1:946595/1‑38 (MQ=255)

                                                                                                                                                                                                                 gCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTa                                                                                                                                                                                       >  2:1076468/1‑75 (MQ=255)

                                                                                                                                                                                                                 gCCCAACCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGAt                                                                                                                                  >  1:1557051/1‑128 (MQ=255)

                                                                                                                                                                                                                      aCCCGCAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGa                                                                                                                                   <  2:584990/122‑1 (MQ=255)

                                                                                                                                                                                                                           cAACCTCAGAAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCa         >  2:188943/1‑239 (MQ=255)

                                                                                                                                                                                                                                 cAGAAAAATTTGAACGCAAACCCATCTACTGGATg                                                                                                                                                                                                               >  2:813073/1‑35 (MQ=255)

                                                                                                                                                                                                                                  aGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTa                                                                                                                                                              >  1:1823483/1‑83 (MQ=255)

                                                                                                                                                                                                                                  aGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTg  >  1:1097989/1‑239 (MQ=255)

GTAACACATTGATCAATCAACATGGCGACATCAATCACACCCATGAGGTGCTTCAGGATGACAGCAAGTGCGAGATGATTGTTGGCATTGACCACTTCATGACGGCCTCGGCTAAGTATTGCGATATCCTGTTGCCCGACCTGATGCCAACAGAGCAAGAAGACCTTATCTCTCATGAATCTGCAGGGAATATGGGCTATGTGATCCTCGCCCAACCCGCAACCTCAGCAAAATTTGAACGCAAACCCATCTACTGGATGCTGAGTGAAGTCGCCAAACGCTTAGGACCAGACGTTTATCAAACCTTTACTGAAGGTCGCAGTCAGCATGAATGGATCAAATATCTCCATGCGAAAACGAAGGAACGTAACCCTGAGATGCCCGACTACGAGGAGATGAAAACGACCGGGATCTTTAAGAAAAAATGCCCGGAAGAACACTACGTCGCTTTCCGCGCATTCCGTG  >  NC_000913/1661970‑1662434

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