BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,447,964	T→A	5.9%	L24* (TTA→TAA)	ytfP →	GGCT‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,447,964	0	T	A	5.9%	591.9 / 18.9	170	L24* (TTA→TAA)	ytfP	GGCT‑like protein
Reads supporting (aligned to +/- strand): ref base T (83/77); new base A (4/6); total (87/83)
Fisher's exact test for biased strand distribution p-value = 5.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 171 total aligned reads displayed.

ATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGG  >  NC_000913/4447817‑4448111

aTGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTa                                                                                                                                                    >  1:712354/1‑149 (MQ=255)

   ccTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTg                                                                                                                                                 <  3:500151/149‑1 (MQ=255)

     tAAGTTATATCTGGAAGCCGTGTCGGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTggg                                                                                                                                               >  1:938326/1‑149 (MQ=255)

     tAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTGTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTggg                                                                                                                                               >  2:943489/1‑149 (MQ=255)

     tAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTggg                                                                                                                                               >  2:508563/1‑149 (MQ=255)

       ggtcctatTTGGAACCCGTGTCTGGGGAAGCCAAGGCCTTGGATTGGCTCTTTCAGTCCAAGTTTAGGCAATCGGAATATTTGCCACCGCCAG‑TTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGa                                                                                                                                            <  4:46354/144‑1 (MQ=255)

        gCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTAACTGGGcat                                                                                                                                            <  3:1075048/149‑3 (MQ=255)

        gCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTAACTGGGcat                                                                                                                                            <  2:1044616/149‑3 (MQ=255)

         cTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAt                                                                                                                                           >  2:770407/1‑149 (MQ=255)

         cTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAt                                                                                                                                           <  1:679912/149‑1 (MQ=255)

             atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCa                                                                                                                                       <  4:138122/149‑1 (MQ=255)

               cTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTAAGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGt                                                                                                                                     <  3:1021695/149‑1 (MQ=255)

                 ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                                                                                                                                   >  3:443446/1‑149 (MQ=255)

                 ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                                                                                                                                   >  3:499590/1‑149 (MQ=255)

                  gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATc                                                                                                                                  <  2:22422/149‑1 (MQ=255)

                    aGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGa                                                                                                                                <  4:838844/149‑1 (MQ=255)

                     gCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGAt                                                                                                                               <  3:693588/149‑1 (MQ=255)

                     gCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGAt                                                                                                                               <  3:607152/149‑1 (MQ=255)

                      ccGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATa                                                                                                                              >  3:843748/1‑149 (MQ=255)

                       cGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATaa                                                                                                                             <  1:396795/149‑1 (MQ=255)

                           tCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTAc                                                                                                                         <  2:449506/149‑1 (MQ=255)

                             tGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTCCGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTACCGATAACTACCa                                                                                                                       >  3:202393/1‑149 (MQ=255)

                             tGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGGTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATGGATAACTACCa                                                                                                                       >  4:118855/1‑149 (MQ=255)

                              ggTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAg                                                                                                                      <  3:464544/149‑1 (MQ=255)

                                tgtAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTGCCAGtt                                                                                                                    <  1:138587/149‑1 (MQ=255)

                                tgtAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGtt                                                                                                                    >  1:743153/1‑149 (MQ=255)

                                tgtAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGtt                                                                                                                    >  3:293832/1‑149 (MQ=255)

                                     aCCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGtata                                                                                                               >  2:35687/1‑149 (MQ=255)

                                     aCCAGACACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGtata                                                                                                               <  4:364616/149‑1 (MQ=255)

                                        aGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGcc                                                                                                            >  2:1104765/1‑149 (MQ=255)

                                        aGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTCGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGcc                                                                                                            <  3:889942/149‑1 (MQ=255)

                                          gCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTg                                                                                                          <  1:1084469/149‑1 (MQ=255)

                                                gATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTAc                                                                                                                         >  4:683294/1‑128 (MQ=255)

                                                  tttGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTa                                                                                                  >  1:334637/1‑149 (MQ=255)

                                                     gCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTgg                                                                                                         <  3:809040/139‑1 (MQ=255)

                                                     gCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATcc                                                                                               >  3:30358/1‑149 (MQ=255)

                                                       tctATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAg                                                                                             >  3:107464/1‑149 (MQ=255)

                                                         tATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGgc                                                                                           <  4:99070/149‑1 (MQ=255)

                                                           tCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGgcgc                                                                                         >  1:768310/1‑149 (MQ=255)

                                                           tCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGgcgc                                                                                         >  1:838695/1‑149 (MQ=255)

                                                              gTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGt                                                                                      >  4:1051715/1‑149 (MQ=255)

                                                              gTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTGTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAGCTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGt                                                                                      >  1:757766/1‑149 (MQ=255)

                                                               ttCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGtt                                                                                     >  2:556344/1‑149 (MQ=255)

                                                               ttCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGtt                                                                                     >  4:634513/1‑149 (MQ=255)

                                                               ttCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGtt                                                                                     >  4:347081/1‑149 (MQ=255)

                                                                 cGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTcc                                                                                   >  2:32470/1‑149 (MQ=255)

                                                                   aGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCgg                                                                                 <  2:768310/149‑1 (MQ=255)

                                                                   aGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCgg                                                                                 <  2:569855/149‑1 (MQ=255)

                                                                     ttttAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTgg                                                                                                                                                <  4:928017/84‑1 (MQ=255)

                                                                         aGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGggaacg                                                                           >  3:504477/1‑149 (MQ=255)

                                                                          gCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGggaacgg                                                                          >  2:185551/1‑149 (MQ=255)

                                                                            aaTGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGggaacggaa                                                                        <  2:655292/149‑1 (MQ=255)

                                                                            aaTGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGggaacggaa                                                                        >  1:321690/1‑149 (MQ=255)

                                                                              tGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTTTCGATAACTACCAGTTGTATCGCCAGGGCCACTATCCAGACGCAGTTCCTGtgaacgcatcg                                                                      <  4:901520/149‑6 (MQ=255)

                                                                                  aaTATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTac                                                                  >  3:932263/1‑149 (MQ=255)

                                                                                     atTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTc                                                                                    >  2:83675/1‑128 (MQ=255)

                                                                                      tttGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACgg                                                              >  1:472145/1‑149 (MQ=255)

                                                                                        tGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTg                                                            <  3:393914/149‑1 (MQ=255)

                                                                                          tCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGaa                                                          <  3:41320/149‑1 (MQ=255)

                                                                                          tCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGaa                                                          <  1:454883/149‑1 (MQ=255)

                                                                                           cTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCGGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAg                                                         <  2:1027591/149‑1 (MQ=255)

                                                                                                gCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTAt                                                    <  2:458170/149‑1 (MQ=255)

                                                                                                gCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTAt                                                    <  2:707411/149‑1 (MQ=255)

                                                                                                 cAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATc                                                   <  3:326632/149‑1 (MQ=255)

                                                                                                     ttACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTAt                                               <  2:360159/149‑1 (MQ=255)

                                                                                                        cGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGa                                            <  4:749554/149‑1 (MQ=255)

                                                                                                          ccACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACa                                          >  4:340320/1‑149 (MQ=255)

                                                                                                           cacaAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGa                                            <  3:221501/146‑1 (MQ=255)

                                                                                                           cacaAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGa                                            >  4:221501/1‑146 (MQ=255)

                                                                                                           cacaAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTAACTGGGcat                                                                                                                                            <  2:35560/50‑3 (MQ=255)

                                                                                                           cacaAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTAACTGGGcat                                                                                                                                            >  1:35560/1‑48 (MQ=255)

                                                                                                                 caaGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCAc                                   <  1:147533/149‑1 (MQ=255)

                                                                                                                  aaGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACg                                  >  2:1061135/1‑149 (MQ=255)

                                                                                                                     gCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTg                               >  2:659387/1‑149 (MQ=255)

                                                                                                                           gTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGaa                         <  1:1082655/149‑1 (MQ=255)

                                                                                                                              aCTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTAc                                                                                                                         <  3:895930/50‑1 (MQ=255)

                                                                                                                              aCTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTAc                                                                                                                         >  4:895930/1‑50 (MQ=255)

                                                                                                                              aCTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACtt                      >  2:306850/1‑149 (MQ=255)

                                                                                                                                tGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGa                    >  3:496227/1‑149 (MQ=255)

                                                                                                                                 ggATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGAt                   <  2:514385/149‑1 (MQ=255)

                                                                                                                                 ggATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGAt                   <  3:579985/149‑1 (MQ=255)

                                                                                                                                 ggATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGAt                   <  1:624800/149‑1 (MQ=255)

                                                                                                                                      aCCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCtt              <  4:30358/149‑1 (MQ=255)

                                                                                                                                        cAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTgc            >  3:383415/1‑149 (MQ=255)

                                                                                                                                        cAATGCCCAGTTACGGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGa                    >  1:827415/1‑141 (MQ=255)

                                                                                                                                         atcgCCCAGTAACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTgcg           >  3:863941/4‑149 (MQ=255)

                                                                                                                                         aaTGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTgcg           >  4:164867/1‑149 (MQ=255)

                                                                                                                                          tcgCCCAGTAACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTgg                                                                                                         <  3:923801/52‑1 (MQ=255)

                                                                                                                                          tcgCCCAGTAACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGAt                   <  1:865680/138‑1 (MQ=255)

                                                                                                                                          tcgCCCAGTAACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGAt                   >  2:865680/3‑140 (MQ=255)

                                                                                                                                          aTGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGggaac                                                                            >  1:552558/1‑83 (MQ=255)

                                                                                                                                            gCCCAGTTACTGGGCGACTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCAc        >  1:117892/1‑149 (MQ=255)

                                                                                                                                              ccAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGa                    <  4:419684/135‑1 (MQ=255)

                                                                                                                                              ccAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCtt              >  4:1050757/1‑141 (MQ=255)

                                                                                                                                              ccAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCtt              <  3:1050757/141‑1 (MQ=255)

                                                                                                                                              ccAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCa      >  3:452632/1‑149 (MQ=255)

                                                                                                                                              ccAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCa      >  1:181854/1‑149 (MQ=255)

                                                                                                                                                aGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCt               >  4:504327/1‑138 (MQ=255)

                                                                                                                                                 gTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAggg   >  3:681539/1‑149 (MQ=255)

                                                                                                                                                  ttACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAgggg  <  2:528187/149‑1 (MQ=255)

ATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGG  >  NC_000913/4447817‑4448111

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: `‑`