BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,339,805	C→T	5.9%	Q239Q (CAG→CAA)	adiA ←	arginine decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,339,805	0	C	T	5.9%	366.9 / 7.4	118	Q239Q (CAG→CAA)	adiA	arginine decarboxylase
Reads supporting (aligned to +/- strand): ref base C (58/53); new base T (3/4); total (61/57)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-02

Only 100 of 121 total aligned reads displayed.

TCCTGCGGATAGATTGGCCCGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCC  >  NC_000913/4339656‑4339951

tCCTGCGGATAGATTGGCCCGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGcc                                                                                                                                                    <  1:1416889/150‑1 (MQ=255)

tCCTGCGGATAGATTGGCCCGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGcc                                                                                                                                                    <  1:80750/150‑1 (MQ=255)

  cTGCGGATAGATTGGCCCGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTg                                                                                                                                                  >  1:921404/1‑150 (MQ=255)

    gCGGATAGATTGGCCCGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCa                                                                                                                                                <  2:1192191/150‑1 (MQ=255)

    gCGGATAGATTGGCCCGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGGCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCa                                                                                                                                                >  1:1639649/1‑150 (MQ=255)

    gCGGATAGATTGGCCCGATAAAGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCa                                                                                                                                                <  1:589844/150‑1 (MQ=255)

      ggATAGATTGGCCCGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCatg                                                                                                                                              <  1:646256/150‑1 (MQ=255)

        aTAGATTGGCCCGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAGGCTTGCatgac                                                                                                                                            <  1:294911/150‑2 (MQ=255)

          aGATTGGCCCGATACTGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAGGCTTGCatgaccg                                                                                                                                          <  1:1305933/150‑4 (MQ=255)

           gATTGGCCCGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGt                                                                                                                                         >  1:438444/1‑150 (MQ=255)

                gCCCGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGt                                                                                                                                    <  1:187147/150‑1 (MQ=255)

                gCCCGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGt                                                                                                                                    <  2:154664/150‑1 (MQ=255)

                 cccGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGtt                                                                                                                                   >  2:703706/1‑150 (MQ=255)

                 cccGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGtt                                                                                                                                   >  2:189255/1‑150 (MQ=255)

                 cccGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGGCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGtt                                                                                                                                   <  2:1639649/150‑1 (MQ=255)

                   cGATAATGCCGTAGCGGTTGCGGCTTTGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTag                                                                                                                                 >  2:307299/1‑150 (MQ=255)

                   cGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTag                                                                                                                                 <  2:886445/150‑1 (MQ=255)

                           ccGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGaa                                                                                                                         <  1:901438/150‑1 (MQ=255)

                           ccGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGaa                                                                                                                         <  2:910939/150‑1 (MQ=255)

                           ccGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGaa                                                                                                                         >  2:624058/1‑150 (MQ=255)

                           ccGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGaa                                                                                                                         <  2:780312/150‑1 (MQ=255)

                            cGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAg                                                                                                                        <  1:1386861/150‑1 (MQ=255)

                                 cGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACcg                                                                                                                   <  2:299002/150‑1 (MQ=255)

                                 cGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGGCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACcg                                                                                                                   <  2:846619/150‑1 (MQ=255)

                                   gTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACcgac                                                                                                                 >  1:1436048/1‑150 (MQ=255)

                                   gTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCGGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACcgac                                                                                                                 <  1:913499/150‑1 (MQ=255)

                                    ttGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACcgacg                                                                                                                <  2:176740/150‑1 (MQ=255)

                                    ttGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACcgacg                                                                                                                <  1:1197812/150‑1 (MQ=255)

                                      gCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCgacgac                                                                                                              <  2:938903/150‑1 (MQ=255)

                                       cGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACt                                                                                                             >  1:473199/1‑150 (MQ=255)

                                        ggCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTa                                                                                                            >  2:376240/1‑150 (MQ=255)

                                         gCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTAc                                                                                                           >  1:1370482/1‑150 (MQ=255)

                                          cTTGGCACAATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGACAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTAcc                                                                                                          >  1:870138/1‑150 (MQ=255)

                                            tttgaCCATATAGACCGGTTTCGCGCCTGGCAGCAGCATACCTTGTTCGATGGATTTATGGCAGTTACGGTAAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTGAGAGCCGGAAGTACCGACGACTACCGa                                                                                                        <  1:1064917/146‑1 (MQ=255)

                                              gCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGAcc                                                                                                      >  1:1693632/1‑150 (MQ=255)

                                                 ccATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAgg                                                                                                   <  2:901058/150‑1 (MQ=255)

                                                  cATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGa                                                                                                  <  2:921404/150‑1 (MQ=255)

                                                    tataGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGc                                                                                                <  1:1241123/150‑1 (MQ=255)

                                                        gACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATc                                                                                            >  2:578836/1‑150 (MQ=255)

                                                         aCCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCCAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCg                                                                                           <  1:190689/150‑1 (MQ=255)

                                                         aCCGGATTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCg                                                                                           <  1:147822/150‑1 (MQ=255)

                                                          ccGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAgg                                                                                                   <  1:414139/141‑1 (MQ=255)

                                                          ccGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAgg                                                                                                   >  2:414139/1‑141 (MQ=255)

                                                          ccGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCgg                                                                                          >  2:917724/1‑150 (MQ=255)

                                                          ccGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCgg                                                                                          <  2:1342221/150‑1 (MQ=255)

                                                            ggTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCa                                                                                        <  1:1667297/150‑1 (MQ=255)

                                                            ggTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTCATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCa                                                                                        <  2:892224/150‑1 (MQ=255)

                                                                   cgcCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATTATGGTGCGGTTAGAGCCGGAAGTACCTACGACTACCGACCAGGAGCGATCGGCACCAAATa                                                                                 >  1:1674423/1‑150 (MQ=255)

                                                                   cgcCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATa                                                                                 >  2:1257070/1‑150 (MQ=255)

                                                                   cgcCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACAGACGACTACCGACCAGGAGCGATCGGCACCAAATa                                                                                 <  2:524040/150‑1 (MQ=255)

                                                                      cTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGCAGCGATCGGCACCAAATAcgc                                                                              >  1:792727/1‑150 (MQ=255)

                                                                           aGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCg                                                                         >  2:1619807/1‑150 (MQ=255)

                                                                             cATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGc                                                                       >  2:188031/1‑150 (MQ=255)

                                                                              aTCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCa                                                                      <  2:339207/150‑1 (MQ=255)

                                                                                  aaCCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATAtt                                                                  >  1:1578327/1‑150 (MQ=255)

                                                                                     cTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATAttttt                                                               >  2:358882/1‑150 (MQ=255)

                                                                                     cTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATAttttt                                                               >  1:1054881/1‑150 (MQ=255)

                                                                                       ttttCGTTTGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGGATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcg                                                             <  2:379049/148‑1 (MQ=255)

                                                                                       tGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcg                                                             <  1:578836/150‑1 (MQ=255)

                                                                                           cGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcgcttt                                                         >  2:558028/1‑150 (MQ=255)

                                                                                             aTGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcgctttcg                                                       >  1:1706500/1‑150 (MQ=255)

                                                                                              tGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTtttcgctttcgc                                                      <  2:807349/150‑1 (MQ=255)

                                                                                                 aTTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCCCCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCaa                                                   >  1:1623477/1‑150 (MQ=255)

                                                                                                    tATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATg                                                >  2:295996/1‑150 (MQ=255)

                                                                                                     aTGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATgc                                               >  1:1274624/1‑150 (MQ=255)

                                                                                                      tGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATgcg                                              <  1:26201/150‑1 (MQ=255)

                                                                                                        gCAGTTACGGTCAACGACCACGACATCGGTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGcc                                            <  2:1128654/150‑1 (MQ=255)

                                                                                                          aGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAg                                          <  1:376240/150‑1 (MQ=255)

                                                                                                             tACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTAt                                       >  1:336909/1‑150 (MQ=255)

                                                                                                              aCGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATg                                      <  2:1669817/150‑1 (MQ=255)

                                                                                                               cGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATgg                                     >  2:1474823/1‑150 (MQ=255)

                                                                                                               cGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATgg                                     >  1:616406/1‑150 (MQ=255)

                                                                                                               cGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAATACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATgg                                     >  2:739575/1‑150 (MQ=255)

                                                                                                                  tCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCa                                  >  2:619395/1‑150 (MQ=255)

                                                                                                                    aaCGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAg                                <  2:336909/150‑1 (MQ=255)

                                                                                                                          cACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGa                          <  1:832656/150‑1 (MQ=255)

                                                                                                                          cACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGa                          <  1:307299/150‑1 (MQ=255)

                                                                                                                           aCGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATAGGCACCAAATACGCGAGCGGCATATTTATCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGaa                         >  1:1086224/1‑150 (MQ=255)

                                                                                                                              aCATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCg                      <  1:703706/150‑1 (MQ=255)

                                                                                                                               cATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGa                     >  2:1293828/1‑150 (MQ=255)

                                                                                                                                 tCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAgg                   <  2:1188584/150‑1 (MQ=255)

                                                                                                                                 tCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAgg                   <  2:1370482/150‑1 (MQ=255)

                                                                                                                                   gTTATCGGTCATGGAAGCCTGCATGATGGTGCGGTGAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCATATTCGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGa                 <  1:1355434/150‑1 (MQ=255)

                                                                                                                                   gTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGa                 <  2:237220/150‑1 (MQ=255)

                                                                                                                                     tATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAg               >  2:109327/1‑150 (MQ=255)

                                                                                                                                        ccatCATGCAGGCTTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAgttc            >  1:1601817/4‑150 (MQ=255)

                                                                                                                                         catCATGCAGGCTTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAgttcg           <  1:1150340/148‑1 (MQ=255)

                                                                                                                                          gTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAgttcgt          >  2:1489902/1‑150 (MQ=255)

                                                                                                                                          gTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACAGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAgttcgt          >  1:266436/1‑150 (MQ=255)

                                                                                                                                           tCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAgttcgtt         <  2:1407351/150‑1 (MQ=255)

                                                                                                                                           tCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAgttcgtt         >  1:1501754/1‑150 (MQ=255)

                                                                                                                                            cATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAgttcgttc        >  2:1580299/1‑150 (MQ=255)

                                                                                                                                              tGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGa      >  1:1377096/1‑150 (MQ=255)

                                                                                                                                                cAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATg    >  1:1248341/1‑150 (MQ=255)

                                                                                                                                                 aaGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACTACTACCGACCAGGAGCGATCGGCACAAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGc   <  2:870138/150‑1 (MQ=255)

                                                                                                                                                 aaGACTGCATGATTGTGCGGTTAGAGCCGGAAGTACCTACGAATACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGc   <  1:122290/150‑1 (MQ=255)

                                                                                                                                                  cgCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGAAGACTACCGACCAGGAGCGATCGGCACCAAATTCGCGTGCGGCATATTTTTCGATTTCGCCAATTGCGGTAGTCTGGTCAAGCAAAGTGCCAAGGGAAGTTCGTACGCTGcc  >  2:1582107/2‑150 (MQ=255)

                                                                                                                                                  aGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGcc  >  2:645658/1‑150 (MQ=255)

                                                                                                                                                  aGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGcc  <  1:109327/150‑1 (MQ=255)

                                                                                                                                                  aGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGcc  <  1:358882/150‑1 (MQ=255)

TCCTGCGGATAGATTGGCCCGATAATGCCGTAGCGGTTGCGGCTTGGCACCATATAGACCGGTTTCGCGCCTGTCAGCATCAAACCTTGTTCGATGGATTTATGGCAGTTACGGTCAACGACCACGACATCGTTATCGGTCATGCAAGCCTGCATGATGGTGCGGTTAGAGCCGGAAGTACCGACGACTACCGACCAGGAGCGATCGGCACCAAATACGCGTGCGGCATATTTTTCGCTTTCGCCAAATGCGCCAGTATGGTCAAGCAAAGAACCGAGGGAAGTTCGTTCGATGCC  >  NC_000913/4339656‑4339951

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