BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,247,843:1	+G	5.9%	coding (71/663 nt)	yqjA →	DedA family protein YqjA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,247,843	1	.	G	5.9%	624.5 / 35.4	152	coding (71/663 nt)	yqjA	DedA family protein YqjA
Reads supporting (aligned to +/- strand): ref base . (72/71); new base G (2/7); total (74/78)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

Only 100 of 162 total aligned reads displayed.

ATACCAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTG  >  NC_000913/3247694‑3247991

 tACCAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCgat‑g                                                                                                                                                     <  2:1474480/150‑1 (MQ=255)

   ccAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCgat‑gat                                                                                                                                                   >  1:2309394/1‑150 (MQ=255)

   ccAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCgat‑gat                                                                                                                                                   >  1:433996/1‑150 (MQ=255)

    cAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GAtt                                                                                                                                                  >  1:610306/1‑150 (MQ=255)

     aaaaCAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTg                                                                                                                                                 >  1:1835280/1‑150 (MQ=255)

       aaCAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAAT‑CATCGATGGATTGGc                                                                                                                                               <  2:2389833/150‑1 (MQ=255)

        aCAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCa                                                                                                                                              <  2:446212/150‑1 (MQ=255)

         cAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAAT‑CATCGATGGATTGGcca                                                                                                                                             <  1:862952/150‑3 (MQ=255)

          aGCATTAACATTGCTGTAGTGAATCTTTTAAAACCAGCGTTGCGTTAACCGATTCACCAGTAATAATGAAGTTAACTTTTGACCCAATTGATGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAAT‑CATCGATGGATTGGccat                                                                                                                                            <  2:1343876/150‑4 (MQ=255)

           gCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAAT‑CATCGATGGATTGGccagg                                                                                                                                           <  2:206139/150‑5 (MQ=255)

           gCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCCGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAAT‑CATCGATGGATTGGccagg                                                                                                                                           <  2:1389855/150‑5 (MQ=255)

           ccAAAAACATTTCCGTACGGCATCCTTAAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAAT‑CATCGATGGATTGGccagg                                                                                                                                           <  1:380076/149‑5 (MQ=255)

                 aCATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTAtt                                                                                                                                     <  1:352238/150‑1 (MQ=255)

                  cATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTTGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTAttt                                                                                                                                    >  1:2272488/1‑150 (MQ=255)

                  cATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTAttt                                                                                                                                    <  1:2031714/150‑1 (MQ=255)

                   tttgggttAGTGATTCTTTTAAACCCTG‑GTGGGGTTAAACGTTTTAACTTGTATAATTTAATGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTg                                                                                                                                  <  2:2252238/143‑1 (MQ=255)

                    ttGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTg                                                                                                                                  <  2:1445472/150‑1 (MQ=255)

                     tGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGt                                                                                                                                 <  2:2456585/150‑1 (MQ=255)

                       cGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCt                                                                                                                               <  1:1735188/150‑1 (MQ=255)

                       cGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCt                                                                                                                               <  1:1563566/150‑1 (MQ=255)

                         gTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCttg                                                                                                                             >  1:763074/1‑150 (MQ=255)

                           aGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCttgtt                                                                                                                           <  2:1855663/150‑1 (MQ=255)

                            gTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCttgttt                                                                                                                          <  1:1585246/150‑1 (MQ=255)

                             tGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCttgtttg                                                                                                                         <  1:1066778/150‑1 (MQ=255)

                             tGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCttgtttg                                                                                                                         <  2:570529/150‑1 (MQ=255)

                              gCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCttgtttgt                                                                                                                        >  2:429061/1‑150 (MQ=255)

                               cATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTa                                                                                                                       <  1:1807465/150‑1 (MQ=255)

                                aTCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTat                                                                                                                      >  2:1570035/1‑149 (MQ=255)

                                aTCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTaa                                                                                                                      >  1:1628506/1‑150 (MQ=255)

                                 tCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAAt                                                                                                                     >  1:1430756/1‑150 (MQ=255)

                                    tttAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAAtttt                                                                                                                  >  2:1674325/1‑150 (MQ=255)

                                    tttAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAAATGCTTCAAGCACTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATAGGCATGTTGTATTTTGTCTTGTTTGTAAtttt                                                                                                                  >  1:1769241/1‑150 (MQ=255)

                                     ttAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTg                                                                                                                 <  2:1593076/150‑1 (MQ=255)

                                      tAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGt                                                                                                                <  2:1023479/150‑1 (MQ=255)

                                      tAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGt                                                                                                                >  2:622185/1‑150 (MQ=255)

                                      tAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGt                                                                                                                >  2:909587/1‑150 (MQ=255)

                                          aCCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCt                                                                                                            <  1:46492/150‑1 (MQ=255)

                                            cAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTg                                                                                                          >  1:287664/1‑150 (MQ=255)

                                                  ggCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACg                                                                                                    <  1:2214587/150‑1 (MQ=255)

                                                   gCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACgg                                                                                                   >  1:1754574/1‑150 (MQ=255)

                                                        aaCCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGgcttgc                                                                                              >  2:2041051/1‑150 (MQ=255)

                                                        aaCCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGgcttgc                                                                                              >  2:2032565/1‑150 (MQ=255)

                                                          ccGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGgcttgctt                                                                                            >  1:1872195/1‑150 (MQ=255)

                                                                cACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCcggcg                                                                                      >  1:1062763/1‑150 (MQ=255)

                                                                  ccAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCggcggc                                                                                    <  2:557097/150‑1 (MQ=255)

                                                                          aatgaatgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTAc                                                                            >  1:795843/1‑150 (MQ=255)

                                                                          aatgaatgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTAc                                                                            >  2:2267902/1‑150 (MQ=255)

                                                                          aatgaatgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTAc                                                                            <  2:106950/150‑1 (MQ=255)

                                                                             gaatgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTTGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCgg                                                                         <  2:454488/150‑1 (MQ=255)

                                                                              aatgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCggg                                                                        <  1:1076109/150‑1 (MQ=255)

                                                                               atgGAACTTTTGACCCCATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGACACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGc                                                                       <  1:1706245/150‑1 (MQ=255)

                                                                                  gAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGAc                                                                    >  2:1130005/1‑150 (MQ=255)

                                                                                   aaCTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACa                                                                   <  2:1779310/150‑1 (MQ=255)

                                                                                    aCTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAg                                                                  <  1:257173/150‑1 (MQ=255)

                                                                                     cTTTTGACCCAATTGCTGCACGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGt                                                                 >  2:2209246/1‑150 (MQ=255)

                                                                                     cTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGt                                                                 >  2:1449324/1‑150 (MQ=255)

                                                                                     cTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGt                                                                 >  1:1453656/1‑150 (MQ=255)

                                                                                      ttttGCCCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGtt                                                                >  1:2399294/1‑150 (MQ=255)

                                                                                      ttttGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGtt                                                                <  1:1674325/150‑1 (MQ=255)

                                                                                         tGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTAc                                                             <  1:909587/150‑1 (MQ=255)

                                                                                         tGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTAc                                                             <  2:763074/150‑1 (MQ=255)

                                                                                           aCCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTg                                                           <  2:1729137/150‑1 (MQ=255)

                                                                                            cccAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTgg                                                          >  2:740291/1‑150 (MQ=255)

                                                                                            cccAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTgg                                                          <  1:469095/150‑1 (MQ=255)

                                                                                            cccAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTgg                                                          <  1:1922787/150‑1 (MQ=255)

                                                                                              cAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTATCGGGCGACAGTTTACTGGTa                                                        <  1:680419/150‑1 (MQ=255)

                                                                                              cAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTa                                                        >  1:576551/1‑150 (MQ=255)

                                                                                                  tgctgcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTgg                                                    <  2:2419194/150‑1 (MQ=255)

                                                                                                     tgcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCg                                                 >  1:999798/1‑150 (MQ=255)

                                                                                                      gcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCgg                                                <  1:1217804/150‑1 (MQ=255)

                                                                                                       cAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGc                                               <  1:1796492/150‑1 (MQ=255)

                                                                                                        aaGCCCTGTGGGCGGAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCg                                              <  2:645983/150‑1 (MQ=255)

                                                                                                            ccTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGtt                                          <  2:1059424/150‑1 (MQ=255)

                                                                                                            ccTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGtt                                          <  2:1158222/150‑1 (MQ=255)

                                                                                                             cTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTg                                         <  2:1962897/150‑1 (MQ=255)

                                                                                                               gtgGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGAt                                       <  2:790474/150‑1 (MQ=255)

                                                                                                               gtgGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGAt                                       <  1:198381/150‑1 (MQ=255)

                                                                                                               gattggccAGGGTTTCAAAATCCTGCGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGa                                        >  1:285534/8‑150 (MQ=255)

                                                                                                                  ggCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGc                                    <  1:2306318/150‑1 (MQ=255)

                                                                                                                       aGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAg                               <  1:1356230/150‑1 (MQ=255)

                                                                                                                        ggATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAgg                              >  2:1987628/1‑150 (MQ=255)

                                                                                                                        ggATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTACTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCCGCGTGTTGATTGCGAAAgg                              >  2:218856/1‑150 (MQ=255)

                                                                                                                          tttttGAAAGCCTC‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTATGTAATTATGATCGTTGAAAAGGGCATCCTTCCGGCGGCCTTTTTACCGGGCGACAGATTACTCGTATTGGTCGGCGTGTTGAATGCGAAAGgcg                            <  1:2188001/149‑1 (MQ=255)

                                                                                                                          aatttGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGgcg                            >  2:1324849/3‑150 (MQ=255)

                                                                                                                           ttttGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGgcgc                           <  1:1978818/150‑1 (MQ=255)

                                                                                                                           ttttGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCAGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGgcgc                           >  2:2392264/1‑150 (MQ=255)

                                                                                                                             ttGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGa                         >  2:1162869/1‑150 (MQ=255)

                                                                                                                               gAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATg                       >  2:305373/1‑150 (MQ=255)

                                                                                                                               gAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATg                       <  1:1297251/150‑1 (MQ=255)

                                                                                                                               gAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGGGTTGATTGCGAAAGGCGCGATg                       >  2:231885/1‑150 (MQ=255)

                                                                                                                               gAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGATGTATTTTGTCTGGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTCGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATg                       >  2:2486838/1‑150 (MQ=255)

                                                                                                                                aaaCCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATgg                      <  1:723355/150‑1 (MQ=255)

                                                                                                                                  aCCCTG‑GCGAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGc                    <  2:2243504/150‑1 (MQ=255)

                                                                                                                                  aCCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGc                    <  1:1056734/150‑1 (MQ=255)

                                                                                                                                  aCCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGc                    >  1:1025645/1‑150 (MQ=255)

                                                                                                                                    aacatgCCAAT‑CATCGATGGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCt                   >  2:342684/6‑150 (MQ=255)

                                                                                                                                     cTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTAt                 >  2:729053/1‑150 (MQ=255)

                                                                                                                                     cTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTAt                 >  2:1065242/1‑150 (MQ=255)

                                                                                                                                       g‑gCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATcc               >  1:236486/1‑150 (MQ=255)

                                                                                                                                              tCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAAc         <  2:2272488/150‑1 (MQ=255)

ATACCAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTG‑GCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTG  >  NC_000913/3247694‑3247991

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