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.5%	628.8 / ‑4.9	204	A67A (GCA→GCG)	hrpA	helicase, ATP‑dependent
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (80/123); total (81/123)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.44e-01

Only 100 of 210 total aligned reads displayed.

TAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAT  >  NC_000913/1481195‑1481374

tAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGcgcg                                                                                         >  2:4076792/1‑93 (MQ=255)

 aaaaaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGa                                                                                        >  1:366662/1‑93 (MQ=255)

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

   aaaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGAcc                                                                                      <  2:2709817/93‑1 (MQ=255)

       ccTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGaa                                                                                  <  2:1010929/93‑1 (MQ=255)

        cTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGaaa                                                                                 >  1:5556281/1‑93 (MQ=255)

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

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

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

         tGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAAt                                                                                >  1:4020904/1‑93 (MQ=255)

           aTGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTa                                                                              <  1:5110457/93‑1 (MQ=255)

           aTGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTa                                                                              <  1:1518824/93‑1 (MQ=255)

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

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

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

                   caGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCt                                                                      <  2:1324886/93‑1 (MQ=255)

                   caGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCt                                                                      <  1:6837927/93‑1 (MQ=255)

                   caGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCt                                                                      <  2:7226457/93‑1 (MQ=255)

                      gCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGAc                                                                   <  1:5013519/93‑1 (MQ=255)

                        cATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGTAAATTACTTATCCTGACaa                                                                 <  1:3365519/93‑1 (MQ=255)

                        cATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACaa                                                                 <  1:5225443/93‑1 (MQ=255)

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

                         cttttCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAAt                                                                <  1:1342821/92‑1 (MQ=255)

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

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

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

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

                               cAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCg                                                          <  1:6600572/93‑1 (MQ=255)

                               cAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCg                                                          <  1:2597435/93‑1 (MQ=255)

                               cAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCg                                                          >  2:3472033/1‑93 (MQ=255)

                                aGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCgg                                                         <  1:6171824/93‑1 (MQ=255)

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

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

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

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

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

                                     aTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGt                                                    >  2:5903965/1‑93 (MQ=255)

                                       ggCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCa                                                  >  1:2835078/1‑93 (MQ=255)

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

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

                                         cGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGa                                                <  2:5685654/93‑1 (MQ=255)

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

                                            aaGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGaaaa                                             <  1:4343520/93‑1 (MQ=255)

                                             agagATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGaaaaa                                            <  1:6623146/93‑1 (MQ=255)

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

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

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

                                                 aTTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAg                                        <  1:3481594/93‑1 (MQ=255)

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

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

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

                                                   tGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGa                                      <  1:7007505/93‑1 (MQ=255)

                                                   tGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGa                                      <  1:3202491/93‑1 (MQ=255)

                                                     aCCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACa                                    <  2:766733/93‑1 (MQ=255)

                                                      ccAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACAt                                   <  1:782554/93‑1 (MQ=255)

                                                      ccAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACAt                                   <  2:1370588/93‑1 (MQ=255)

                                                       cAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACAtt                                  >  1:6432898/1‑93 (MQ=255)

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

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

                                                         ggcggcAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATtct                                >  2:6894929/1‑93 (MQ=255)

                                                           cggcAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCg                              <  2:2153218/93‑1 (MQ=255)

                                                            ggcAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGCCAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGa                             >  2:6079164/1‑93 (MQ=255)

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

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

                                                             gcAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGaa                            >  2:4623167/1‑93 (MQ=255)

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

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

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

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

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

                                                                  gAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGAt                       >  2:4335902/1‑93 (MQ=255)

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

                                                                    aaGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTc                     <  2:331716/93‑1 (MQ=255)

                                                                    aaGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTc                     <  1:2435468/93‑1 (MQ=255)

                                                                    aaGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTc                     >  1:5733225/1‑93 (MQ=255)

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

                                                                     aGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTAC‑TATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGt                   <  2:6103992/93‑1 (MQ=255)

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

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

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

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

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

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

                                                                        ccTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGa                 >  1:1598192/1‑93 (MQ=255)

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

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

                                                                            ctgcGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTCAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCAc             <  2:7058679/93‑1 (MQ=255)

                                                                             tgcGTGAAGGGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCAcc            <  1:2944983/93‑1 (MQ=255)

                                                                             tgcGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCAcc            >  2:6839665/1‑93 (MQ=255)

                                                                              gcGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCa           <  1:5419939/93‑1 (MQ=255)

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

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

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

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

                                                                                gTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAgg         >  2:7735025/1‑93 (MQ=255)

                                                                                gTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGAGATTCGTGGTCACCAgg         >  2:5846021/1‑93 (MQ=255)

                                                                                 tGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAggt        >  2:55999/1‑93 (MQ=255)

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

                                                                                    aGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAggtggt     <  1:6497079/93‑1 (MQ=255)

                                                                                     gcggcgCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGgtggtg    >  2:5945015/1‑93 (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` `< 30 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`