BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	955,785	+A	coding (71/1761 nt)	ycaO ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	955,782	1	.	A	98.4%	447.3 / ‑0.5	125	coding (74/1761 nt)	ycaO	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (2/0); new base A (67/56); total (69/56)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

Only 100 of 136 total aligned reads displayed.

AGTGCGCACTCTTTGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTG‑AAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAATCACTTTA  >  NC_000913/955692‑955868

aGTGCGCACTCTTTGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTG‑aa                                                                                      <  1:3776204/93‑1 (MQ=255)

aGTGCGCACTCTTTGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTG‑aa                                                                                      >  1:4593986/1‑93 (MQ=255)

aGTGCGCACTCTTTGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTG‑aa                                                                                      >  2:444296/1‑93 (MQ=255)

       aCTCTTTGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGtttt                                                                                <  2:1908044/93‑1 (MQ=255)

        ctctTTGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGttttt                                                                               >  2:3776633/1‑93 (MQ=255)

         tctTTGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTg                                                                              <  1:3906492/93‑1 (MQ=255)

            ttGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTg                                                                           >  2:430562/1‑93 (MQ=255)

              gTCGCGAATATGTACAGACCAGCCGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGa                                                                         <  1:3672675/93‑1 (MQ=255)

              gTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGa                                                                         <  1:4676330/93‑1 (MQ=255)

              gTCGCGAATAGGTACAGACCAGACGTTAGGCTCGGGATTCAGCCAGGAGGCCTCTTGAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGa                                                                         <  1:598025/93‑1 (MQ=255)

               tCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGaa                                                                        <  2:4420402/93‑1 (MQ=255)

                cgcgAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAag                                                                       >  2:1518743/1‑93 (MQ=255)

                  cgAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAagcg                                                                     >  1:1636419/1‑93 (MQ=255)

                   gAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAagcga                                                                    <  1:2957833/93‑1 (MQ=255)

                    aaTATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAagcgag                                                                   >  2:4552460/1‑93 (MQ=255)

                    aaTATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAagcgag                                                                   <  1:2863883/93‑1 (MQ=255)

                      tatGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAagtgagcg                                                                 >  1:2786224/1‑93 (MQ=255)

                      tatGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAagcgcgcg                                                                 >  1:4178268/1‑93 (MQ=255)

                      tatGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAagcgagcg                                                                 >  2:4504342/1‑93 (MQ=255)

                       atGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAgcgagcga                                                                <  1:3221852/93‑1 (MQ=255)

                        tGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGAt                                                               <  1:4470864/93‑1 (MQ=255)

                        tGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGAt                                                               >  1:1728619/1‑93 (MQ=255)

                        tGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGAt                                                               >  1:3323997/1‑93 (MQ=255)

                        tGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGAt                                                               >  1:4523168/1‑93 (MQ=255)

                         gTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATg                                                              <  2:2024898/93‑1 (MQ=255)

                           aCAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGa                                                            <  1:3649357/93‑1 (MQ=255)

                           aCAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGa                                                            <  2:596178/93‑1 (MQ=255)

                            cagaccagacGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGaa                                                           >  1:1824917/1‑93 (MQ=255)

                            cagaccagacGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGaa                                                           >  1:611777/1‑93 (MQ=255)

                              gaccagacGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATc                                                         >  1:4468428/1‑93 (MQ=255)

                               accagacGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCt                                                        >  1:2609558/1‑93 (MQ=255)

                               accagacGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCt                                                        >  2:2642600/1‑93 (MQ=255)

                                ccagacGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCtt                                                       >  2:677779/1‑93 (MQ=255)

                                ccagacGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCtt                                                       >  2:2777007/1‑93 (MQ=255)

                                 cagacGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGGGAGCGATGGAATCTTc                                                      >  2:4122610/1‑93 (MQ=255)

                                 cagacGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTc                                                      <  2:1242109/93‑1 (MQ=255)

                                  agacGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCGATCTGAAAGCCGAGGGCTGTAGAGTTTTTGCTGGGAGCGAGCGGTGGAATCTTcg                                                     >  2:495856/1‑92 (MQ=255)

                                  agacGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTcc                                                     >  1:2132432/1‑93 (MQ=255)

                                   gacGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCa                                                    <  1:1048668/93‑1 (MQ=255)

                                   gacGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCa                                                    >  2:2283457/1‑93 (MQ=255)

                                   gacGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCa                                                    >  2:3153104/1‑93 (MQ=255)

                                    acGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCag                                                   >  2:2923259/1‑93 (MQ=255)

                                       ttAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGGTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGc                                                >  2:3519784/1‑93 (MQ=255)

                                         aGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGAAATCTTCCAGAGCgg                                              <  2:2531451/93‑1 (MQ=255)

                                               gggATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCtt                                        >  1:665277/1‑93 (MQ=255)

                                               gggATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCtt                                        <  2:2926538/93‑1 (MQ=255)

                                               gggATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCtt                                        <  1:44215/93‑1 (MQ=255)

                                                ggATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGCAGCGAGCGATGGAATCTTCCAGAGCGGCATCttt                                       >  1:1167823/1‑93 (MQ=255)

                                                 gATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTG‑AAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGc                                     >  1:1485947/1‑93 (MQ=255)

                                                  aTTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGc                                     >  1:997012/1‑93 (MQ=255)

                                                   ttCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGcc                                    >  1:568115/1‑93 (MQ=255)

                                                     cAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCgg                                  <  2:1953676/93‑1 (MQ=255)

                                                        ccAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGa                               >  1:1092593/1‑93 (MQ=255)

                                                        ccAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGa                               <  2:4212040/93‑1 (MQ=255)

                                                           ggaggCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATaa                            <  1:795733/93‑1 (MQ=255)

                                                           ggaggCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATaa                            >  1:1855897/1‑93 (MQ=255)

                                                           ggaggCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATaa                            >  2:1500375/1‑93 (MQ=255)

                                                           ggaggCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATaa                            <  2:2809126/93‑1 (MQ=255)

                                                            gaggCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATaaa                           <  1:1995068/93‑1 (MQ=255)

                                                            gaggCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATaaa                           <  1:3345144/93‑1 (MQ=255)

                                                             aggCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAAt                          <  2:3627452/93‑1 (MQ=255)

                                                             aggCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAAt                          <  1:4139376/93‑1 (MQ=255)

                                                             aggCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAAt                          >  1:578405/1‑93 (MQ=255)

                                                              ggCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATg                         <  2:815508/93‑1 (MQ=255)

                                                              ggCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATg                         >  1:2920065/1‑93 (MQ=255)

                                                               gCCTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGt                        >  1:3916704/1‑93 (MQ=255)

                                                                ccTCTTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGtt                       <  2:2759528/93‑1 (MQ=255)

                                                                 ctctTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGttt                      <  1:2970662/93‑1 (MQ=255)

                                                                  tctTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTg                     >  1:1504693/1‑93 (MQ=255)

                                                                  tctTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTg                     >  2:3884708/1‑93 (MQ=255)

                                                                  tctTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTg                     <  2:1288882/93‑1 (MQ=255)

                                                                  tctTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTg                     >  2:1680392/1‑93 (MQ=255)

                                                                  tctTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTg                     >  1:2086893/1‑93 (MQ=255)

                                                                   ctTCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGc                    >  1:4620849/1‑93 (MQ=255)

                                                                    ttCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCg                   <  2:2096460/93‑1 (MQ=255)

                                                                     tCAATCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGt                  >  1:4394934/1‑93 (MQ=255)

                                                                       aaTCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCa                <  2:2146225/93‑1 (MQ=255)

                                                                       aaTCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGGTTGCGTCa                <  2:4557609/93‑1 (MQ=255)

                                                                        aTCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCAt               <  1:741645/93‑1 (MQ=255)

                                                                        aTCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCAt               >  1:244936/1‑93 (MQ=255)

                                                                        aTCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCAt               <  1:2463533/93‑1 (MQ=255)

                                                                        aTCTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCAt               >  2:1558130/1‑93 (MQ=255)

                                                                          cTGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGa             <  2:1844589/93‑1 (MQ=255)

                                                                          cTGAAAGCCGAGGTCTGAAAAGTGTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGAGATAAATGTTTGCGTCATGa             <  1:2611952/93‑1 (MQ=255)

                                                                           tGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGaa            <  1:2254407/93‑1 (MQ=255)

                                                                           tGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGaa            >  1:3578672/1‑93 (MQ=255)

                                                                           tGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGaa            >  2:1760318/1‑93 (MQ=255)

                                                                           tGAAAGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGaa            >  1:4412724/1‑93 (MQ=255)

                                                                             aaaGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGaaaa          <  1:3626431/93‑1 (MQ=255)

                                                                              aaGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAAt         <  2:4011477/93‑1 (MQ=255)

                                                                              aaGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAAt         >  1:2046842/1‑93 (MQ=255)

                                                                               aGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAATc        >  2:127869/1‑93 (MQ=255)

                                                                               aGCCGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAATc        >  2:2794266/1‑93 (MQ=255)

                                                                                 ccGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAATCAc      >  2:109539/1‑93 (MQ=255)

                                                                                  cGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGCAAATCACt     >  2:2412757/1‑93 (MQ=255)

                                                                                  cGAGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAATCACt     <  2:3954233/93‑1 (MQ=255)

                                                                                    aGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAATCACttt   >  2:795348/1‑93 (MQ=255)

                                                                                    aGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAATCACttt   <  1:727068/93‑1 (MQ=255)

                                                                                    aGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAATCACttt   <  1:3122007/93‑1 (MQ=255)

                                                                                    aGGTCTGAAAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAATCACttt   <  1:574719/93‑1 (MQ=255)

AGTGCGCACTCTTTGTCGCGAATATGTACAGACCAGACGTTAGGCACGGGATTCAGCCAGGAGGCCTCTTCAATCTGAAAGCCGAGGTCTG‑AAAGTTTTTGCTGGAAGCGAGCGATGGAATCTTCCAGAGCGGCATCTTTGCCGGGGATAAATGTTTGCGTCATGAAAATCACTTTA  >  NC_000913/955692‑955868

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