Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,499,547	G→T	5.4%	D34Y (GAC→TAC)	yhfL →	small lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,499,547	0	G	T	5.4%	255.9 / 4.7	74	D34Y (GAC→TAC)	yhfL	small lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (50/20);  new base T (2/2);  total (52/22)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

TAATCAATGAATATTAATTAATATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCGAGGGCGTATTTGTTTCTATGT  >  NC_000913/3499404‑3499695                                                                                                                                                |                                                                                                                                                     taATCAATGAATATTAATTAATATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTAc                                                                                                                                                 <  2:532208/149‑1 (MQ=255) taATCAATGAATATTAATTAATATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTAc                                                                                                                                                 <  4:318234/149‑1 (MQ=255)                 attaatATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCaa                                                                                                                                 >  1:119182/1‑149 (MQ=255)                  ttaatATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAat                                                                                                                                >  2:568664/1‑149 (MQ=255)                          taataCATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAatttctattt                                                                                                                        >  4:333092/1‑149 (MQ=255)                           aataCATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAatttctatttc                                                                                                                       >  3:306074/1‑149 (MQ=255)                                    ggAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCAt                                                                                                              >  4:845576/1‑149 (MQ=255)                                    ggAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCAt                                                                                                              <  3:333092/149‑1 (MQ=255)                                       aTGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTgg                                                                                                           <  2:567344/149‑1 (MQ=255)                                          tAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCgg                                                                                                        >  1:575160/1‑149 (MQ=255)                                                aaCAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCgg                                                                                                  >  2:875807/1‑149 (MQ=255)                                                 aCAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGt                                                                                                 >  3:849825/1‑149 (MQ=255)                                                   attgaatttaAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTcc                                                                                               >  4:23362/8‑149 (MQ=255)                                                   aaaTTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTcc                                                                                               >  4:60633/1‑149 (MQ=255)                                                   aaaTTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTcc                                                                                               >  4:437107/1‑149 (MQ=255)                                                          ttaAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGca                                                                                        >  4:762171/1‑149 (MQ=255)                                                          ttaAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGca                                                                                        >  1:121201/1‑149 (MQ=255)                                                          ttaAAGTTGCACCGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGca                                                                                        >  3:66334/1‑149 (MQ=255)                                                            aaaGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGcaca                                                                                      <  2:331439/149‑1 (MQ=255)                                                              aGTTGCACTGGTAGGTGCAGTACTCGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGt                                                                                    >  4:678032/1‑149 (MQ=255)                                                                ttGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTaa                                                                                  >  4:379457/1‑149 (MQ=255)                                                                ttGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTaa                                                                                  <  1:653458/149‑1 (MQ=255)                                                                        gTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGa                                                                          >  2:543982/1‑149 (MQ=255)                                                                             tGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGaa                                                                     >  2:869417/1‑149 (MQ=255)                                                                              gCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGaaa                                                                    >  1:866843/1‑149 (MQ=255)                                                                                    cTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCt                                                              >  2:207817/1‑149 (MQ=255)                                                                                     tGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACTACTACCTGCTGCACCCGGCAGTTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTa                                                             >  4:522048/1‑149 (MQ=255)                                                                                     tGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACTACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTa                                                             >  4:831650/1‑149 (MQ=255)                                                                                         tACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAacac                                                         >  3:496839/1‑149 (MQ=255)                                                                                         tACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAacac                                                         >  4:520702/1‑149 (MQ=255)                                                                                         tACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAacac                                                         >  3:414997/1‑149 (MQ=255)                                                                                           cGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCg                                                       >  2:136903/1‑149 (MQ=255)                                                                                            gTTAACTGCATGTACTGGTCATATTGGAAACCGTGGTAAAACCTGATCTTACGACTACCTGCTGCACTCTGCGAGTTCTATTTCTAAAATCATTGGCGGGTGGGGGCCTACTACAGAGTAAGTCTTTGTCAGAAATCGGCTAAGACCgg                                                      >  4:194440/1‑149 (MQ=255)                                                                                              tAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGtt                                                    <  4:552537/149‑1 (MQ=255)                                                                                                   gCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAAtt                                               >  1:371298/1‑149 (MQ=255)                                                                                                     aTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGc                                             <  1:207817/149‑1 (MQ=255)                                                                                                     aTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGc                                             <  4:386658/149‑1 (MQ=255)                                                                                                     aTGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGc                                             >  1:693025/1‑149 (MQ=255)                                                                                                       gTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGccc                                           <  4:428227/149‑1 (MQ=255)                                                                                                         aCTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTAc                                                                                            <  3:224777/98‑1 (MQ=255)                                                                                                         aCTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTAc                                                                                            >  4:224777/1‑98 (MQ=255)                                                                                                         aCTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGa                                         >  4:415302/1‑149 (MQ=255)                                                                                                          cTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGAc                                        >  1:51966/1‑149 (MQ=255)                                                                                                          cTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGAc                                        <  1:568664/149‑1 (MQ=255)                                                                                                              tCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGccc                                    >  2:713859/1‑149 (MQ=255)                                                                                                                atatTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAAt                                                >  4:830543/1‑135 (MQ=255)                                                                                                                atatTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAAt                                                <  3:830543/135‑1 (MQ=255)                                                                                                                   ttGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTgg                               >  4:423052/1‑149 (MQ=255)                                                                                                                    tGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGa                              >  3:418646/1‑149 (MQ=255)                                                                                                                      aaaaCCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGAtt                            >  4:189432/1‑149 (MQ=255)                                                                                                                      aaaaCCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGAtt                            >  1:810616/1‑149 (MQ=255)                                                                                                                         aCCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGATTTTc                         >  1:722694/1‑149 (MQ=255)                                                                                                                          ccGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGATTTTcc                        >  4:296838/1‑149 (MQ=255)                                                                                                                           cGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCg                       >  2:519651/1‑149 (MQ=255)                                                                                                                           cGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCg                       >  2:636108/1‑149 (MQ=255)                                                                                                                              gATAAGAACTGCTCTTACTACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCGAgg                    <  3:522048/149‑1 (MQ=255)                                                                                                                              gATAAGAACTGCTCTTACTACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCGAgg                    <  3:831650/149‑1 (MQ=255)                                                                                                                               aTAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCGAggg                   >  3:427652/1‑149 (MQ=255)                                                                                                                               aTAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCGAggg                   >  3:359837/1‑149 (MQ=255)                                                                                                                                tAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTctct                                                                            >  1:596388/1‑91 (MQ=255)                                                                                                                                tAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTctct                                                                            <  2:596388/91‑1 (MQ=255)                                                                                                                                tAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCGAGGGc                  <  1:543982/149‑1 (MQ=255)                                                                                                                                 aaGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCGAGGGCg                 >  1:62501/1‑149 (MQ=255)                                                                                                                                   gAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGtt                                                    <  2:215861/112‑1 (MQ=255)                                                                                                                                   gAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGtt                                                    >  1:215861/1‑112 (MQ=255)                                                                                                                                   gAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGAtttt                          <  2:219644/138‑1 (MQ=255)                                                                                                                                   gAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGAtttt                          >  1:219644/1‑138 (MQ=255)                                                                                                                                    aaCTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCGAGGGCGTAt              <  3:189432/149‑1 (MQ=255)                                                                                                                                          tctTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCGAGGGCGTATTTGttt        >  3:624626/1‑149 (MQ=255)                                                                                                                                          tctTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCGAGGGCGTATTTGttt        >  1:61253/1‑149 (MQ=255)                                                                                                                                          tctTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCGAGGGCGTATTTGttt        >  1:566483/1‑149 (MQ=255)                                                                                                                                            ttACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCGAGGGCGTATTTGTTTCt      >  1:881315/1‑149 (MQ=255)                                                                                                                                              aCGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCGAGGGCGTATTTGTTTCTAt    <  3:845576/149‑1 (MQ=255)                                                                                                                                              aCGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCGAGGGCGTATTTGTTTCTAt    <  4:306074/149‑1 (MQ=255)                                                                                                                                                gACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCGAGGGCGTATTTGTTTCTATGt  >  3:293058/1‑149 (MQ=255)                                                                                                                                                gACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCGAGGGCGTATTTGTTTCTATGt  >  2:462115/1‑149 (MQ=255)                                                                                                                                                gACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCGAGGGCGTATTTGTTTCTATGt  >  1:387285/1‑149 (MQ=255)                                                                                                                                                gACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCGAGGGCGTATTTGTTTCTATGt  >  3:511702/1‑149 (MQ=255)                                                                                                                                                gACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTCAATTGCCCGACGCCCTCTTGATTTTCCGAGGGCGTATTTGTTTCTATGt  >  3:33014/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                     TAATCAATGAATATTAATTAATATATAATACATAGGGAATGTAAATGAACAAATTTATTAAAGTTGCACTGGTAGGTGCAGTACTGGCTACGTTAACTGCATGTACTGGTCATATTGAAAACCGTGATAAGAACTGCTCTTACGACTACCTGCTGCACCCGGCAATTTCTATTTCTAAAATCATTGGCGGTTGCGGTCCTACTGCACAGTAAGTCTCTGACAGAAATCGGCTAACACCGGTTTAATTGCCCGACGCCCTCTGGATTTTCCGAGGGCGTATTTGTTTCTATGT  >  NC_000913/3499404‑3499695

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