BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,902,495	A→C	5.5%	G200G (GGA→GGC) ‡	ygcE →	putative kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,902,495	0	A	C	5.5%	378.8 / 12.7	110	G200G (GGA→GGC) ‡	ygcE	putative kinase
Reads supporting (aligned to +/- strand): ref base A (51/53); new base C (4/2); total (55/55)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.81e-01

Only 100 of 111 total aligned reads displayed.

TATAAGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACT  >  NC_000913/2902349‑2902643

tataAGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGAc                                                                                                                                                     >  2:249147/1‑148 (MQ=255)

 ataAGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTg                                                                                                                                                   >  3:909143/1‑149 (MQ=255)

  taAGAGCTGGGCATGGAGCGAAGACGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATTTGCAAATGCCTGGCACCGTCCTCGGTCATATCACACCACAAGCCGCTCTGGCGACACATTTCCCGGCAGGACTGc                                                                                                                                                  >  3:599402/1‑149 (MQ=255)

   aaGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGcc                                                                                                                                                 >  4:791073/1‑149 (MQ=255)

   aaGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGc                                                                                                                                                  >  2:953724/1‑148 (MQ=255)

    agagCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGcc                                                                                                                                                 >  1:638408/1‑148 (MQ=255)

     gagCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCgg                                                                                                                                               <  3:480575/149‑1 (MQ=255)

       gCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgt                                                                                                                                              <  1:882686/148‑1 (MQ=255)

          gggCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgggaaat                                                                                                                                          <  4:792407/149‑6 (MQ=255)

          gggCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgttgtt                                                                                                                                          <  4:214545/149‑1 (MQ=255)

          gggCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgttgt                                                                                                                                           <  2:203275/148‑1 (MQ=255)

           ggCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgttgtt                                                                                                                                          <  2:412179/148‑1 (MQ=255)

            gCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGTCCTGCCgggaaatgt                                                                                                                                        <  4:475811/149‑8 (MQ=255)

                          tgctgcGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGAt                                                                                                                           >  1:702233/1‑148 (MQ=255)

                            ctgcGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATaa                                                                                                                         <  2:644090/148‑1 (MQ=255)

                              gcGGTAATGGATAAGTTTATTATCCCCCGTCATATGCTGTTTGATGAGCAAATGGCTGGCAGCTTCCTCGCACATATCACACCACATGCCGCACTGGCGGCACATTTCCCGGCAGGACTGCCGGTTGTTTTTACCACTAGTGAAAAAcc                                                                                                                      <  4:186614/149‑1 (MQ=255)

                              gcGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAAc                                                                                                                       >  2:471272/1‑148 (MQ=255)

                               cGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCg                                                                                                                     <  3:587876/149‑1 (MQ=255)

                                         tAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGct                                                                                                            <  2:855774/148‑1 (MQ=255)

                                             tttAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTggg                                                                                                       >  3:875582/1‑149 (MQ=255)

                                             tttAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTggg                                                                                                       <  3:791073/149‑1 (MQ=255)

                                             tttAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTggg                                                                                                       >  4:558307/1‑149 (MQ=255)

                                             tttAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTggg                                                                                                       >  4:185732/1‑149 (MQ=255)

                                              ttAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTg                                                                                                                             >  2:835523/1‑126 (MQ=255)

                                              ttAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTg                                                                                                                             <  1:835523/126‑1 (MQ=255)

                                              ttAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTgggg                                                                                                      <  4:728360/149‑1 (MQ=255)

                                                  tatCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCg                                                                                                  >  3:428966/1‑149 (MQ=255)

                                                    tCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgttgtt                                                                                                                                          <  1:124609/107‑1 (MQ=255)

                                                    tCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgttgtt                                                                                                                                          >  2:124609/1‑107 (MQ=255)

                                                     cccccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGa                                                                                                >  1:43418/1‑148 (MQ=255)

                                                     cccccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGa                                                                                                >  2:542090/1‑148 (MQ=255)

                                                     cccccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGa                                                                                                >  2:438378/1‑148 (MQ=255)

                                                        ccGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTa                                                                                             <  2:681465/148‑1 (MQ=255)

                                                         cGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTAc                                                                                            >  2:194151/1‑148 (MQ=255)

                                                         cGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACt                                                                                           <  4:701630/149‑1 (MQ=255)

                                                         cGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGAGAAACCGGTAGAAGCTCTGGGGGCCGGATTAc                                                                                            <  2:168446/148‑1 (MQ=255)

                                                            cATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTgg                                                                                         >  2:744188/1‑148 (MQ=255)

                                                             atatGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGga                                                                                        <  1:827048/148‑1 (MQ=255)

                                                              tatGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGgat                                                                                       >  1:47852/1‑148 (MQ=255)

                                                               atGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGgatg                                                                                      <  2:134989/148‑1 (MQ=255)

                                                                tGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGgatga                                                                                     >  2:628024/1‑148 (MQ=255)

                                                                 gCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGgatgat                                                                                    >  2:595808/1‑148 (MQ=255)

                                                                  cTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGatgatg                                                                                   >  2:759604/1‑148 (MQ=255)

                                                                  cTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGatgatg                                                                                   >  1:855479/1‑148 (MQ=255)

                                                                    gTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGaa                                                                                 <  1:817505/148‑1 (MQ=255)

                                                                     tttGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGaatc                                                                               >  4:819541/1‑147 (MQ=255)

                                                                     tttGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAAc                                                                               >  4:517034/1‑149 (MQ=255)

                                                                     tttGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTACCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAAc                                                                               >  3:665274/1‑149 (MQ=255)

                                                                     tttGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAACCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAAc                                                                               >  4:908870/1‑149 (MQ=255)

                                                                         aTGTGCAAATGCCTGGCACCGTCCTGGGACATATCACACCACAAGCCGCACTGTCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAAcggc                                                                            >  1:501617/1‑148 (MQ=255)

                                                                         aTGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAAcggcg                                                                           <  4:476575/149‑1 (MQ=255)

                                                                         aTGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAAcggc                                                                            <  1:542090/148‑1 (MQ=255)

                                                                                aaTGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAAttt                                                                    <  4:290298/149‑1 (MQ=255)

                                                                                aaTGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAAtt                                                                     >  1:220730/1‑148 (MQ=255)

                                                                                    ccTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTtctc                                                                 <  1:471272/148‑1 (MQ=255)

                                                                                     cTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTctct                                                                <  1:595808/148‑1 (MQ=255)

                                                                                     cTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTctct                                                                <  1:907918/148‑1 (MQ=255)

                                                                                      tGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCtt                                                               <  2:649798/148‑1 (MQ=255)

                                                                                       ggCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTg                                                              >  1:264615/1‑148 (MQ=255)

                                                                                         cACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGc                                                            >  2:496081/1‑148 (MQ=255)

                                                                                              tCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCTCGGCAGGACTGCCGGttgttt                                                                                                                                         >  4:587943/1‑66 (MQ=255)

                                                                                              tCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCTCGGCAGGACTGCCGGttgttt                                                                                                                                         <  3:587943/66‑1 (MQ=255)

                                                                                              tCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgttgtt                                                                                                                                          >  1:297334/1‑65 (MQ=255)

                                                                                              tCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGgttgtt                                                                                                                                          <  2:297334/65‑1 (MQ=255)

                                                                                                 tCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACAt                                                    >  1:409902/1‑148 (MQ=255)

                                                                                                   ggACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGc                                                 >  4:843095/1‑149 (MQ=255)

                                                                                                    gACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGc                                                 <  2:332112/148‑1 (MQ=255)

                                                                                                      cATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAATCCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCAtt                                              <  4:768832/149‑1 (MQ=255)

                                                                                                        tatCACACCACAAGCCGCACTGGCGACACAGTTCCCGGCAGGACTGCCGGTTGTTTGTACCTCCAGTGATAAACCGGTAGAAGCTCTGGGGGGCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtga                                            <  3:843095/149‑1 (MQ=255)

                                                                                                         atCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgat                                           <  3:355044/149‑1 (MQ=255)

                                                                                                         atCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgat                                           <  4:387019/149‑1 (MQ=255)

                                                                                                         atCACACCACAAGCCGCACTGGCGACACAATTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgat                                           <  3:542513/149‑1 (MQ=255)

                                                                                                          tcaatccacAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgat                                           >  2:65709/6‑148 (MQ=255)

                                                                                                           cacaCCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGga                                                                                        >  1:732945/1‑102 (MQ=255)

                                                                                                           cacaCCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGga                                                                                        <  2:732945/102‑1 (MQ=255)

                                                                                                            acacCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatgat                                        <  3:908870/149‑1 (MQ=255)

                                                                                                            acacCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatgat                                        >  3:587461/1‑149 (MQ=255)

                                                                                                            acacCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatga                                         <  1:700292/148‑1 (MQ=255)

                                                                                                             caccacAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATtgatgat                                        >  1:927631/1‑148 (MQ=255)

                                                                                                               ccacAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGaa                                     <  4:52434/149‑1 (MQ=255)

                                                                                                                cacaAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAAc                                    >  4:219659/1‑149 (MQ=255)

                                                                                                                cacaAGCCGCACTGGCGACACATATCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAAc                                    >  4:61780/1‑149 (MQ=255)

                                                                                                                 acaAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAAc                                    >  1:465999/1‑148 (MQ=255)

                                                                                                                      ccGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCaa                               <  1:5296/148‑1 (MQ=255)

                                                                                                                                  cacaTTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGaa                   <  1:744188/148‑1 (MQ=255)

                                                                                                                                  cacaTTTCCCGGCAGGACTGCCGATTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGaaa                  <  4:665274/149‑1 (MQ=255)

                                                                                                                                    aaacaaccGGCAGTCCTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGa                >  4:613584/7‑149 (MQ=255)

                                                                                                                                     aTTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGa                <  2:820631/148‑1 (MQ=255)

                                                                                                                                      tttCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGAt               >  2:695315/1‑148 (MQ=255)

                                                                                                                                       ttCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGGACGGCAAAGCACTGCCGAAAGATc              >  1:662438/1‑148 (MQ=255)

                                                                                                                                        tCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATcc             <  2:702233/148‑1 (MQ=255)

                                                                                                                                        aaccGGCAGTCCTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATcc             >  2:147619/3‑148 (MQ=255)

                                                                                                                                        aaccGGCAGTCCTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCa            >  3:398738/3‑149 (MQ=255)

                                                                                                                                         cccGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCa            <  2:488697/148‑1 (MQ=255)

                                                                                                                                         cccGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCa            <  1:438378/148‑1 (MQ=255)

                                                                                                                                         accGGCAGTCCTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCa            >  1:659958/2‑148 (MQ=255)

                                                                                                                                          ccGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGATCTGGGGGCCGGATTAGTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAAGGGCAAAGCACTGCCGAAAGATCCAGt          >  4:18740/1‑149 (MQ=255)

                                                                                                                                             gCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAg        <  2:43418/148‑1 (MQ=255)

                                                                                                                                                ggACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGCGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTa    >  3:114171/1‑149 (MQ=255)

                                                                                                                                                  aCTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACt  >  4:233182/1‑149 (MQ=255)

TATAAGAGCTGGGCATGGAGCGAAGATGCTGCGGTAATGGATAAGTTTAATATCCCCCGTCATATGCTGTTTGATGTGCAAATGCCTGGCACCGTCCTCGGACATATCACACCACAAGCCGCACTGGCGACACATTTCCCGGCAGGACTGCCGGTTGTTTGTACCACCAGTGATAAACCGGTAGAAGCTCTGGGGGCCGGATTACTGGATGATGAAACGGCGGTAATTTCTCTTGGCACTTACATCGCATTGATGATGAACGGCAAAGCACTGCCGAAAGATCCAGTAGCGTACT  >  NC_000913/2902349‑2902643

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