Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,367,513	C→A	10.5%	G134C (GGC→TGC)	pspF ←	psp operon transcriptional activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,367,513	0	C	A	10.5%	201.6 / 8.8	65	G134C (GGC→TGC)	pspF	psp operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (28/29);  new base A (5/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAAC  >  NC_000913/1367373‑1367661                                                                                                                                             |                                                                                                                                                     gcgcAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAAcgc                                                                                                                                              >  4:433194/1‑149 (MQ=255) gcgcAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCGGCAGCCAAcca                                                                                                                                              >  4:795669/1‑147 (MQ=255)  cgcAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCGGCAGCCAAccat                                                                                                                                             <  2:369030/149‑4 (MQ=255)   gcAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACgcgc                                                                                                                                            >  2:137782/1‑149 (MQ=255)     aGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCGGCAGCCAAccattgc                                                                                                                                          <  3:795669/149‑7 (MQ=255)          tggCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGtt                                                                                                                                     >  1:59615/1‑149 (MQ=255)            gCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCa                                                                                                                                   >  4:76074/1‑149 (MQ=255)            gCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCa                                                                                                                                   <  2:745947/149‑1 (MQ=255)               gTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCg                                                                                                                                <  4:819743/149‑1 (MQ=255)                 tGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTa                                                                                                                              <  1:889251/149‑1 (MQ=255)                 tGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTa                                                                                                                              <  2:195958/149‑1 (MQ=255)                   tACAACAGCAAAAGCCCGTCGGTTGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACTATCGCCGGGAGATCGGCATTCGTCGAGCATACCAACCGCACATTCACCTGCAAGGGTTGGCTGCTGCCAACGCGCTCCAGTTCACCGTACt                                                                                                                            >  4:848456/1‑149 (MQ=255)                        caTCAAAAGACAGTCAGTCGAGCAGGTCAACGCGAAAAGTGCCTTCATTGACCATCGCCGGGAAATCGGCCTACGTCGGCCATACCAACCGCACATTCACCTGCAATGTTTGACTACCGCCAACGCGCCCCAGTTCACCGTACTCAATc                                                                                                                       >  4:269903/1‑149 (MQ=255)                          tCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCAc                                                                                                                     >  4:171155/1‑149 (MQ=255)                            aaaaGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCAcgc                                                                                                                   <  3:521064/149‑1 (MQ=255)                            aaaaGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCAcgc                                                                                                                   <  3:691516/149‑1 (MQ=255)                             aaaGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCAcgcg                                                                                                                  >  2:432274/1‑149 (MQ=255)                                gCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCaa                                                                                                               <  2:123585/149‑1 (MQ=255)                                     tcggtcgAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAAtt                                                                                                          <  2:754542/149‑1 (MQ=255)                                        gtcgAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTtct                                                                                                       <  3:66477/149‑1 (MQ=255)                                         tcgAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCCCCGTACTCACTCCCGCGCAATAATTcctc                                                                                                      <  4:459031/149‑1 (MQ=255)                                         tcgAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCACCGGGAGATCGGCATCCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTtctc                                                                                                      <  4:346236/149‑1 (MQ=255)                                           gAGCAGGTCAGTGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCt                                                                                                    >  2:781348/1‑149 (MQ=255)                                                     gcgcGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACcatcat                                                                                          <  2:2148/149‑1 (MQ=255)                                                      cgcgAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCatcatc                                                                                         <  2:74841/149‑1 (MQ=255)                                                       gcgAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCg                                                                                        >  2:829238/1‑149 (MQ=255)                                                        cgAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCgg                                                                                       >  2:850384/1‑149 (MQ=255)                                                           aaaGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGc                                                                                    <  1:833579/149‑1 (MQ=255)                                                            aaGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGcc                                                                                   <  3:273238/149‑1 (MQ=255)                                                               tGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCAcgc                                                                                                                   >  4:399009/1‑114 (MQ=255)                                                               tGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCAcgc                                                                                                                   <  3:399009/114‑1 (MQ=255)                                                                 ccTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGc                                                                              >  1:654391/1‑149 (MQ=255)                                                                   ttCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCGATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGa                                                                            >  4:508915/1‑149 (MQ=255)                                                                      aTTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGtt                                                                         >  4:98298/1‑149 (MQ=255)                                                                            cATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCaa                                                                   >  4:183487/1‑149 (MQ=255)                                                                             aTCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAg                                                                  >  4:813075/1‑149 (MQ=255)                                                                              tcgtcgGGAGATCGGCATTCGCCGCGCATACCAACCGCACATTCCCCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCAGCGGTGCCGTAGCGAGTTCATCAAGa                                                                 <  4:241699/149‑1 (MQ=255)                                                                               cGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGaa                                                                <  4:707583/149‑1 (MQ=255)                                                                                 ccGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAAt                                                              <  1:468543/149‑1 (MQ=255)                                                                                    ggAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGc                                                           <  4:327239/149‑1 (MQ=255)                                                                                        aTCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTAc                                                       >  2:139891/1‑149 (MQ=255)                                                                                         tCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTAcc                                                      >  1:154250/1‑149 (MQ=255)                                                                                                           aCCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGc                                                                                    >  3:715509/1‑101 (MQ=255)                                                                                                           aCCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGc                                                                                    <  4:715509/101‑1 (MQ=255)                                                                                                           aCCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCaaa                                    <  3:813075/149‑1 (MQ=255)                                                                                                            ccAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAAt                                   <  3:653029/149‑1 (MQ=255)                                                                                                            ccAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAAt                                   >  2:149799/1‑149 (MQ=255)                                                                                                             cAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAAtc                                  >  1:685357/1‑149 (MQ=255)                                                                                                             cAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAAtc                                  >  3:104969/1‑149 (MQ=255)                                                                                                                  gCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCt                             <  4:104969/149‑1 (MQ=255)                                                                                                                   cacaTTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCt                             >  2:269748/1‑148 (MQ=255)                                                                                                                   cacaTTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCt                             <  1:269748/148‑1 (MQ=255)                                                                                                                   cacaTTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTg                            <  3:183487/149‑1 (MQ=255)                                                                                                                      aTTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGAt                         >  1:765336/1‑149 (MQ=255)                                                                                                                          aCCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACg                     >  4:640696/1‑149 (MQ=255)                                                                                                                            cTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGtt                   <  3:98298/149‑1 (MQ=255)                                                                                                                               agcgcgTTGGCGGCAGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGttttt                >  3:391438/6‑149 (MQ=255)                                                                                                                                gcgcgTTGGCGGCAGCCCACGCGCTACAGTTCCCCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGATCAAATCTCCCTGGATGACGTTTTTg               >  2:34803/5‑149 (MQ=255)                                                                                                                                gcgcgTTGGCGGCAGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTg               >  4:77911/5‑149 (MQ=255)                                                                                                                                aaTGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTg               >  4:206862/1‑149 (MQ=255)                                                                                                                                    gTTGGCGGCAGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCa           >  1:281978/1‑149 (MQ=255)                                                                                                                                     ttGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCAc          >  4:109323/1‑149 (MQ=255)                                                                                                                                        gCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCTCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTcc                               >  3:576723/1‑125 (MQ=255)                                                                                                                                        gCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTcc                               <  4:576723/125‑1 (MQ=255)                                                                                                                                         cTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGt      <  1:137782/149‑1 (MQ=255)                                                                                                                                             cgccAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAAc  <  1:635397/149‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                     GCGCAGTGGTGGCAGTTGTACAACATCAAAAGCCAGTCGGTCGAGCAGGTCAGCGCGAAAAGTGCCTTCATTGACCATCGCCGGGAGATCGGCATTCGTCGCGCATACCAACCGCACATTCACCTGCAATGGTTGGCTGCCGCCAACGCGCTCCAGTTCACCGTACTCAATCACGCGCAATAATTTCTCCTGCACCATCATCGGTGCCGTAGCGAGTTCATCAAGAAATAGCGTACCGCCGTCGGCACGTTCAAATCTCCCTGGATGACGTTTTTGCGCACCGGTAAAC  >  NC_000913/1367373‑1367661

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