BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,247,836	Δ1 bp	5.6%	coding (64/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,836	0	C	.	5.6%	555.0 / 22.7	162	coding (64/663 nt)	yqjA	DedA family protein YqjA
Reads supporting (aligned to +/- strand): ref base C (78/74); new base . (2/7); total (80/82)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 163 total aligned reads displayed.

GCAGAACATACCAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAAC  >  NC_000913/3247687‑3247984

   gAACATACCAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCAt                                                                                                                                                    >  2:1297251/1‑150 (MQ=255)

   gAACATACCAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCAt                                                                                                                                                    >  2:1796492/1‑150 (MQ=255)

    aaCATACCAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATc                                                                                                                                                   <  1:857309/150‑1 (MQ=255)

        tACCAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCgat‑g                                                                                                                                              >  1:446212/1‑150 (MQ=255)

        tACCAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCgat‑g                                                                                                                                              <  1:1791887/150‑1 (MQ=255)

        tACCAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCgat‑g                                                                                                                                              >  2:1356230/1‑150 (MQ=255)

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

          ccAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCACTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCgat‑gat                                                                                                                                            >  2:2054226/1‑150 (MQ=255)

           cAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATc                                                                                                                                                   <  1:2227886/143‑1 (MQ=255)

           cAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATc                                                                                                                                                   >  2:2227886/1‑143 (MQ=255)

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

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

               aCAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCa                                                                                                                                       <  1:2308862/150‑1 (MQ=255)

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

                 aGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAAT‑CATCGATGGATTGGccag                                                                                                                                     <  1:1831812/150‑4 (MQ=255)

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

                        aCATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTAtt                                                                                                                              <  1:271019/150‑1 (MQ=255)

                        aCATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTAtt                                                                                                                              >  2:1413852/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

                                     gCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCttgtttgt                                                                                                                 >  2:1594799/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

                                                   cAGCGTGTCGTTAACCGATTCACCAGGAATATTGAATGGAAGTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGTCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTg                                                                                                   <  1:785700/150‑1 (MQ=255)

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

                                                       gTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGaaaa                                                                                               <  1:2015418/150‑1 (MQ=255)

                                                        tGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAAc                                                                                              >  1:524148/1‑150 (MQ=255)

                                                          gCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACgg                                                                                            >  1:1729137/1‑150 (MQ=255)

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

                                                              tAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGgcttg                                                                                        >  2:1500681/1‑150 (MQ=255)

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

                                                                aCCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGTGCGCAGGATTTTGAAACCCTGGCAAATACATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGgcttgct                                                                                      <  2:563769/150‑1 (MQ=255)

                                                                aCCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGgcttgct                                                                                      >  1:428146/1‑150 (MQ=255)

                                                                     ttCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCcgg                                                                                 <  2:2337892/150‑1 (MQ=255)

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

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

                                                                          cAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGcc                                                                            >  1:850779/1‑150 (MQ=255)

                                                                             gAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCttt                                                                         <  1:582738/150‑1 (MQ=255)

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

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

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

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

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

                                                                                       tgGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCg                                                               >  2:226664/1‑150 (MQ=255)

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

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

                                                                                           aCTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAg                                                           <  1:1954980/150‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                                                          gctgcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGt                                            <  2:2412078/150‑1 (MQ=255)

                                                                                                          gctgcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGt                                            >  1:1843304/1‑150 (MQ=255)

                                                                                                          gctgcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGt                                            >  2:2017729/1‑150 (MQ=255)

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

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

                                                                                                             gcAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCgg                                         >  1:683280/1‑150 (MQ=255)

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

                                                                                                               aaGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGAGTCAGTTTACTGGTATTGGTCGGCg                                       <  2:1463887/150‑1 (MQ=255)

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

                                                                                                                   ccTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGtt                                   <  1:1449324/150‑1 (MQ=255)

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

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

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

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

                                                                                                                              aGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAg                        <  2:2309394/150‑1 (MQ=255)

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

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

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

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

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

                                                                                                                                   tttGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGcgcg                   <  2:1004038/150‑1 (MQ=255)

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

                                                                                                                                     tGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGAt                 <  2:1754574/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                                                                                                             atgCCAAT‑CATCGATGGATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATc         >  2:1820807/3‑150 (MQ=255)

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

GCAGAACATACCAAAACAGCAATAACATTGCGGTAGTGCATCTTTTAAAACCAGCGTGGCGTTAACCGATTCACCAGGAATAATGAATGGAACTTTTGACCCAATTGCTGCAAGCCCTGTGGGCGCAGGATTTTGAAACCCTGGCCAATCCATCGAT‑GATTGGCATGTTGTATTTTGTCTTGTTTGTAATTTTGTTCCTTGAAAACGGCTTGCTTCCGGCGGCCTTTTTACCGGGCGACAGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAAC  >  NC_000913/3247687‑3247984

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