BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,242,000	T→A	5.0%	*158L (TAA→TTA)	bfr‑II ←	bacterioferritin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,242,000	0	T	A	5.0%	429.7 / 6.1	127	*158L (TAA→TTA)	bfr‑II	bacterioferritin
Reads supporting (aligned to +/- strand): ref base T (63/57); new base A (5/2); total (68/59)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

AGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCACGCGAGCCGAAGTCGCCAGCGGTTTCA  >  NC_002947/1241852‑1242149

aGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTa                                                                                                                                                      <  1:568427/150‑1 (MQ=255)

 gCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTAc                                                                                                                                                     <  2:2465370/150‑1 (MQ=255)

  cGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTAct                                                                                                                                                    <  2:1187241/150‑1 (MQ=255)

  cGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTAct                                                                                                                                                    >  2:1618018/1‑150 (MQ=255)

      atcatGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTc                                                                                                                                                <  1:3884385/150‑1 (MQ=255)

      atcatGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTc                                                                                                                                                >  1:4939599/1‑150 (MQ=255)

      atcatGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCCGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTc                                                                                                                                                <  1:2643466/150‑1 (MQ=255)

           gCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCa                                                                                                                                           >  1:4072469/1‑150 (MQ=255)

             tGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATc                                                                                                                                         >  2:4301090/1‑150 (MQ=255)

               gAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAATCCCGC‑GATCAGCGGGTTTTTTTATGGGGCGACTAGTTACTCTTCGCCCATCTg                                                                                                                                       >  1:4432956/1‑150 (MQ=255)

               gAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTGCGCCCATCTg                                                                                                                                       >  2:2900617/1‑150 (MQ=255)

                aaaTGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCCAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTGCGCCCATCTGc                                                                                                                                      >  1:3088701/1‑150 (MQ=255)

                aaaTGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGc                                                                                                                                      <  2:1785187/150‑1 (MQ=255)

                aaaTGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGc                                                                                                                                      <  2:4072469/150‑1 (MQ=255)

                    ggtCGAGAAGGGTGAGCTGCAAAACCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAATAACCCGC‑GAGCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGAca                                                                                                                                  >  1:1449581/3‑149 (MQ=255)

                      tCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGc                                                                                                                                <  1:4086799/150‑1 (MQ=255)

                       cGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATACAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCa                                                                                                                               <  1:3774557/150‑1 (MQ=255)

                        gagaAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAg                                                                                                                              >  1:1631666/1‑150 (MQ=255)

                             gggTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCCTCGCCCATCTGCGACTGCAGGTAGt                                                                                                                         >  2:2644147/1‑150 (MQ=255)

                                tGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAAGAGTAACTATTCGcccc                                                                                                                                           <  1:964485/129‑2 (MQ=255)

                                tGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAAGAGTAACTATTCGcccc                                                                                                                                           >  2:964485/1‑128 (MQ=255)

                                 gAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTtctc                                                                                                                     <  2:543710/150‑1 (MQ=255)

                                 gAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTtctc                                                                                                                     <  2:3350861/150‑1 (MQ=255)

                                   gCTGCAATCCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGa                                                                                                                   >  1:3202284/1‑150 (MQ=255)

                                   gCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGa                                                                                                                   <  1:3330612/150‑1 (MQ=255)

                                   gCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTGTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGa                                                                                                                   >  2:3762673/1‑150 (MQ=255)

                                    cTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGAt                                                                                                                  >  1:3175561/1‑150 (MQ=255)

                                        aaaCCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCg                                                                                                              <  2:3236421/150‑1 (MQ=255)

                                           cccTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATc                                                                                                           <  2:3591592/150‑1 (MQ=255)

                                            ccTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCt                                                                                                          >  1:2846087/1‑150 (MQ=255)

                                               tgctggggACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGt                                                                                                       <  2:2695951/144‑1 (MQ=255)

                                                aTCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTc                                                                                                      >  1:2979438/1‑150 (MQ=255)

                                                aTCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTc                                                                                                      <  1:4155881/150‑1 (MQ=255)

                                                  cAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGa                                                                                                    <  1:5034638/150‑1 (MQ=255)

                                                  cAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGa                                                                                                    <  1:5034546/150‑1 (MQ=255)

                                                   aaGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGAt                                                                                                   >  1:2237552/1‑150 (MQ=255)

                                                     ggACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCa                                                                                                 <  1:4301090/150‑1 (MQ=255)

                                                     ggACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCa                                                                                                 <  1:1845329/150‑1 (MQ=255)

                                                       aCGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGc                                                                                               <  2:2237552/150‑1 (MQ=255)

                                                        cgcgGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGcc                                                                                              <  1:4054782/150‑1 (MQ=255)

                                                           ggCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCa                                                                                           <  2:126643/150‑1 (MQ=255)

                                                              tGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGtt                                                                                        >  1:2661861/1‑150 (MQ=255)

                                                              tGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGtt                                                                                        <  2:3998034/150‑1 (MQ=255)

                                                                ccaaggccaaggCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTgg                                                                                      <  1:875806/150‑1 (MQ=255)

                                                                 caaggccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTggg                                                                                     >  2:3631481/1‑150 (MQ=255)

                                                                 caaggccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTggg                                                                                     <  2:1486172/150‑1 (MQ=255)

                                                                 caaggccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTggg                                                                                     >  1:866449/1‑150 (MQ=255)

                                                                   aggccaaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGtt                                                                                   >  1:1926279/1‑150 (MQ=255)

                                                                       caaggcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGACTAGTTACTCTTCGCCCATCTGCGCCTGCAGGTAGTTCTCGTTACCGATCTTGTCGATCAGCCCCAGTTGGGTTTcca                                                                               >  2:3638319/1‑150 (MQ=255)

                                                                           gcTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGAAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTccagccc                                                                           >  2:2125740/1‑149 (MQ=255)

                                                                              tCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTc                                                                        <  1:4745190/150‑1 (MQ=255)

                                                                                 gAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGAt                                                                     <  2:5163696/150‑1 (MQ=255)

                                                                                 gAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGAt                                                                     <  2:3322991/150‑1 (MQ=255)

                                                                                 gAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGAt                                                                     >  1:5163696/1‑150 (MQ=255)

                                                                                   aGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGAtgt                                                                   <  1:546257/150‑1 (MQ=255)

                                                                                    gCCTGCTCGTCCGGCACGCCTAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGAtgtg                                                                  <  1:2899351/150‑1 (MQ=255)

                                                                                     ccTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGc                                                                 >  1:647992/1‑150 (MQ=255)

                                                                                       tGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGctc                                                               >  1:3754180/1‑150 (MQ=255)

                                                                                          tCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTc                                                            >  1:1563065/1‑150 (MQ=255)

                                                                                                ggCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGa                                                      >  2:2501846/1‑150 (MQ=255)

                                                                                                    agaaTAAAAAAACCCGC‑GATCAGCGGTTTTTTTTAT‑GGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGa                                                 <  2:736866/147‑1 (MQ=255)

                                                                                                      cATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGACAGTCGATGTGCTCTTCCTCGGACTCGAg                                                >  2:3664210/1‑150 (MQ=255)

                                                                                                       aTAACAAAACCCGC‑GAGCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGCTCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAgg                                               <  2:4432956/150‑1 (MQ=255)

                                                                                                       aTAAAAAAACCCGCTGATC‑GCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGGTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAgg                                               >  1:5389211/1‑150 (MQ=255)

                                                                                                       aTAAAAAAACCCGC‑GATCAGCGGGTTTTTTTTTGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCATGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAgg                                               <  1:3344665/150‑1 (MQ=255)

                                                                                                        tAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTGTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGa                                              >  2:2981245/1‑150 (MQ=255)

                                                                                                          aaaaaaCCCGCTGATC‑GCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCtt                                                                                                         <  2:1517609/89‑1 (MQ=38)

                                                                                                          aaaaaaCCCGCTGATC‑GCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCtt                                                                                                         >  1:1517609/1‑89 (MQ=38)

                                                                                                          aaaaaaCCCGCTGATC‑GCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGata                                            >  1:1644263/1‑150 (MQ=255)

                                                                                                          aaaaaaCCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGata                                            >  1:4732004/1‑150 (MQ=255)

                                                                                                           aaaaaCCCGCTGATC‑GCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGc                                                                                               >  1:3702165/1‑98 (MQ=38)

                                                                                                           aaaaaCCCGCTGATC‑GCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGc                                                                                               <  2:3702165/98‑1 (MQ=38)

                                                                                                              aaCCCGCTGATC‑GCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCtt                                        >  2:3793781/1‑150 (MQ=255)

                                                                                                              aaCCCGCTGATC‑GCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCtt                                        >  1:1008030/1‑150 (MQ=255)

                                                                                                              aaCCCGCTGATC‑GCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCtt                                        >  1:1007938/1‑150 (MQ=255)

                                                                                                               aCCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTc                                       <  2:3991928/150‑1 (MQ=255)

                                                                                                               aCCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTc                                       <  2:1693248/150‑1 (MQ=255)

                                                                                                                 ccGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGa                                     >  1:1021487/1‑150 (MQ=255)

                                                                                                                 ccGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGa                                     <  2:1344062/150‑1 (MQ=255)

                                                                                                                   gc‑gATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGc                                   >  1:5094627/1‑150 (MQ=255)

                                                                                                                           gCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCAc                            >  1:2995654/1‑150 (MQ=255)

                                                                                                                           gCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCAc                            >  2:939841/1‑150 (MQ=255)

                                                                                                                           gCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCAc                            >  1:1209979/1‑150 (MQ=255)

                                                                                                                            cGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAgg                                               >  1:3325498/1‑130 (MQ=255)

                                                                                                                            cGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAgg                                               <  2:3325498/130‑1 (MQ=255)

                                                                                                                            cGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCAcg                           >  1:2410223/1‑150 (MQ=255)

                                                                                                                             gggTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCAcgc                          >  2:295548/1‑150 (MQ=255)

                                                                                                                             gggTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCAcgc                          >  1:529815/1‑150 (MQ=255)

                                                                                                                             gggTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCAcgc                          >  1:529794/1‑150 (MQ=255)

                                                                                                                               gTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCACGCGa                        >  2:1055867/1‑150 (MQ=255)

                                                                                                                                      agatgggCGAAGAGTAACTATTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCACGCGAGCCGAAg                 >  1:5253342/5‑150 (MQ=255)

                                                                                                                                         tgggCGAAGAGTAACTATTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCACGCGAGCCGAAGTCg              >  2:3299575/2‑150 (MQ=255)

                                                                                                                                         tgggCGAAGAGTAACTATTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCACGCGAGCCGAAGTCg              >  2:3299484/2‑150 (MQ=255)

                                                                                                                                          gggCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCACGCGAGCCGAAGTCGc             <  1:3223368/150‑1 (MQ=255)

                                                                                                                                           ggCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCACGCGAGCCGAAGTCGcc            >  2:4618300/1‑150 (MQ=255)

                                                                                                                                            gCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCACGCGAGCCGAAGTCGCCa           <  2:2661861/150‑1 (MQ=255)

                                                                                                                                                aTAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCACGCGAGCCGAAGTCGCCAGCgg       >  1:976963/1‑150 (MQ=255)

                                                                                                                                                    ttACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCACGCGAGCCGAAGTCGCCATCgcccat   >  2:4977668/1‑145 (MQ=255)

                                                                                                                                                    ttACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCACGCGAGCCGAAGTCGCCAGCGGTTTc   >  2:1235728/1‑150 (MQ=255)

                                                                                                                                                     tACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCACGCGAGCCGAAGTCGCCAGCGGTTTca  <  2:2937344/150‑1 (MQ=255)

AGCGACATCATGCTGGAAATGTTCGAGAAGGGTGAGCTGCAAACCCTGATCAAGGACGCGGCTGCCAAGGCCAAGGCTTCTGAAGCCTGATCTTCCGGCAAGCATAAAAAAACCCGC‑GATCAGCGGGTTTTTTTATGGGGCGAATAGTTACTCTTCGCCCATCTGCGACTGCAGGTAGTTCTCGATACCGATCTTGTCGATCAGCCCCAGTTGGGTTTCCAGCCAGTCGATGTGCTCTTCCTCGGACTCGAGGATATCTTCGAGCACGTCACGCGAGCCGAAGTCGCCAGCGGTTTCA  >  NC_002947/1241852‑1242149

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