Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,096,388	T→C	5.5%	N2N (AAT→AAC)	yggU →	UPF0235 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,096,388	0	T	C	5.5%	234.7 / 3.1	72	N2N (AAT→AAC)	yggU	UPF0235 family protein
Reads supporting (aligned to +/- strand):  ref base T (35/33);  new base C (2/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GCCTGCTACCGGCAATGGGTGGGATTGATTTCTCGCCGATGATCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTCGCCATTACCGCGCCGCCGGTTGACGGCCAGGCCAACAG  >  NC_000913/3096247‑3096528                                                                                                                                              |                                                                                                                                             gccTGCTACGGGCAATGGGTGGGATTGATGTCTCGCCGATGATCCTCGTTCTGGTGCTGTATTTCATGAATATGGGTGTCGCAGAAGTATTATAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTGATGAATGCCGTaa                                                                                                                                       <  2:355511/149‑1 (MQ=255)     gCTACCGGCAATGGGTGGGATTGATTTCTCGCCGATGATCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGt                                                                                                                                   >  2:335012/1‑149 (MQ=255)     gCTACCGGCAATGGGTGGGATTGATTTCTCGCCGATGATCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGt                                                                                                                                   <  4:428516/149‑1 (MQ=255)        aCCGGCAATGGGTGGGATTGATTTCTCGCCGATGATCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTaa                                                                                                                                <  2:527483/149‑1 (MQ=255)              aaTGGGTGGGATTGATTTCTCGCCGATGATCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAAtgatga                                                                                                                          <  3:457580/149‑1 (MQ=255)               aTGGGTGGGATTGATTTCTCGCCGATGATCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGTAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGAc                                                                                                                         >  1:416378/1‑149 (MQ=255)               aTGGGTGGGATTGATTTCTCGCCGATGATCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGAc                                                                                                                         >  4:430152/1‑149 (MQ=255)                    tgggATTGATTTCTCGCCGATGATCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCt                                                                                                                    >  3:475066/1‑149 (MQ=255)                    tgggATTGATTTCTCGCCGATGATCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCt                                                                                                                    >  3:287441/1‑149 (MQ=255)                       gattgattTCTCGCCGATGATCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGt                                                                                                                 <  3:526348/149‑1 (MQ=255)                         ttgattTCTCGCCGATGATCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGttt                                                                                                               >  1:344486/1‑149 (MQ=255)                         ttgattTCTCGCCGATGATCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGttt                                                                                                               >  1:410114/1‑149 (MQ=255)                         ttgattTCTCGCCGATGATCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGttt                                                                                                               >  3:315245/1‑149 (MQ=255)                           gattTCTCGCCGATGATCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGACGTAGTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTa                                                                                                             <  2:416378/149‑1 (MQ=255)                           gattTCTCGCCGATGATCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTa                                                                                                             <  4:329897/149‑1 (MQ=255)                           gattTCTCGCCGATGATCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTa                                                                                                             <  2:25251/149‑1 (MQ=255)                            attTCTCGCCGATGATCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTAc                                                                                                            >  1:177151/1‑149 (MQ=255)                             tttCTCGCCGATGATCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACg                                                                                                           <  3:356083/149‑1 (MQ=255)                                ctcGCCGATGATCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGct                                                                                                        <  2:364020/149‑1 (MQ=255)                                   gCCGATGATCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCta                                                                                                     <  1:549464/149‑1 (MQ=255)                                     cGATGATCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCtata                                                                                                   <  1:265108/149‑1 (MQ=255)                                     cGATGATCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCtata                                                                                                   <  2:111355/149‑1 (MQ=255)                                          atCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAg                                                                                              >  3:154806/1‑149 (MQ=255)                                           tCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGc                                                                                             >  1:268009/1‑149 (MQ=255)                                           tCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGc                                                                                             <  3:430152/149‑1 (MQ=255)                                             cTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCg                                                                                           >  2:322044/1‑149 (MQ=255)                                             cTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCg                                                                                           >  2:359356/1‑149 (MQ=255)                                                 ttCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAag                                                                                       >  3:405342/1‑149 (MQ=255)                                                       tgctgTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAagccagcc                                                                                 >  3:221612/1‑149 (MQ=255)                                                           gTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGa                                                                             <  2:557992/149‑1 (MQ=255)                                                             aTGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGAtt                                                                           <  4:475066/149‑1 (MQ=255)                                                              tGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTc                                                                          >  2:526256/1‑149 (MQ=255)                                                               gTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCt                                                                         >  3:387732/1‑149 (MQ=255)                                                                  atcaATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTAtt                                                                      >  3:203077/1‑149 (MQ=255)                                                                  atcaATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTAtt                                                                      >  2:553399/1‑149 (MQ=255)                                                                  atcaATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTAtt                                                                      >  4:414225/1‑149 (MQ=255)                                                                  atcaATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTAtt                                                                      >  1:183107/1‑149 (MQ=255)                                                                   tcaATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTg                                                                     >  1:209668/1‑149 (MQ=255)                                                                     aaTATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTc                                                                   <  2:183107/149‑1 (MQ=255)                                                                      atatGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGAt                                                                            >  2:395348/1‑139 (MQ=255)                                                                      atatGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGAt                                                                            <  1:395348/139‑1 (MQ=255)                                                                              gtCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACAt                                                          <  1:553399/149‑1 (MQ=255)                                                                                cGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATgg                                                        <  4:73383/149‑1 (MQ=255)                                                                                    gAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGcgac                                                    >  1:337548/1‑149 (MQ=255)                                                                                       gTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGaa                                                 <  4:42188/149‑1 (MQ=255)                                                                                        tATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAg                                                >  3:151103/1‑149 (MQ=255)                                                                                          ttACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGtt                                              >  3:160448/1‑149 (MQ=255)                                                                                           tACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTa                                             <  2:177151/149‑1 (MQ=255)                                                                                           tACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTa                                             <  4:151103/149‑1 (MQ=255)                                                                                                gCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTc                                        >  2:514030/1‑149 (MQ=255)                                                                                                 cAACCGGAAATACGCTGCTGCCAGGGCTGTGCATGGCGTTATGAATGCCGTAACATTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTCg                                       <  2:451172/149‑1 (MQ=255)                                                                                                     cGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTCGCCAt                                   >  3:395471/1‑149 (MQ=255)                                                                                                       gAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTCGCCATTa                                 >  1:306372/1‑149 (MQ=255)                                                                                                                tgctgcCAGGGCTGTGGATGGCGTTATGAACGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTCGCCATTACCGcgccgc                        >  2:140552/1‑149 (MQ=255)                                                                                                                tgctgcCAGGGCTGTGGATGGCGTTATGAACGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTCGCCATTACCGcgccgc                        >  2:241626/1‑149 (MQ=255)                                                                                                                   tgcCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTCGCCATTACCGCGCCGCCgg                     <  2:306372/149‑1 (MQ=255)                                                                                                                   tgcCAGGGCTGCGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTCGCCATTACCGCGCCGCCgg                     <  1:569189/149‑1 (MQ=255)                                                                                                                     ccAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTCGCCATTACCGCGCCGCCGGtt                   <  3:415887/149‑1 (MQ=255)                                                                                                                      cAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTCGCCATTACCGCGCCGCCGGTTg                  >  4:292639/1‑149 (MQ=255)                                                                                                                        gggCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTCGCCATTACCGCGCCGCCGGTTGAc                <  4:387732/149‑1 (MQ=255)                                                                                                                        gggCTGCGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTCGCCATTACCGCGCCGCCGGTTGAc                <  2:50421/149‑1 (MQ=255)                                                                                                                         ggCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTCGCCATTACCGCGCCGCCGGTTGACg               >  3:263248/1‑149 (MQ=255)                                                                                                                             gtgGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTCGCCATTACCGCGCCGCCGGTTGACggcca           <  2:337548/149‑1 (MQ=255)                                                                                                                               ggATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTCGCCATTACCGCGCCGCCGGTTGACggccagg         <  1:206835/149‑1 (MQ=255)                                                                                                                               ggATGGCGTTATGAACGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTCGCCATTACCGCGCCGCCGGTTGACggccagg         <  1:241626/149‑1 (MQ=255)                                                                                                                               ggATGGCGTTATGAACGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTCGCCATTACCGCGCCGCCGGTTGACggccagg         <  1:140552/149‑1 (MQ=255)                                                                                                                                     cGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTCGCCATTACCGCGCCGCCGGTTGACGGCCAGGCCAACa   <  2:566987/149‑1 (MQ=255)                                                                                                                                      gTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTCGCCATTACCGCGCCGCCGGTTGACGGCCAGGCCAACAg  >  2:166910/1‑149 (MQ=255)                                                                                                                                      gTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTCGCCATTACCGCGCCGCCGGTTGACGGCCAGGCCAACAg  >  1:504475/1‑149 (MQ=255)                                                                                                                                      gTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTCGCCATTACCGCGCCGCCGGTTGACGGCCAGGCCAACAg  >  3:339348/1‑149 (MQ=255)                                                                                                                                       ttATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTCGCCATTACCGCGCCGCCGGTTGAc                >  2:100041/1‑134 (MQ=255)                                                                                                                                       ttATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTCGCCATTACCGCGCCGCCGGTTGAc                <  1:100041/134‑1 (MQ=255)                                                                                                                                              |                                                                                                                                             GCCTGCTACCGGCAATGGGTGGGATTGATTTCTCGCCGATGATCCTCGTTCTGCTGCTGTATGTCATCAATATGGGTGTCGCAGAAGTATTACAGGCAACCGGAAATATGCTGCTGCCAGGGCTGTGGATGGCGTTATGAATGCCGTAACAGTTAATGATGACGGTCTGGTTTTACGGCTCTATATTCAGCCGAAAGCCAGCCGTGATTCTATTGTCGGTTTACATGGCGACGAAGTTAAAGTCGCCATTACCGCGCCGCCGGTTGACGGCCAGGCCAACAG  >  NC_000913/3096247‑3096528

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