Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,304,469	C→G	5.0%	R212R (CGC→CGG)	oppC →	oligopeptide ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,304,469	0	C	G	5.0%	351.5 / 8.4	99	R212R (CGC→CGG)	oppC	oligopeptide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (48/46);  new base G (3/2);  total (51/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

CTGCTGATTTTCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTGATGGC  >  NC_000913/1304326‑1304616                                                                                                                                                |                                                                                                                                                    ctgctgATTTTCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACAt                                                                                                                                                <  4:35521/149‑1 (MQ=255)  tgctgATTTTCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACAtt                                                                                                                                               >  1:468401/1‑149 (MQ=255)  tgctgATTTTCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACAtt                                                                                                                                               >  1:802338/1‑149 (MQ=255)   gctgATTTTCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATtg                                                                                                                                              >  2:499696/1‑149 (MQ=255)   gctgATTTTCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATtg                                                                                                                                              >  1:472073/1‑149 (MQ=255)    ctgATTTTCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATtgt                                                                                                                                             <  4:550593/149‑1 (MQ=255)    ctgATTTTCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATtgt                                                                                                                                             <  2:68405/149‑1 (MQ=255)          ttCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGaa                                                                                                                                       <  1:47758/149‑1 (MQ=255)          ttCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGaa                                                                                                                                       <  2:239002/149‑1 (MQ=255)                 cGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTc                                                                                                                                <  4:609248/149‑1 (MQ=255)                  gATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCg                                                                                                                               >  2:359749/1‑149 (MQ=255)                  gATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCg                                                                                                                               >  3:578054/1‑149 (MQ=255)                   aTAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCgg                                                                                                                              >  2:756846/1‑149 (MQ=255)                       gCATGGTTTCCTGGCTGGATCTGGCTCGGATTGTGCGTGGGCACACCCTGAGTCTGCAGCGCAAAGGGTTTATTGAGGCGGCACAAGTTGGCGGAGTATCGACCTCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGtgtg                                                                                                                          >  2:285917/1‑149 (MQ=255)                        cATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTgtgg                                                                                                                         <  2:1394/149‑1 (MQ=255)                         aTGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTgtggt                                                                                                                        >  3:746854/1‑149 (MQ=255)                          tGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTgtggtg                                                                                                                       >  2:535985/1‑149 (MQ=255)                           ggTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTgtggtgg                                                                                                                      >  1:154652/1‑149 (MQ=255)                            gTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTgtggtggt                                                                                                                     <  4:269183/149‑1 (MQ=255)                            gTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTgtggtggt                                                                                                                     <  4:78296/149‑1 (MQ=255)                             tttCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTgtggtggtg                                                                                                                    <  1:361593/149‑1 (MQ=255)                                   ggctggATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTAc                                                                                                              <  4:384247/149‑1 (MQ=255)                                   ggctggATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTAc                                                                                                              <  1:535985/149‑1 (MQ=255)                                     ctggATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGc                                                                                                            <  2:309164/149‑1 (MQ=255)                                       ggATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCAt                                                                                                          >  2:505085/1‑149 (MQ=255)                                       ggATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGCATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCAt                                                                                                          >  4:44551/1‑149 (MQ=255)                                                tCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTgg                                                                                                 >  4:221360/1‑149 (MQ=255)                                                    aTTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGcc                                                                                             >  2:82467/1‑149 (MQ=255)                                                     ttGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGccc                                                                                            <  2:672628/149‑1 (MQ=255)                                                          cGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGCCTACGCATCACTTCTGGTGCCCAGCAt                                                                                       >  2:710952/1‑149 (MQ=255)                                                              ggCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATc                                                                                   >  3:40591/1‑149 (MQ=255)                                                              ggCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATc                                                                                   <  2:472073/149‑1 (MQ=255)                                                               gCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGc                                                                                                            <  4:79487/123‑1 (MQ=255)                                                               gCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGc                                                                                                            >  3:79487/1‑123 (MQ=255)                                                                cAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCct                                                                                 >  4:145072/1‑149 (MQ=255)                                                                cAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCct                                                                                 >  2:450966/1‑149 (MQ=255)                                                                 aaaCCCTGAGTCTGAAGCGCAAAGAGTTTATTGACGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGCGCCGAATGTACTCGGTGTGGTGGTAGTCTACGCATCACTTCTGGTGCCCAGCATCATCctc                                                                                >  3:532726/1‑149 (MQ=255)                                                                  aaCCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCg                                                                                                                               <  2:224393/101‑1 (MQ=255)                                                                  aaCCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCg                                                                                                                               >  1:224393/1‑101 (MQ=255)                                                                   aCCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCtt                                                                              <  2:154652/149‑1 (MQ=255)                                                                   aCCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCtt                                                                              <  2:422267/149‑1 (MQ=255)                                                                   aCCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCtt                                                                              <  1:499696/149‑1 (MQ=255)                                                                   aCCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCtt                                                                              <  2:243864/149‑1 (MQ=255)                                                                    cccTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCttt                                                                             <  3:449750/149‑1 (MQ=255)                                                                    cccTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCttt                                                                             >  2:272739/1‑149 (MQ=255)                                                                     ccTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTg                                                                            >  3:18751/1‑149 (MQ=255)                                                                     ccTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTg                                                                            >  2:812428/1‑149 (MQ=255)                                                                             tGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTc                                                                    <  2:306036/149‑1 (MQ=255)                                                                             tGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTc                                                                    <  1:756846/149‑1 (MQ=255)                                                                              gAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTcc                                                                   >  2:13756/1‑149 (MQ=255)                                                                               aaGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCt                                                                  >  1:282222/1‑149 (MQ=255)                                                                               aaGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGGCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCt                                                                  >  3:660451/1‑149 (MQ=255)                                                                                aGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCtt                                                                 <  3:658524/149‑1 (MQ=255)                                                                                    cAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCt                                                             >  4:95/1‑149 (MQ=255)                                                                                      aaGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCAt                                                                                       <  2:756998/121‑1 (MQ=255)                                                                                      aaGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCAt                                                                                       >  1:756998/1‑121 (MQ=255)                                                                                         agTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTg                                                        <  1:82467/149‑1 (MQ=255)                                                                                           tttATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTggg                                                      >  2:506448/1‑149 (MQ=255)                                                                                           tttATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTggg                                                      >  3:188786/1‑149 (MQ=255)                                                                                              aTTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGtt                                                   >  4:83803/1‑149 (MQ=255)                                                                                               ttGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTg                                                  <  1:13756/149‑1 (MQ=255)                                                                                                tGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTgg                                                 <  3:673267/149‑1 (MQ=255)                                                                                                tGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTgg                                                 <  4:18751/149‑1 (MQ=255)                                                                                                tGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTgg                                                 <  3:145072/149‑1 (MQ=255)                                                                                                tGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTgg                                                 <  1:450966/149‑1 (MQ=255)                                                                                                tGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATACGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTgg                                                 <  3:514090/149‑1 (MQ=255)                                                                                                   ggcggcACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTgg                                                 <  1:309098/146‑1 (MQ=255)                                                                                                   ggcggcACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTgg                                                 >  2:309098/1‑146 (MQ=255)                                                                                                    gcggcACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTAc                                             >  2:479154/1‑149 (MQ=255)                                                                                                    gcggcACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTAc                                             >  3:785822/1‑149 (MQ=255)                                                                                                       gcACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCa                                          >  1:657556/1‑149 (MQ=255)                                                                                                       gcACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCa                                          >  2:121092/1‑149 (MQ=255)                                                                                                       gcACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCa                                          >  2:302138/1‑149 (MQ=255)                                                                                                          caAGTTGGGGGTGGATCGAACTTGGGCATTGTTTTTCGCCACATTGTGCCGGATTTACTCGGTGTGGTTGTTGTCTACGCATCACTTCTGGTTGCCAGCATGAACCTCTTTTAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAaga                                       <  4:92721/149‑1 (MQ=255)                                                                                                          caAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAaga                                       <  2:282222/149‑1 (MQ=255)                                                                                                          caAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAaga                                       <  4:315580/149‑1 (MQ=255)                                                                                                          caAGTTGGCGGTGTATCGACATAGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAaga                                       <  4:660451/149‑1 (MQ=255)                                                                                                                ggCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGtt                                 <  4:78209/149‑1 (MQ=255)                                                                                                                   ggTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAag                              <  4:746854/149‑1 (MQ=255)                                                                                                                   ggTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGCCTACGCATCACTTCTGGTGCCCAGCATGCTCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAag                              <  1:710952/149‑1 (MQ=255)                                                                                                                        aTCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCt                         <  3:95/149‑1 (MQ=255)                                                                                                                         tggACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTg                        >  2:159804/3‑149 (MQ=255)                                                                                                                           gACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTggg                      <  3:376795/149‑1 (MQ=255)                                                                                                                            aCATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTgg                                                 <  3:662953/121‑1 (MQ=255)                                                                                                                            aCATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTgg                                                 >  4:662953/1‑121 (MQ=255)                                                                                                                                 gggCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCAt                >  3:108854/1‑149 (MQ=255)                                                                                                                                     aTTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCt            >  2:677879/1‑149 (MQ=255)                                                                                                                                      ttGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTg           >  2:61702/1‑149 (MQ=255)                                                                                                                                        gTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAg         <  3:221360/149‑1 (MQ=255)                                                                                                                                        gTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAg         >  1:343607/1‑149 (MQ=255)                                                                                                                                         ttATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGt        <  1:285917/149‑1 (MQ=255)                                                                                                                                         ttATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGt        >  3:126292/1‑149 (MQ=255)                                                                                                                                         acaTTCGGCACAATGTGGCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGt        >  2:704582/3‑149 (MQ=255)                                                                                                                                          tATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATACTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTg       >  4:656267/1‑149 (MQ=255)                                                                                                                                           aTTCGGCACAATGTGGCGTATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTa                                              <  2:624946/109‑1 (MQ=255)                                                                                                                                           aTTCGGCACAATGTGGCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGTATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTa                                              >  1:624946/1‑109 (MQ=255)                                                                                                                                           aTTCGGCACAATGTGGCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGt                    >  4:360409/1‑135 (MQ=255)                                                                                                                                           aTTCGGCACAATGTGGCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGt                    <  3:360409/135‑1 (MQ=255)                                                                                                                                               gCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTGATGgc  >  1:613522/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                    CTGCTGATTTTCGTGGCGATAGGCATGGTTTCCTGGCTGGATATGGCTCGTATTGTGCGTGGGCAAACCCTGAGTCTGAAGCGCAAAGAGTTTATTGAGGCGGCACAAGTTGGCGGTGTATCGACATCGGGCATTGTTATTCGCCACATTGTGCCGAATGTACTCGGTGTGGTGGTGGTCTACGCATCACTTCTGGTGCCCAGCATGATCCTCTTTGAATCTTTCCTTAGCTTCCTGGGGTTGGGTACGCAAGAGCCGTTAAGCAGCTGGGGTGCATTGCTGAGTGATGGC  >  NC_000913/1304326‑1304616

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