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.7%	922.2 / ‑5.0	290	A67A (GCA→GCG)	hrpA	helicase, ATP‑dependent
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (134/155); total (135/155)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.58e-01

Only 100 of 300 total aligned reads displayed.

TAAAAATCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCACGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAT  >  NC_000913/1481195‑1481374

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

 aaaaaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGa                                                                                        <  2:1869851/93‑1 (MQ=255)

 aaaaaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGAGGCGCGa                                                                                        >  1:4028574/1‑93 (MQ=255)

  aaaaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGAc                                                                                       >  1:9691287/1‑93 (MQ=255)

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

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

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

   aaaTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGAcc                                                                                      >  1:7504276/1‑93 (MQ=255)

     aTCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCgg                                                                                    >  2:7805397/1‑93 (MQ=255)

      tCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGa                                                                                   >  2:2210714/1‑93 (MQ=255)

      tCCTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGa                                                                                   >  2:1140035/1‑93 (MQ=255)

       ccTGATGCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGaa                                                                                  >  1:220305/1‑93 (MQ=255)

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

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

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

             gCACAACAGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACt                                                                            >  1:5936526/1‑93 (MQ=255)

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

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

                caacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGGCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTAt                                                                         >  2:124613/1‑93 (MQ=255)

                 aacaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATc                                                                        <  2:2941057/93‑1 (MQ=255)

                  acaGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATcc                                                                       >  1:5821305/1‑93 (MQ=255)

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

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

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

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

                   caGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCt                                                                      >  2:5438129/1‑93 (MQ=255)

                   caGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCt                                                                      >  1:2896375/1‑93 (MQ=255)

                    aGGCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTg                                                                     >  2:1468195/1‑93 (MQ=255)

                     ggCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGa                                                                    <  1:7164342/93‑1 (MQ=255)

                      gCCATTTTCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGAc                                                                   >  2:223299/1‑93 (MQ=255)

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

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

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

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

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

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

                           tttCCAGGAGATGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAAttt                                                              >  2:7616266/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

                                   agaTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTa                                                      <  1:4848420/93‑1 (MQ=255)

                                    gaTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACCATTTACCGGTTAg                                                     >  2:1014956/1‑93 (MQ=255)

                                     agggCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGt                                                    <  1:8544624/91‑1 (MQ=255)

                                     aTGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGt                                                    <  2:9946840/93‑1 (MQ=255)

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

                                      tGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGCCCGGCAATTACTTATCCTGCCAATTTACCGGTTAGTc                                                   >  2:1673559/1‑93 (MQ=255)

                                      tGGCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTc                                                   >  1:1772026/1‑93 (MQ=255)

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

                                        gCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAg                                                 >  1:2315932/1‑93 (MQ=255)

                                        gCGAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAg                                                 <  2:936348/93‑1 (MQ=255)

                                          gAAAGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCTGaa                                               >  1:3947884/1‑93 (MQ=255)

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

                                           aaaGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGaaa                                              >  2:6992155/1‑93 (MQ=255)

                                           aaaGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGaaa                                              >  1:473334/1‑93 (MQ=255)

                                            aaGAGATTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGaaaa                                             >  2:5832393/1‑93 (MQ=255)

                                              gagaTTGACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGaaaaca                                           >  2:5613205/1‑91 (MQ=255)

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

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

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

                                                    gACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGAc                                     >  2:636604/1‑93 (MQ=255)

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

                                                    gACCAGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGAc                                     >  2:7315428/1‑93 (MQ=255)

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

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

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

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

                                                        aGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATtc                                 <  2:2833353/93‑1 (MQ=255)

                                                        aGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATtc                                 <  1:551767/93‑1 (MQ=255)

                                                        aGGCGGCAGGGAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATtc                                 <  2:1436467/93‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                  gAAAGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCGGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGAt                       <  2:3817246/93‑1 (MQ=255)

                                                                     aGTCCTGCTGCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCg                    >  2:6576790/1‑93 (MQ=255)

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

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

                                                                      gTCCTGCTGGGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGt                   <  1:1382618/93‑1 (MQ=255)

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

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

                                                                        ccTGCTGCGTGAAGCGGCGCGACCGGAAATTACGTAGCCTGACAATTTACCGGTTAGTCAGAAAAGACAGGACATTCACGAAGCGATTCGTGa                 <  2:9833885/93‑1 (MQ=255)

                                                                         ctgctgCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGAt                >  2:8073581/1‑93 (MQ=255)

                                                                         ctgctgCGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGAt                >  1:1547415/1‑93 (MQ=255)

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

                                                                            ctgcGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCAc             <  1:4235329/93‑1 (MQ=255)

                                                                             tgcGTGAAGCGGCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCAcc            <  1:3740526/93‑1 (MQ=255)

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

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

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

                                                                                     gcggcgCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGgtggtg    >  1:183680/1‑93 (MQ=255)

                                                                                       ggCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTTGTGATCACCAGGTGGTGAt  >  2:7528538/1‑93 (MQ=255)

                                                                                       ggCGCGACCGGAAATTACTTATCCTGACAATTTACCGGTTAGTCAGAAAAAACAGGACATTCTCGAAGCGATTCGTGATCACCAGGTGGTGAt  >  1:7341055/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: `‑`