BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,113,475	G→A	6.0%	L61L (TTG→TTA) ‡	uspD →	stress‑induced protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,113,475	0	G	A	6.0%	368.1 / 7.8	116	L61L (TTG→TTA) ‡	uspD	stress‑induced protein
Reads supporting (aligned to +/- strand): ref base G (52/57); new base A (3/4); total (55/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

Only 100 of 123 total aligned reads displayed.

GGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCTGTGGTCA  >  NC_000913/4113327‑4113624

ggAATGAAGAAGTTGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgt                                                                                                                                                      <  1:321642/150‑1 (MQ=255)

  aaTGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACTCTCACCTGACGTTAATTCATATTGATTATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                                                                    <  1:156183/150‑1 (MQ=255)

  aaTGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCGGTGTTTCTACTTCCCTGCAACAGTAGATATTCTTCAAtttttg                                                                                                                                                    <  2:767487/150‑1 (MQ=255)

  aaTCAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGTCGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                                                                    <  2:75932/150‑1 (MQ=255)

   aTGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTga                                                                                                                                                   >  1:1332986/1‑150 (MQ=255)

     gtagaagaTGCGTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTAACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACTGAAGATATTCTTCAATTGTTgaag                                                                                                                                                 <  2:854166/148‑1 (MQ=255)

     gaagaaGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaat                                                                                                                                                 >  2:1335171/1‑149 (MQ=255)

     gaagaaGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaag                                                                                                                                                 <  2:162248/150‑1 (MQ=255)

     gaagaaGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaag                                                                                                                                                 <  1:1201633/150‑1 (MQ=255)

      aagaagATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaaga                                                                                                                                                <  1:240088/150‑1 (MQ=255)

       agaagaTGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaagaa                                                                                                                                               <  2:1557919/150‑1 (MQ=255)

           gaTGGCTTCCTGGTCAATAAAGCCCTGGCGCTCGCCAGACATACTGCCGCTCACCAGACGGTATTTCATCTTGCTGATGGCTTAAGCGAGGTGTATCCGGTTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAg                                                                                                                                           <  2:261575/150‑1 (MQ=255)

           gaTGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAg                                                                                                                                           >  1:1726218/1‑150 (MQ=255)

            aTGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatct                                                                                                                                          <  2:1636925/150‑4 (MQ=255)

             tGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGAAGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTc                                                                                                                                         <  2:1387539/150‑1 (MQ=255)

               ccTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCgg                                                                                                                                       <  1:444095/150‑1 (MQ=255)

                 ttACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGAt                                                                                                                                     <  1:770627/150‑1 (MQ=255)

                  tACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGATGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATa                                                                                                                                    <  2:654740/150‑1 (MQ=255)

                     tGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACa                                                                                                                                 <  1:404027/150‑1 (MQ=255)

                     tGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACa                                                                                                                                 <  2:1702995/150‑1 (MQ=255)

                      ggTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTTCAACAGAAGATATTCTTCAATTGTTGAAGAAAAAGTCGGATAACaa                                                                                                                                <  1:13589/150‑1 (MQ=255)

                      ggTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACaa                                                                                                                                >  1:1538499/1‑150 (MQ=255)

                       gTGAATAAAGCCCTGGAGCTCTCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAg                                                                                                                               <  2:746035/150‑1 (MQ=255)

                        tGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACatga                                                                                                                              >  1:906668/1‑147 (MQ=255)

                        tGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGc                                                                                                                              >  1:900646/1‑150 (MQ=255)

                        tGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGc                                                                                                                              <  2:524767/150‑1 (MQ=255)

                         gAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCt                                                                                                                             >  2:734814/1‑150 (MQ=255)

                         gAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATa                                                                                                                                             <  1:571305/134‑1 (MQ=255)

                         gAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATa                                                                                                                                             >  2:571305/1‑134 (MQ=255)

                          aaTAAAGGCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATATCAAGCTg                                                                                                                            <  1:1107322/150‑1 (MQ=255)

                          aaTAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATAATGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTg                                                                                                                            <  1:63980/150‑1 (MQ=255)

                             aaaGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCTGATAACAAGCTGtat                                                                                                                         >  1:282726/1‑150 (MQ=255)

                                gCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATaaa                                                                                                                      >  2:1064390/1‑150 (MQ=255)

                                  ccgcgAGCTCGCCAGACATAAAGCCGCTCACCTGACGTTAATTCATATTGATGCTGGCTTAAGCGACATGGATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACt                                                                                                                    <  2:1049822/146‑1 (MQ=255)

                                  ccTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACt                                                                                                                    <  2:685731/150‑1 (MQ=255)

                                  ccTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACt                                                                                                                    <  1:756511/150‑1 (MQ=255)

                                   cTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACtt                                                                                                                   >  1:102000/1‑149 (MQ=255)

                                      gAGCTCGCCAGACATAATGACTCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACg                                                                                                                <  1:1170654/150‑1 (MQ=255)

                                           cGCCAGACATAATGCCGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaaaa                                                                                                           >  2:1303812/1‑150 (MQ=255)

                                             ccAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCTAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAata                                                                                                         >  2:1122571/1‑150 (MQ=255)

                                                 aCATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCAGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCa                                                                                                     <  2:1096611/150‑1 (MQ=255)

                                                  cATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCaa                                                                                                    >  1:362195/1‑150 (MQ=255)

                                                     aaTGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATgg                                                                                                 <  2:163094/150‑1 (MQ=255)

                                                     aaTGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATgg                                                                                                 >  1:461428/1‑150 (MQ=255)

                                                        gACGCTCACCTGAAGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATACCAAACTGTATAAACTGACGAAAAATATACAATGGCCg                                                                                              >  1:439801/1‑150 (MQ=255)

                                                          cGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGaa                                                                                            >  2:1337644/1‑150 (MQ=255)

                                                              cACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACa                                                                                        >  2:1042114/1‑150 (MQ=255)

                                                               aCCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACaa                                                                                       <  2:1538499/150‑1 (MQ=255)

                                                                  tGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAAc                                                                                    <  2:1682099/150‑1 (MQ=255)

                                                                     cGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGc                                                                                 >  1:1202397/1‑150 (MQ=255)

                                                                             cATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGaa                                                                         >  1:490810/1‑150 (MQ=255)

                                                                             cATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGaa                                                                         >  1:1613461/1‑150 (MQ=255)

                                                                             cATATTGATCATGGCTAAACCGAGTTGAATCCGGGTATCAACTTCCCTGCAACAGAAGATATTCTTAAATTGTTGAAGAATAATTAGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGTCCGAAGACAAAACTGCGTATAGaa                                                                         >  2:13933/1‑150 (MQ=255)

                                                                                atTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgc                                                                      >  2:421652/1‑150 (MQ=255)

                                                                                  tgatgaTGGCTTAAGCGAGTTGTATCCGGGTATCTACTTACCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                                                                    >  1:1032061/1‑150 (MQ=255)

                                                                                    atgatgGCTTAATGGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAATATATTCTTCAATTGTTGAAGAATAAGTCGGATATCAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGgag                                                                  <  2:287336/150‑1 (MQ=255)

                                                                                     tgatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGgaga                                                                 <  2:144622/150‑1 (MQ=255)

                                                                                       atgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGaaa                                                               >  1:1427671/1‑150 (MQ=255)

                                                                                        tgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAAt                                                              <  1:87814/150‑1 (MQ=255)

                                                                                        tgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAAt                                                              <  1:987151/150‑1 (MQ=255)

                                                                                         ggCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATg                                                             <  1:1064390/150‑1 (MQ=255)

                                                                                          gCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGc                                                            >  1:870712/1‑150 (MQ=255)

                                                                                           cTTAATCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGcc                                                           <  1:1225075/150‑1 (MQ=255)

                                                                                           cTTAAGCGAGTTGTATCCGGGTCTCTACTTCCCTGCCACAGAAGATCTTCTTCAATTGTTGAAGAATAAGTCGTATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGcc                                                           >  1:1047022/1‑150 (MQ=255)

                                                                                           cTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGcc                                                           >  1:695196/1‑150 (MQ=255)

                                                                                           cTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGcc                                                           >  1:382126/1‑150 (MQ=255)

                                                                                              aaGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGa                                                        <  2:853537/150‑1 (MQ=255)

                                                                                              aaGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGa                                                        >  1:536784/1‑150 (MQ=255)

                                                                                                 cGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAac                                                     >  2:245461/1‑150 (MQ=255)

                                                                                                 cGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTCTATAAACTGACGAAAAATATTCAATGGCAGAAGACAAAACTGCGAATTGAACGCGGAGAAATGCCGGAAac                                                     >  1:120874/1‑150 (MQ=255)

                                                                                                    gTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACAct                                                  >  2:1316522/1‑150 (MQ=255)

                                                                                                      tGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACActgc                                                <  1:245461/150‑1 (MQ=255)

                                                                                                       gTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACActgct                                               >  1:1282643/1‑150 (MQ=255)

                                                                                                           ccGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATg                                                             >  2:1525942/1‑132 (MQ=255)

                                                                                                           ccGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATg                                                             <  1:1525942/132‑1 (MQ=255)

                                                                                                              ggTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATTCCGGAAACACTGCTGGAAAtt                                        <  2:102000/150‑1 (MQ=255)

                                                                                                                  tCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                                                                    <  2:1483852/118‑1 (MQ=255)

                                                                                                                  tCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                                                                    >  1:1483852/1‑118 (MQ=255)

                                                                                                                  tCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                                                                    <  2:1477173/118‑1 (MQ=255)

                                                                                                                  tCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                                                                    >  1:1477173/1‑118 (MQ=255)

                                                                                                                   ctactaCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCa                                   >  2:1036739/1‑150 (MQ=255)

                                                                                                                        tCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAag                              >  2:453458/1‑150 (MQ=255)

                                                                                                                          ccTGCAACAGAAGATATTCATCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAagag                            <  2:175779/150‑1 (MQ=255)

                                                                                                                           cTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGc                           <  2:1726218/150‑1 (MQ=255)

                                                                                                                              cAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGt                        >  1:1184402/1‑150 (MQ=255)

                                                                                                                                     acttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGAcct                 >  1:809620/4‑150 (MQ=255)

                                                                                                                                     acttatTCTTCAACAATTGAAGAATAAATCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCTGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGAcaa                 >  2:1242212/4‑148 (MQ=255)

                                                                                                                                      tatatTATTCAAATGTTGAATAATAAGTCGTATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCTTTTTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGAcctc                <  1:1588197/149‑1 (MQ=255)

                                                                                                                                           tcttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGt           >  1:488699/1‑150 (MQ=255)

                                                                                                                                            cttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTc          <  1:1478575/150‑1 (MQ=255)

                                                                                                                                            cttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTc          <  2:230499/150‑1 (MQ=255)

                                                                                                                                            cttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTAGTc          >  2:222309/1‑150 (MQ=255)

                                                                                                                                             ttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCt         >  1:1212727/1‑150 (MQ=255)

                                                                                                                                              tcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTCGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCtg        >  1:102549/1‑150 (MQ=255)

                                                                                                                                               cAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGACGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCtgt       <  2:323172/150‑1 (MQ=255)

                                                                                                                                                 aTTGTTGAAGAATAAGTCTGATAACAAGCTTTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCTTATTGAACGCGGAGAAATTCCGGAAACACTGCTTGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCTGTgg     <  1:463598/150‑1 (MQ=255)

                                                                                                                                                 aTTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTTCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCTGTgg     <  1:429788/150‑1 (MQ=255)

                                                                                                                                                 aTTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCTGTgg     <  2:405990/150‑1 (MQ=255)

                                                                                                                                                  ttgttgAAGAATAAGTCGGATAACAATCTGTATAAACTGACGAAAAATATTCAATGGCCTAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCTGTGGt    <  1:1153697/150‑1 (MQ=255)

                                                                                                                                                    gttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGACTGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGAGAACTCCTTGTCTGTGGtca  >  1:1271447/1‑150 (MQ=255)

GGAATGAAGAAGATGCCTTACTGGTGAATAAAGCCCTGGAGCTCGCCAGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACACTGCTGGAAATTATGCAAAAAGAGCAGTGCGACCTCCTTGTCTGTGGTCA  >  NC_000913/4113327‑4113624

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 21 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`