BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,411,241	C→A	5.4%	A38S (GCA→TCA)	yfbS ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,411,241	0	C	A	5.4%	3455.9 / 23.2	1027	A38S (GCA→TCA)	yfbS	putative transporter
Reads supporting (aligned to +/- strand): ref base C (616/317); new base A (3/91); total (619/408)
Fisher's exact test for biased strand distribution p-value = 1.49e-35
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1061 total aligned reads displayed.

CGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCT  >  NC_000913/2411092‑2411390

cGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGc                                                                                                                                                       <  1:2542767/150‑1 (MQ=255)

 gCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCa                                                                                                                                                      >  1:8632984/1‑150 (MQ=255)

    aCTTTTACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAAc                                                                                                                                                   <  1:4506051/150‑1 (MQ=255)

      tttGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAAcgc                                                                                                                                                 >  2:2082728/1‑150 (MQ=255)

      tttGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAAcgc                                                                                                                                                 >  2:8749008/1‑150 (MQ=255)

      tttGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAAcgc                                                                                                                                                 <  2:3950321/150‑1 (MQ=255)

        tGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGa                                                                                                                                               <  2:2763779/150‑1 (MQ=255)

        tGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGa                                                                                                                                               <  2:8597643/150‑1 (MQ=255)

              gCCATGTTCCCATTCCGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATaa                                                                                                                                         <  1:3598874/150‑1 (MQ=255)

              gCCATGTTCCCATTACGGTGGCAACACCGGTACGTACCAAACCATCGCAAATAATAAACAAGGCGGCAATCAGGACAACTTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATaa                                                                                                                                         <  2:7492161/150‑1 (MQ=255)

              gCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATaa                                                                                                                                         >  1:11140865/1‑150 (MQ=255)

                  ttttCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATaacaaa                                                                                                                                     <  1:4652818/148‑1 (MQ=255)

                  tGTTCGCATTTCGGTGGCTACAACTGTACTTACAAAACCATCGCCAATAATAAACAAGGCGTCAATCTGGACAACGTTAGGATCATAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCAGCTTAATGCAAACGCGACAATaacaaa                                                                                                                                     <  2:5176643/150‑1 (MQ=255)

                  tGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACCACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATaacaaa                                                                                                                                     >  1:8040290/1‑150 (MQ=255)

                  tGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATaacaaa                                                                                                                                     <  2:227794/150‑1 (MQ=255)

                    ttCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATaacaaaca                                                                                                                                   >  1:1314781/1‑150 (MQ=255)

                      cccATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAAcaaacaaa                                                                                                                                 >  2:5004410/1‑150 (MQ=255)

                      cccATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAAcaaacaaa                                                                                                                                 >  1:10161503/1‑150 (MQ=255)

                       ccATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAg                                                                                                                                >  1:1741141/1‑150 (MQ=255)

                       ccATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAg                                                                                                                                <  2:10380725/150‑1 (MQ=255)

                        cATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGc                                                                                                                               <  2:7725995/150‑1 (MQ=255)

                         aTTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCg                                                                                                                              >  2:1070106/1‑150 (MQ=255)

                          ttACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGa                                                                                                                             <  2:6769944/150‑1 (MQ=255)

                           tACGGTGGCAACACCGGTACGGACCAAACCATCGCAAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGAc                                                                                                                            <  2:10431163/150‑1 (MQ=255)

                            aCGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGAcc                                                                                                                           >  2:1935086/1‑150 (MQ=255)

                               gTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGGTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCa                                                                                                                        >  2:622579/1‑150 (MQ=255)

                                tGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCAt                                                                                                                       >  2:3363710/1‑150 (MQ=255)

                                 ggaaaCACGGGTAGGGACCAACCCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGTATCAGAAAAGCCGTAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTATTGCAAACGCGAAAATAACAAAAAACGCGACCGCATc                                                                                                                      <  2:5897673/147‑1 (MQ=255)

                                 ggaaaCACCGGTACGTACCAAACCATCGCAAATAATAAAAAAGGCGTAATTCTGTAAAACGTTTGTATCATAAAAGCCGTAAAATAATTTTGGTACT‑TCATGCTTTCCGATTAATGAAAACGCTAAAATAACAAACAAAGCGACCGCATc                                                                                                                      <  2:4279515/147‑1 (MQ=255)

                                  gCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATcc                                                                                                                     >  2:7050710/1‑150 (MQ=255)

                                    aaCACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCACACGCGACAATAACAAACAAAGCGACCGCATCCAt                                                                                                                   >  1:10669632/1‑150 (MQ=255)

                                    aaCACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCAt                                                                                                                   >  2:9273187/1‑150 (MQ=255)

                                         cGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCa                                                                                                              >  1:5387064/1‑150 (MQ=255)

                                               gtaaaaaaCCATCTAAAATAATAAACATGGCGTAAATCTGTAAAACGTTTGTATCATAAATGCCTTAAAATACTTCTGGGAGTGTCA‑GCTTTCGGTTTATTGCAAACGCGAAAATAACAAACAACGCGACCGCATCCATACGCACTCTGc                                                                                                        <  2:4127594/145‑1 (MQ=255)

                                                gACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGcc                                                                                                       >  1:398516/1‑150 (MQ=255)

                                                 aCCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGccc                                                                                                      >  1:6398966/1‑150 (MQ=255)

                                                  ccAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCg                                                                                                     >  1:1912301/1‑150 (MQ=255)

                                                   cAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGt                                                                                                    >  1:6555175/1‑150 (MQ=255)

                                                   cAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGt                                                                                                    >  1:5278156/1‑150 (MQ=255)

                                                   cAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGt                                                                                                    >  1:10512372/1‑150 (MQ=255)

                                                      aCCATCGCAAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGc                                                                                                 <  2:5562422/150‑1 (MQ=255)

                                                       ccATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCa                                                                                                >  2:3507351/1‑150 (MQ=255)

                                                       ccATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCa                                                                                                >  1:9245888/1‑150 (MQ=255)

                                                       ccATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCa                                                                                                >  1:9480079/1‑150 (MQ=255)

                                                        cATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCaa                                                                                               >  2:4476304/1‑150 (MQ=255)

                                                        cATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCaa                                                                                               >  2:6061620/1‑150 (MQ=255)

                                                         aTCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCaaa                                                                                              <  2:3733538/150‑1 (MQ=255)

                                                             ccAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAg                                                                                          >  1:4862787/1‑150 (MQ=255)

                                                                 taataaACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGa                                                                                      >  1:9797017/1‑150 (MQ=255)

                                                                 taataaACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGa                                                                                      >  1:3294876/1‑150 (MQ=255)

                                                                   ataaCCACGGCGTCAATCCTTAAAACTTTCGTATCATAAATGCCGTAAAATACTTGTGGTAGTTTCT‑GGTTTCCGATTAATCAAAACGCTAAAATAACAAAAAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATg                                                                                    <  1:2374457/150‑1 (MQ=255)

                                                                   ataaACATGGCGTCAATCAGTACAACGTTAGTATCATAAAAGCCGTAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGAAAACGCGACAATAACAAACAACGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATg                                                                                    <  2:1511620/150‑1 (MQ=255)

                                                                   ataaACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATg                                                                                    <  2:8184444/150‑1 (MQ=255)

                                                                   ataaACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATg                                                                                    >  1:9001135/1‑150 (MQ=255)

                                                                   ataaACAACCCGCAAATCAGCACAACGTTAGCAAAACAAAACCCGCAAAATACTTCTGGGACTGACA‑GCTTTCCCCTTAATCAAAACGCTAAAATAAAAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATg                                                                                    <  2:9639425/150‑1 (MQ=255)

                                                                     aaaCAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGc                                                                                  >  2:1652312/1‑150 (MQ=255)

                                                                      aaCAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCa                                                                                 >  2:1252330/1‑150 (MQ=255)

                                                                         aaGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACg                                                                              >  2:6886768/1‑150 (MQ=255)

                                                                                 aaTCTGTACAACGTTAGCATAACAAAAGCCGTAAAATACTCCTGGTACTGTCA‑GCGTCCCGATTAATGCAAACGCTAAAATAACAAACAACGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTaa                                                                      <  2:2399047/150‑1 (MQ=255)

                                                                                 aaTCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTaa                                                                      >  2:7193454/1‑150 (MQ=255)

                                                                                  aTCTGGACAAGGTTAGTATCAGAAAAGCCGTAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAAt                                                                     <  2:8922118/150‑1 (MQ=255)

                                                                                      ggACAACGTTTGGATCAGAAAAGCCGTAAAATACTTCTGGGACTTTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGaaa                                                                 <  2:2583524/150‑1 (MQ=255)

                                                                                       gACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAg                                                                >  1:9965115/1‑150 (MQ=255)

                                                                                       gACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAg                                                                >  2:9054972/1‑150 (MQ=255)

                                                                                         cAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGaa                                                              >  1:7522701/1‑150 (MQ=255)

                                                                                          aaCGTTAGGATCAGAAAAGCCGTAAAATACTTCTGGTACTGTCA‑CCGTTCCGATTAATGCAAACGCGAAAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAAc                                                             <  2:8026787/150‑1 (MQ=255)

                                                                                              ttAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCaa                                                         >  1:784609/1‑150 (MQ=255)

                                                                                                  gATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCa                                                     >  1:6203422/1‑150 (MQ=255)

                                                                                                        aaaaGCCGTAAAATAGTTGTGGGAATGTCT‑GGTTTCCTCTTATTGAAAACGCTAAAATAACAAAAAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCAc                                               <  1:1518455/150‑1 (MQ=255)

                                                                                                               ggAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTcac                                        >  1:4199511/1‑150 (MQ=255)

                                                                                                               ggAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTcac                                        >  1:2818261/1‑150 (MQ=255)

                                                                                                                 aaaaTACTTGTGGGAGTTTCA‑GGGTTCCGCTTAATGCAAACGCGAAAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACaa                                      <  1:2221291/150‑1 (MQ=255)

                                                                                                                   aaTACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACt                                    >  1:7365679/1‑150 (MQ=255)

                                                                                                                    aTACTTCTGGGACTGTCA‑GCGTTCCGATTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACtt                                   <  2:10173306/150‑1 (MQ=255)

                                                                                                                      aCTTCTGGGACTGTCA‑GCGTTCCGCTTAATGAAAACGCGACAATAACAAAAAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTAt                                 <  2:11100179/150‑1 (MQ=255)

                                                                                                                       cTTCTGGTACTGTCA‑GCGTTCCGATTATTGCAAACGCGACAATAACAAACAACGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATc                                <  1:5969567/150‑1 (MQ=255)

                                                                                                                          cTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCtt                             >  2:10913936/1‑150 (MQ=255)

                                                                                                                          cTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGAAAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCtt                             <  1:10011370/150‑1 (MQ=255)

                                                                                                                           tGGTACTGTCT‑GCGTTCCTATTAATGCAAACGCTACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTg                            <  2:7908165/150‑1 (MQ=255)

                                                                                                                           tGGTACTGTCA‑GCTTTCCTATTAATGAAAACCTGAAAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTg                            <  2:5330438/150‑1 (MQ=255)

                                                                                                                              gACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTa                         >  1:303214/1‑150 (MQ=255)

                                                                                                                                 tGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAAtt                      >  2:2930539/1‑150 (MQ=255)

                                                                                                                                   tCA‑GCGTTCCTCTTAATGCAAACGCGAAAATAAAAAACAAAGCGACCGCATCCATACGCACTCTGCACGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGa                    <  2:9830856/150‑1 (MQ=255)

                                                                                                                                   tCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGa                    >  2:4783018/1‑150 (MQ=255)

                                                                                                                                   tCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGa                    >  1:5645599/1‑150 (MQ=255)

                                                                                                                                    cA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAg                   <  2:6472820/150‑1 (MQ=255)

                                                                                                                                     ctgCGTTCGGCTTAATGCAAACGCTAAAATAACAAACAACGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAg                   <  2:3228842/148‑1 (MQ=255)

                                                                                                                                     caggtttCGGCTTAATGCAACCGCGACAATAACAAACACCGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAg                   <  1:10380972/145‑1 (MQ=255)

                                                                                                                                        cGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAgggg                >  2:4874848/1‑150 (MQ=255)

                                                                                                                                         gTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAggggg               >  1:6932727/1‑150 (MQ=255)

                                                                                                                                          ttCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGa              >  1:5841354/1‑150 (MQ=255)

                                                                                                                                             cgattAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGa           <  1:8298748/147‑1 (MQ=255)

                                                                                                                                              gCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGAc          <  2:8373230/150‑1 (MQ=255)

                                                                                                                                               catAATGCAAACGCGACAATAACAAACAACGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACt         <  1:10907508/148‑1 (MQ=255)

                                                                                                                                               cTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACt         >  1:5096599/1‑150 (MQ=255)

                                                                                                                                                  aaTGCAAACGCGACAATAACAAACAACGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGa      <  2:5677468/150‑1 (MQ=255)

                                                                                                                                                  aaTCCACACGCTAAAATAACAAACAACGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACACGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGa      <  1:3777096/150‑1 (MQ=255)

                                                                                                                                                   aTGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGAt     >  1:8150976/1‑150 (MQ=255)

                                                                                                                                                      cAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCt  >  2:7353757/1‑150 (MQ=255)

                                                                                                                                                      cAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCt  >  2:9472610/1‑150 (MQ=255)

CGCAACTTTGACCAGCCATGTTCCCATTACGGTGGCAACACCGGTACGGACCAAACCATCGCCAATAATAAACAAGGCGGCAATCAGGACAACGTTAGGATCAGAAAAGCCGGAAAATACTTCTGGGACTGTCA‑GCGTTCCGCTTAATGCAAACGCGACAATAACAAACAAAGCGACCGCATCCATACGCACTCTGCCCGTCGCAAACAAGACGATGGCAACGGCCAGTAATGAAAGAACCCAAATCAATTCACCGTTCACAACTTATCCTTGTTAATTGAGGGGGATGACTTGATTCT  >  NC_000913/2411092‑2411390

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