BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	656,969	A→G	W65R (TGG→CGG)	crcB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	656,969	0	A	G	99.6%	1564.1 / ‑4.1	481	W65R (TGG→CGG)	crcB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (250/229); total (251/230)
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.05e-01

Only 100 of 508 total aligned reads displayed.

GGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCA  >  NC_000913/656879‑657058

ggCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCcgca                                                                                         <  1:8785276/93‑1 (MQ=255)

 gCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCgcac                                                                                        >  1:9205010/1‑93 (MQ=255)

 gCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCgcac                                                                                        >  1:6094769/1‑93 (MQ=255)

  cccTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACt                                                                                       <  2:9754581/93‑1 (MQ=255)

  cccTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACt                                                                                       >  1:3670151/1‑93 (MQ=255)

    ctctTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCCCAAAATCCGGGGGTGATTAATACTTTCCGCACTgg                                                                                     >  2:1508291/1‑93 (MQ=255)

    ctctTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTgg                                                                                     >  1:2326468/1‑93 (MQ=255)

    ctctTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTgg                                                                                     <  1:5300820/93‑1 (MQ=255)

     tctTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGa                                                                                    <  1:9404166/93‑1 (MQ=255)

     tctTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGa                                                                                    <  1:7570116/93‑1 (MQ=255)

     tctTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGa                                                                                    <  1:5122105/93‑1 (MQ=255)

       ttGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATc                                                                                  <  2:7859037/93‑1 (MQ=255)

        tGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCa                                                                                 >  2:724069/1‑93 (MQ=255)

        tGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCa                                                                                 >  2:1898690/1‑93 (MQ=255)

          tAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAAt                                                                               >  1:9731250/1‑93 (MQ=255)

            aCAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGt                                                                             >  1:6652872/1‑93 (MQ=255)

            aCAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGt                                                                             >  2:3623840/1‑93 (MQ=255)

            aCAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGt                                                                             <  2:344217/93‑1 (MQ=255)

              aaaaaCCCCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTc                                                                           <  2:8317770/93‑1 (MQ=255)

               aaaaCACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCg                                                                          <  2:5540253/93‑1 (MQ=255)

                aaaCCCCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGt                                                                         <  1:9549221/93‑1 (MQ=255)

                aaaCACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGt                                                                         >  1:5404393/1‑93 (MQ=255)

                  acacCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCa                                                                       >  1:1530251/1‑93 (MQ=255)

                    cccacTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATc                                                                     <  2:7431424/92‑1 (MQ=255)

                    cccacTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATc                                                                     <  2:7970106/92‑1 (MQ=255)

                    accacTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATc                                                                     >  1:188103/1‑93 (MQ=255)

                    accacTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATc                                                                     >  2:4721627/1‑93 (MQ=255)

                     ccacTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATcc                                                                    <  2:5242979/93‑1 (MQ=255)

                     ccacTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACAGGATCAATGTTCGTCATct                                                                    >  1:3198637/1‑92 (MQ=255)

                       ccTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCctg                                                                  <  2:176075/92‑1 (MQ=255)

                       ccTTCTGCCAAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCctg                                                                  <  1:3381992/92‑1 (MQ=255)

                       acTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCctg                                                                  >  2:311589/1‑93 (MQ=255)

                       acTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCctg                                                                  >  1:992772/1‑93 (MQ=255)

                        cTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCctgc                                                                 >  1:8354411/1‑93 (MQ=255)

                         ttCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCctgct                                                                <  2:6821146/93‑1 (MQ=255)

                           cTGCCGAAAATGTGGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGa                                                              <  2:4122400/93‑1 (MQ=255)

                            tGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGaa                                                             >  1:9289154/1‑93 (MQ=255)

                              ccGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAAcc                                                           <  2:4864696/93‑1 (MQ=255)

                               cGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCa                                                          <  1:7275803/93‑1 (MQ=255)

                                gAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCAt                                                         <  1:5724225/93‑1 (MQ=255)

                                 aaaaTGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATg                                                        <  1:8267955/93‑1 (MQ=255)

                                 aaaaTGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATg                                                        >  1:8681077/1‑93 (MQ=255)

                                  aaaTGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGc                                                       <  2:8782231/93‑1 (MQ=255)

                                    aTGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGa                                                     <  1:9018649/93‑1 (MQ=255)

                                      gTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAt                                                   >  1:9990915/1‑93 (MQ=255)

                                         gagaAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtcct                                                >  1:5882751/1‑93 (MQ=255)

                                         gagaAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtcct                                                >  2:8468160/1‑93 (MQ=255)

                                             aGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtcctattc                                            >  2:2889812/1‑93 (MQ=255)

                                             aGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtcctattc                                            >  1:3099533/1‑93 (MQ=255)

                                              ggTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtcctattcc                                           <  1:978996/93‑1 (MQ=255)

                                               gttgttAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtcctattcct                                          >  1:1765358/1‑93 (MQ=255)

                                               gttgttAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtcctattcct                                          >  1:4190650/1‑93 (MQ=255)

                                                ttgttAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtcctattccta                                         <  1:3975594/93‑1 (MQ=255)

                                                 tgttAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtcctattcctat                                        >  2:5035145/1‑93 (MQ=255)

                                                 tgttAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAAtcctattcctat                                        >  1:7118277/1‑93 (MQ=255)

                                                  gttAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTatg                                       >  1:1349283/1‑93 (MQ=255)

                                                   ttAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTatga                                      >  1:4825285/1‑93 (MQ=255)

                                                   ttAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTatga                                      >  2:708267/1‑93 (MQ=255)

                                                     aGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTatgatg                                    <  1:2616275/93‑1 (MQ=255)

                                                      gACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTAtgatga                                   <  1:1115210/93‑1 (MQ=255)

                                                      gACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTAtgatga                                   <  2:9650838/93‑1 (MQ=255)

                                                      gACCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTAtgatga                                   >  1:7836383/1‑93 (MQ=255)

                                                       aCCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGaa                                  <  2:2530267/93‑1 (MQ=255)

                                                       aCCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGaa                                  <  2:8619592/93‑1 (MQ=255)

                                                       aCCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGaa                                  >  2:3212710/1‑93 (MQ=255)

                                                       aCCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGaa                                  <  2:7614045/93‑1 (MQ=255)

                                                       aCCGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGaa                                  <  2:6631224/93‑1 (MQ=255)

                                                        ccGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAAt                                 <  2:1353461/93‑1 (MQ=255)

                                                         cGCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATg                                <  1:7990241/93‑1 (MQ=255)

                                                          gCCACACAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGCACCATGCGAATACTATTCCTATGATCAATGc                               >  2:598975/1‑93 (MQ=255)

                                                          gCCACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGc                               >  2:4597992/1‑93 (MQ=255)

                                                           ccACAAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGcc                              >  1:7448523/1‑93 (MQ=255)

                                                             acaAAATCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGcccc                            >  2:3892477/1‑93 (MQ=255)

                                                               aaaaTCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCaa                          <  2:2768880/93‑1 (MQ=255)

                                                                  aTCCGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCa                       <  2:9219868/93‑1 (MQ=255)

                                                                    ccGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAgg                     <  2:3337877/93‑1 (MQ=255)

                                                                     cGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGt                    >  2:1233686/1‑93 (MQ=255)

                                                                     cGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGt                    >  1:4462118/1‑93 (MQ=255)

                                                                     cGGTGGTGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGt                    <  2:7215166/93‑1 (MQ=255)

                                                                      ggtggtGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGtt                   <  2:1069004/93‑1 (MQ=255)

                                                                      ggtggtGATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGtt                   <  1:6299175/93‑1 (MQ=255)

                                                                       gtggtgATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGttt                  >  2:4218508/1‑93 (MQ=255)

                                                                         ggtgATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGc                <  1:4308847/93‑1 (MQ=255)

                                                                          gtgATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCt               <  1:3040164/93‑1 (MQ=255)

                                                                          gtgATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCt               >  2:4369809/1‑93 (MQ=255)

                                                                          gtgATTAATACTTTCCGCACTGGATCAATGGTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCt               >  2:6293727/1‑93 (MQ=255)

                                                                           tgATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTg              <  2:2211554/93‑1 (MQ=255)

                                                                            gATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGt             >  1:2948408/1‑93 (MQ=255)

                                                                            gATTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCACAATCAGGTTTGCTGt             >  1:847069/1‑93 (MQ=255)

                                                                             aTTAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTc            >  1:6011108/1‑93 (MQ=255)

                                                                               tAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAg          >  1:4984556/1‑93 (MQ=255)

                                                                               tAATACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAg          <  1:452592/93‑1 (MQ=255)

                                                                                aaTACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGc         >  1:10034832/1‑93 (MQ=255)

                                                                                aaTACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGc         <  2:1637322/93‑1 (MQ=255)

                                                                                aaTACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGc         >  2:5807017/1‑93 (MQ=255)

                                                                                aaTACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTCAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGc         <  1:3796708/93‑1 (MQ=255)

                                                                                aaTACTTTCCGCACTCGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGc         >  1:642237/1‑93 (MQ=255)

                                                                                 aTACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCg        <  1:5622195/93‑1 (MQ=255)

                                                                                  tACTTTCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGt       >  1:660984/1‑93 (MQ=255)

                                                                                     tttCCGCACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTccc    >  2:8921547/1‑93 (MQ=255)

GGCCCTCTTGTAACAAAAACACCACTTCTGCCGAAAATGTTGAGAAGGTTGTTAGACCGCCACAAAATCCGGTGGTGATTAATACTTTCCACACTGGATCAATGTTCGTCATCCTGCTGAACCATGCGAATCCTATTCCTATGATGAATGCCCCAATCAGGTTTGCTGTCAGCGTCCCCA  >  NC_000913/656879‑657058

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: `‑`