BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,481,285	A→G	A67A (GCA→GCG)	hrpA →	helicase, ATP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,481,285	0	A	G	99.4%	552.7 / ‑5.4	173	A67A (GCA→GCG)	hrpA	helicase, ATP‑dependent
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (64/108); total (64/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 176 total aligned reads displayed.

TAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAT  >  NC_000913/1481195‑1481374

tAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGcgcg                                                                                         <  2:2641568/93‑1 (MQ=255)

  aaaaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGAc                                                                                       <  1:1711763/93‑1 (MQ=255)

        cTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGaaa                                                                                 <  2:1557963/93‑1 (MQ=255)

              cacaACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACtt                                                                           <  2:4523527/93‑1 (MQ=255)

              caaaacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACtt                                                                           <  2:1686696/90‑1 (MQ=255)

               acaacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGAGCGGAAATTACTTa                                                                          <  2:4072163/93‑1 (MQ=255)

               acaacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGGAATTACTTa                                                                          >  2:4236150/1‑93 (MQ=255)

               acaacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTa                                                                          <  1:4313288/93‑1 (MQ=255)

               acaacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTa                                                                          <  1:1518214/93‑1 (MQ=255)

               acaacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTa                                                                          >  2:3460953/1‑93 (MQ=255)

               acaacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTa                                                                          <  2:3802775/93‑1 (MQ=255)

                caacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTAt                                                                         <  2:1756936/93‑1 (MQ=255)

                caacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTAt                                                                         <  2:2971301/93‑1 (MQ=255)

                 aacaGGCCATTTTCCAGGAGATGGCGAACGAGATTGACCAGGCGGCAGGGAAAGTCCTGCAGCGTGAAGCGGCGCGACCGGAAATTACTTATc                                                                        <  1:4214485/93‑1 (MQ=255)

                 aacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATc                                                                        >  1:2872408/1‑93 (MQ=255)

                 aacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATc                                                                        >  2:674638/1‑93 (MQ=255)

                 aacaGACCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATc                                                                        >  2:1392071/1‑93 (MQ=255)

                  acaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATcc                                                                       <  2:1162337/93‑1 (MQ=255)

                    aGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTg                                                                     <  2:3991505/93‑1 (MQ=255)

                     ggCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTgc                                                                    >  2:973050/1‑92 (MQ=255)

                     ggCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGa                                                                    >  2:2302921/1‑93 (MQ=255)

                        cATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACaa                                                                 >  2:4327271/1‑93 (MQ=255)

                         aTTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAAt                                                                <  2:796622/93‑1 (MQ=255)

                         aTTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAAt                                                                >  1:1210525/1‑93 (MQ=255)

                          ttttCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAAtt                                                               >  2:1243130/1‑93 (MQ=255)

                          ttttCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAAtt                                                               <  1:3732683/93‑1 (MQ=255)

                           tttCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAAttt                                                              <  2:3580401/93‑1 (MQ=255)

                           tttCCAGGAGATGGCGAAAGAGATGGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAAttt                                                              <  1:363099/93‑1 (MQ=255)

                            ttCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTa                                                             <  1:715744/93‑1 (MQ=255)

                             tCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTAc                                                            <  1:3974296/93‑1 (MQ=255)

                             tCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCC‑GACAATTTAcc                                                           >  2:1045515/1‑93 (MQ=255)

                              ccAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTAcc                                                           <  1:1242699/93‑1 (MQ=255)

                                aGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCgg                                                         >  2:1278799/1‑93 (MQ=255)

                                aGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCgg                                                         >  1:2155428/1‑93 (MQ=255)

                                 ggAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGt                                                        <  2:1086164/93‑1 (MQ=255)

                                 ggAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGt                                                        >  1:4520206/1‑93 (MQ=255)

                                  gagaTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGtt                                                       >  1:4726506/1‑93 (MQ=255)

                                  gagaTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGtt                                                       <  2:4470330/93‑1 (MQ=255)

                                  gagaTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGtt                                                       <  1:96709/93‑1 (MQ=255)

                                  gagaTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGtt                                                       <  2:715108/93‑1 (MQ=255)

                                   agaTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTa                                                      >  2:559028/1‑93 (MQ=255)

                                    gaTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAg                                                     >  1:3125207/1‑93 (MQ=255)

                                    gaTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAg                                                     >  1:531042/1‑93 (MQ=255)

                                    gaTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAg                                                     <  1:3722843/93‑1 (MQ=255)

                                     aTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCCGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGt                                                    <  2:3237231/93‑1 (MQ=255)

                                      tGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTc                                                   <  1:1350393/93‑1 (MQ=255)

                                       ggCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCa                                                  <  2:4733608/93‑1 (MQ=255)

                                       ggCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCa                                                  <  1:614688/93‑1 (MQ=255)

                                       ggCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGCGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCa                                                  >  1:2572642/1‑93 (MQ=255)

                                          gAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGaa                                               >  2:2140433/1‑93 (MQ=255)

                                          gAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGaa                                               >  1:225607/1‑93 (MQ=255)

                                            aaGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTTAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGaaca                                             >  1:2281943/1‑91 (MQ=255)

                                                gaTTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACa                                         <  1:4368935/93‑1 (MQ=255)

                                                gaTTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACa                                         >  2:4675214/1‑93 (MQ=255)

                                                 aTTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAg                                        >  2:3981233/1‑93 (MQ=255)

                                                  ttGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAgg                                       <  2:1708345/93‑1 (MQ=255)

                                                    gACCAGGCGGCAGGGAAAGTCCTGCTGCGTGACGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGAc                                     <  2:1343771/93‑1 (MQ=255)

                                                    gACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGAc                                     <  2:692854/93‑1 (MQ=255)

                                                    gACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGAc                                     <  1:1393157/93‑1 (MQ=255)

                                                       cAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACAtt                                  <  1:2685487/93‑1 (MQ=255)

                                                        aGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATtc                                 >  2:2854900/1‑93 (MQ=255)

                                                        aGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATtc                                 <  1:4014471/93‑1 (MQ=255)

                                                          gcggcAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATtctc                               >  1:2941613/1‑93 (MQ=255)

                                                          gcggcAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATtctc                               <  1:1106962/93‑1 (MQ=255)

                                                            ggcAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGa                             <  1:812653/93‑1 (MQ=255)

                                                             gcAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGaa                            <  1:468469/93‑1 (MQ=255)

                                                              cAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAg                           <  1:1465917/93‑1 (MQ=255)

                                                              cAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAg                           >  2:3643132/1‑93 (MQ=255)

                                                               aGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCATAAAAAACAGGACATTCTCGAAGc                          <  2:1187332/93‑1 (MQ=255)

                                                               aGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGc                          >  1:3167301/1‑93 (MQ=255)

                                                               aGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGc                          >  1:3553712/1‑93 (MQ=255)

                                                                gggAAAGTCCTGCTGCGTGAAGCGGTGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCg                         >  2:3880860/1‑93 (MQ=255)

                                                                gggAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCg                         >  2:4323643/1‑93 (MQ=255)

                                                                gggAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCg                         <  2:3394175/93‑1 (MQ=255)

                                                                gggAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCg                         >  2:3202069/1‑93 (MQ=255)

                                                                 ggAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGa                        >  1:661500/1‑93 (MQ=255)

                                                                  gAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGAt                       <  1:1169098/93‑1 (MQ=255)

                                                                  gAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGAt                       <  2:1897630/93‑1 (MQ=255)

                                                                   aaaGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGAtt                      <  1:2829172/93‑1 (MQ=255)

                                                                   aaaGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGAtt                      <  1:3581932/93‑1 (MQ=255)

                                                                     aGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCg                    <  2:1491397/93‑1 (MQ=255)

                                                                     aGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCg                    <  1:3355499/93‑1 (MQ=255)

                                                                      gTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGt                   >  2:2787426/1‑93 (MQ=255)

                                                                      gTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGt                   >  2:1011848/1‑93 (MQ=255)

                                                                      gTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGt                   <  1:1870556/93‑1 (MQ=255)

                                                                      gTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGt                   <  2:623775/93‑1 (MQ=255)

                                                                       tCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTg                  <  1:3854199/93‑1 (MQ=255)

                                                                       tCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTg                  <  1:3097621/93‑1 (MQ=255)

                                                                       tCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTg                  >  1:932010/1‑93 (MQ=255)

                                                                        ccTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGa                 <  2:1566688/93‑1 (MQ=255)

                                                                         ctgctgCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGAt                <  2:239762/93‑1 (MQ=255)

                                                                          tgctgcGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATc               >  1:3822456/1‑93 (MQ=255)

                                                                           gctgcGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCa              >  1:2795520/1‑93 (MQ=255)

                                                                               cGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAg          <  1:1720267/93‑1 (MQ=255)

                                                                               cGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAg          <  2:2001071/93‑1 (MQ=255)

                                                                                gTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAgg         <  1:3274201/93‑1 (MQ=255)

                                                                                 tGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAggt        <  2:2222205/93‑1 (MQ=255)

                                                                                 tGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAggt        <  2:3081458/93‑1 (MQ=255)

                                                                                  gAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAggtg       <  1:75979/93‑1 (MQ=255)

                                                                                       ggCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAt  <  2:1767340/93‑1 (MQ=255)

TAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAT  >  NC_000913/1481195‑1481374

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