BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,995,724	C→T	5.2%	R48C (CGC→TGC)	yigA →	conserved protein YigA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,995,724	0	C	T	5.2%	305.1 / 4.0	95	R48C (CGC→TGC)	yigA	conserved protein YigA
Reads supporting (aligned to +/- strand): ref base C (46/44); new base T (3/2); total (49/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

Only 100 of 101 total aligned reads displayed.

TTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACC  >  NC_000913/3995576‑3995873

ttCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTAcg                                                                                                                                                      >  1:1611517/1‑150 (MQ=255)

 tCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACTTAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTAcgc                                                                                                                                                     <  1:598230/150‑1 (MQ=255)

  cATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTAcgc                                                                                                                                                     <  2:1612620/149‑1 (MQ=255)

  cATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTAcgc                                                                                                                                                     >  1:1612620/1‑149 (MQ=255)

    tgtATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGc                                                                                                                                                  <  2:508593/148‑1 (MQ=255)

    tCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGc                                                                                                                                                  <  2:1611517/150‑1 (MQ=255)

     cTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCa                                                                                                                                                 >  1:967967/1‑150 (MQ=255)

         gAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCGTACCGGATGCGGCAcacg                                                                                                                                             <  2:448331/150‑4 (MQ=255)

         gAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCGTACCGGATGCGGCAcacg                                                                                                                                             <  2:1226584/150‑4 (MQ=255)

           aGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGttt                                                                                                                                           >  1:491425/1‑150 (MQ=255)

           aGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGttt                                                                                                                                           >  1:334283/1‑150 (MQ=255)

            gCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTc                                                                                                                                          >  1:1544999/1‑150 (MQ=255)

            gCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTc                                                                                                                                          <  1:1412987/150‑1 (MQ=255)

            gCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTc                                                                                                                                          <  1:29506/150‑1 (MQ=255)

               aCCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGtt                                                                                                                                       <  1:1516274/150‑1 (MQ=255)

               aCCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGtt                                                                                                                                       >  1:81261/1‑150 (MQ=255)

                  aGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGt                                                                                                                                    <  1:1471781/150‑1 (MQ=255)

                   ggggAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTTCCGCATCCGGTACGCTGAACCGTTTCGTTGTTc                                                                                                                                   <  2:961700/150‑1 (MQ=255)

                   ggggAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTc                                                                                                                                   <  1:1090285/150‑1 (MQ=255)

                    gggAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCg                                                                                                                                  <  2:447936/150‑1 (MQ=255)

                          aaCTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGc                                                                                                                            <  1:1042219/150‑1 (MQ=255)

                           aCTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGca                                                                                                                           <  2:1557073/150‑1 (MQ=255)

                            cTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGcaa                                                                                                                          >  1:618468/1‑149 (MQ=255)

                              gCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGATGGTCGAGTGGCACAt                                                                                                                        >  1:706824/1‑150 (MQ=255)

                                      cacTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCcgcg                                                                                                                >  2:231497/1‑150 (MQ=255)

                                      cacTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCcgcg                                                                                                                >  2:885927/1‑150 (MQ=255)

                                      cacTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCcgcg                                                                                                                <  2:789185/150‑1 (MQ=255)

                                        cTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCa                                                                                                              <  1:308053/150‑1 (MQ=255)

                                        cTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCa                                                                                                              <  2:1282127/150‑1 (MQ=255)

                                         tCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCAc                                                                                                             >  1:1012172/1‑150 (MQ=255)

                                          cACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACg                                                                                                            >  1:1524712/1‑150 (MQ=255)

                                                gCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCa                                                                                                      <  2:334283/150‑1 (MQ=255)

                                                gCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCa                                                                                                      <  2:276989/150‑1 (MQ=255)

                                                gCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCa                                                                                                      >  2:1209626/1‑150 (MQ=255)

                                                 cTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGAGCGCGCAGTAGAAGCGATACGTGTGCCGAATCCGGTACGCGGCAACGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCaa                                                                                                     >  2:160397/1‑149 (MQ=255)

                                                   tgatgaCCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCatat                                                                                                   >  2:1164199/1‑150 (MQ=255)

                                                   tgatgaCCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCAGGTACGAGGCACCGTTTAGTTTGTATAGTGGCACATGGACCGCGCACGTAATCatat                                                                                                   >  1:553139/1‑150 (MQ=255)

                                                        aCCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATg                                                                                              >  1:1645582/1‑150 (MQ=255)

                                                        aCCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATg                                                                                              <  1:231497/150‑1 (MQ=255)

                                                         ccGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGt                                                                                             >  1:943995/1‑150 (MQ=255)

                                                          cGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGtt                                                                                            >  2:830389/1‑150 (MQ=255)

                                                                gTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGaa                                                                                      >  2:95069/1‑150 (MQ=255)

                                                                  tGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAaga                                                                                    <  2:967967/150‑1 (MQ=255)

                                                                    tCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGGGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTTATCATATTCATGTTCTGGAAgaga                                                                                  >  1:714681/1‑150 (MQ=255)

                                                                       aTTATCTGATTATAAATCCTGAGTTTTTTTTCCTTAATGCGGGCGCAGTAGAAGCGATACGTGTGCGGCATCCGGTACGCGGCACCGTTTCGTTGGTCGTGTGGCAGATGGCCGGCGCACGTAATCATATTCATGTGCTGGTAGAGAACa                                                                               <  2:849725/150‑1 (MQ=255)

                                                                       aTTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACa                                                                               <  1:248794/150‑1 (MQ=255)

                                                                         tATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATg                                                                             >  1:1207729/1‑150 (MQ=255)

                                                                           tCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGgc                                                                           >  2:602941/1‑150 (MQ=255)

                                                                                 tAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGtt                                                                     <  2:1326304/150‑1 (MQ=255)

                                                                                  aaaaaTCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTg                                                                    <  1:219364/150‑1 (MQ=255)

                                                                                   aaaaTCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGa                                                                   >  1:1631678/1‑150 (MQ=255)

                                                                                        ccTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGaa                                                              >  2:1150499/1‑150 (MQ=255)

                                                                                         cTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAAc                                                             >  2:857670/1‑150 (MQ=255)

                                                                                            aGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAgg                                                          <  1:683011/150‑1 (MQ=255)

                                                                                             gTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGTCACCGTTTCGTTGGGCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACTGGc                                                         <  2:1544999/150‑1 (MQ=255)

                                                                                             gTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGc                                                         <  2:943995/150‑1 (MQ=255)

                                                                                             gTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGc                                                         <  2:831197/150‑1 (MQ=255)

                                                                                                ttttATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCt                                                                                          >  1:21505/1‑114 (MQ=255)

                                                                                                ttttATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCt                                                                                          <  2:21505/114‑1 (MQ=255)

                                                                                                   tATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGc                                                   <  1:95069/150‑1 (MQ=255)

                                                                                                      ccGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCaa                                                >  2:1517906/1‑150 (MQ=255)

                                                                                                        gTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACg                                              >  1:441947/1‑150 (MQ=255)

                                                                                                            tGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCTCCAACGAAgg                                          <  2:618468/150‑1 (MQ=255)

                                                                                                            tGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAgg                                          <  1:538220/150‑1 (MQ=255)

                                                                                                            tGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAgg                                          >  1:981782/1‑150 (MQ=255)

                                                                                                                  cgcAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGtt                                    >  2:645903/1‑150 (MQ=255)

                                                                                                                  cgcAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGtt                                    <  2:81261/150‑1 (MQ=255)

                                                                                                                   gcAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTTGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGttt                                   <  2:491425/150‑1 (MQ=255)

                                                                                                                   gcAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGttt                                   <  1:1517906/150‑1 (MQ=255)

                                                                                                                     aGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTa                                 >  2:212659/1‑150 (MQ=255)

                                                                                                                     aGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTa                                 <  1:602941/150‑1 (MQ=255)

                                                                                                                     aGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTa                                 <  1:1447407/150‑1 (MQ=255)

                                                                                                                      gTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTAt                                >  2:513814/1‑150 (MQ=255)

                                                                                                                          aaGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGcc                            >  2:32952/1‑150 (MQ=255)

                                                                                                                           aGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCt                           >  1:142185/1‑150 (MQ=255)

                                                                                                                            gCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTa                          >  1:1596378/1‑150 (MQ=255)

                                                                                                                             cGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTAc                         <  2:1207729/150‑1 (MQ=255)

                                                                                                                              gATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAAc                                                                                <  1:1666260/94‑1 (MQ=255)

                                                                                                                              gATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAAc                                                                                >  2:1666260/1‑94 (MQ=255)

                                                                                                                                tACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTActct                      >  2:1078387/1‑150 (MQ=255)

                                                                                                                                  acggtgCCGCGTACCGGATGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTAc                    >  1:635892/4‑150 (MQ=255)

                                                                                                                                   gtgtGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGTGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTAcc                   <  1:1352358/150‑1 (MQ=255)

                                                                                                                                   cggtgCCGCGTACCGGATGCGGCACCGTTTCGTTCGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTAcc                   >  2:331222/3‑150 (MQ=255)

                                                                                                                                    tgtgCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCt                  >  2:1024755/1‑150 (MQ=255)

                                                                                                                                    ggtgCCGCGTACCGGATGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCt                  >  1:1203999/2‑150 (MQ=255)

                                                                                                                                     gtgCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTg                 >  2:1651402/1‑150 (MQ=255)

                                                                                                                                     gtgCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTg                 >  1:408057/1‑150 (MQ=255)

                                                                                                                                       gCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCa               >  1:59348/1‑150 (MQ=255)

                                                                                                                                        ccGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTActct                      <  2:30395/142‑1 (MQ=255)

                                                                                                                                        ccGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTActct                      >  1:30395/1‑142 (MQ=255)

                                                                                                                                        ccGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAg              >  1:914140/1‑150 (MQ=255)

                                                                                                                                           cATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAgcgc           <  2:59348/150‑1 (MQ=255)

                                                                                                                                             tCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAg         <  1:1651402/150‑1 (MQ=255)

                                                                                                                                             tCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAg         <  2:1145155/150‑1 (MQ=255)

                                                                                                                                             tCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCAAGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAg         <  2:519426/150‑1 (MQ=255)

                                                                                                                                              ccGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGt        >  2:816871/1‑150 (MQ=255)

                                                                                                                                               cGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGtc       >  1:589108/1‑150 (MQ=255)

                                                                                                                                                   aCGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCAc   >  1:1623474/1‑150 (MQ=255)

                                                                                                                                                    cgcgGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCATACTCTACCTGCAGCGCAGTCTCAcc  >  2:719834/1‑150 (MQ=255)

                                                                                                                                                    cgcgGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGACTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCAcc  >  2:1403194/1‑150 (MQ=255)

TTCATCTATGAAGCAACCAGGGGAAGAACTGCAGGAAACACTCACGGAGCTTGATGACCGGGCGGTTGTCGATTATCTGATTAAAAATCCTGAGTTTTTTATCCGTAATGCGCGCGCAGTAGAAGCGATACGTGTGCCGCATCCGGTACGCGGCACCGTTTCGTTGGTCGAGTGGCACATGGCCCGCGCACGTAATCATATTCATGTTCTGGAAGAGAACATGGCGCTGTTGATGGAACAGGCTATCGCCAACGAAGGCCTGTTTTATCGCCTACTCTACCTGCAGCGCAGTCTCACC  >  NC_000913/3995576‑3995873

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