BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	37,052	A→G	V258A (GTA→GCA)	caiC ←	putative crotonobetain

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	37,052	0	A	G	99.1%	676.3 / ‑3.7	215	V258A (GTA→GCA)	caiC	putative crotonobetain
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (121/92); total (122/93)
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.78e-01

Only 100 of 222 total aligned reads displayed.

TGTTGCTGATCGTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGG  >  NC_000913/36962‑37141

tgtTGCTGATCGTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGcc                                                                                         <  2:2158709/93‑1 (MQ=255)

 gtTGCTGATCGTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCt                                                                                        <  1:3070329/93‑1 (MQ=255)

 gtTGCTGATCGTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCt                                                                                        <  1:1152243/93‑1 (MQ=255)

  ttGCTGATCGTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTg                                                                                       <  2:2888580/93‑1 (MQ=255)

  ttGCTGATCGTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTg                                                                                       <  1:360023/93‑1 (MQ=255)

  gtGCTGATCGTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTg                                                                                       <  1:3693828/92‑1 (MQ=255)

      tGATCGTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTccc                                                                                   >  2:3956982/1‑93 (MQ=255)

      tGATCGTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTccc                                                                                   >  1:2113623/1‑93 (MQ=255)

        aTCGTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCa                                                                                 >  1:4538477/1‑93 (MQ=255)

        aTCGTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCa                                                                                 <  1:504778/93‑1 (MQ=255)

        aTCGTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCa                                                                                 >  2:1234707/1‑93 (MQ=255)

           gTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGaa                                                                              <  2:1606336/93‑1 (MQ=255)

           gTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGaa                                                                              <  1:3623877/93‑1 (MQ=255)

           gTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGaa                                                                              <  1:2729638/93‑1 (MQ=255)

           gTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGaa                                                                              <  1:4178771/93‑1 (MQ=255)

           gTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGaa                                                                              <  1:4582076/93‑1 (MQ=255)

            ttCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAg                                                                             <  2:438644/93‑1 (MQ=255)

            ttCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAg                                                                             >  1:1976813/1‑93 (MQ=255)

            ttCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAg                                                                             >  2:2084945/1‑93 (MQ=255)

             tCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAgg                                                                            <  1:3328264/93‑1 (MQ=255)

              cGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGgc                                                                           <  2:3088556/93‑1 (MQ=255)

              cGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGgc                                                                           <  1:472261/93‑1 (MQ=255)

               gCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGgcg                                                                          >  1:3219509/1‑93 (MQ=255)

               gCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGgcg                                                                          >  2:2648626/1‑93 (MQ=255)

               gCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGgcg                                                                          >  1:877654/1‑93 (MQ=255)

                 tGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGcgcg                                                                        >  1:2063019/1‑93 (MQ=255)

                  gAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCgg                                                                       <  1:3031110/93‑1 (MQ=255)

                  gAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCgg                                                                       <  2:1671100/93‑1 (MQ=255)

                    aGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCgggg                                                                     >  2:2669277/1‑92 (MQ=255)

                    aGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGgc                                                                     >  2:3027752/1‑93 (MQ=255)

                      gcggcTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGgcgc                                                                   >  1:2883931/1‑93 (MQ=255)

                        ggcTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTg                                                                 <  1:1166522/93‑1 (MQ=255)

                         gcTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGt                                                                >  2:2458295/1‑93 (MQ=255)

                         gcTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGt                                                                >  2:1065894/1‑93 (MQ=255)

                          cTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTa                                                               >  2:3272697/1‑93 (MQ=255)

                          cTGTACCATCAACGTACGGATCATCAAAGGAATACATTCGGTTTAGGTGACGCGGGCCTTCTGTGACTGTCTCCAGAAAGCGCGGGCACTggt                                                               >  2:3402825/1‑91 (MQ=255)

                           tGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTAt                                                              >  2:392917/1‑93 (MQ=255)

                            gTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTAtt                                                             <  1:76774/93‑1 (MQ=255)

                                 aTCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCg                                                        >  1:4097771/1‑93 (MQ=255)

                                 aTCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCg                                                        >  1:4441807/1‑93 (MQ=255)

                                 aTCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCg                                                        >  1:1076933/1‑93 (MQ=255)

                                  tCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCgc                                                       >  2:3620818/1‑92 (MQ=255)

                                  tCAACGTACGCATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGa                                                       <  2:4613029/93‑1 (MQ=255)

                                  tCAACGGACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCgc                                                       >  2:1108680/1‑92 (MQ=255)

                                   cAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGAc                                                      <  2:928484/93‑1 (MQ=255)

                                     aCGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCa                                                    <  1:631280/93‑1 (MQ=255)

                                     aCGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCa                                                    <  2:1800337/93‑1 (MQ=255)

                                     aCGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCa                                                    <  2:4159716/93‑1 (MQ=255)

                                       gTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGc                                                  >  2:2011766/1‑93 (MQ=255)

                                       gTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGc                                                  >  1:1681831/1‑93 (MQ=255)

                                       gTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGc                                                  >  2:2788199/1‑93 (MQ=255)

                                          cGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGcaca                                               <  2:3176449/93‑1 (MQ=255)

                                          cGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGcaca                                               <  1:2591929/93‑1 (MQ=255)

                                          cGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGcaca                                               >  1:1407256/1‑93 (MQ=255)

                                           ggATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACaa                                              >  2:1460485/1‑93 (MQ=255)

                                            gATCATCATCGGAATATATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACaaa                                             <  2:3017377/93‑1 (MQ=255)

                                            gATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACaaa                                             >  2:1038923/1‑93 (MQ=255)

                                            gATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACaaa                                             >  2:46951/1‑93 (MQ=255)

                                             atcatcatcGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAg                                            <  2:592363/93‑1 (MQ=255)

                                              tcatcatcGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAgg                                           <  2:343639/93‑1 (MQ=255)

                                              tcatcatcGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAgg                                           <  1:12382/93‑1 (MQ=255)

                                                atcatcGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTg                                         >  1:4342532/1‑93 (MQ=255)

                                                atcatcGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTg                                         <  1:813748/93‑1 (MQ=255)

                                                 tcatcGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTgg                                        >  1:1049043/1‑93 (MQ=255)

                                                 tcatcGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTgg                                        >  1:855832/1‑93 (MQ=255)

                                                 tcatcGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTgg                                        >  1:3348324/1‑93 (MQ=255)

                                                  catcGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGc                                       >  1:2428441/1‑93 (MQ=255)

                                                   atcGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGcc                                      >  2:2404292/1‑93 (MQ=255)

                                                   atcGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGcc                                      >  1:1938691/1‑93 (MQ=255)

                                                        aaTACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGc                                 >  2:1894203/1‑93 (MQ=255)

                                                         aTACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCa                                >  2:962267/1‑93 (MQ=255)

                                                           aCATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGa                              >  1:644035/1‑93 (MQ=255)

                                                            cATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGaa                             <  1:1728473/93‑1 (MQ=255)

                                                              ttCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGaaaa                           <  2:1439522/93‑1 (MQ=255)

                                                              ttCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGaaaa                           >  2:3692630/1‑93 (MQ=255)

                                                              ttCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGAAGaaaa                           >  2:2541322/1‑93 (MQ=255)

                                                                 ggTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAAcgc                        >  1:3552309/1‑93 (MQ=255)

                                                                 ggTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAAcgc                        <  1:1609496/93‑1 (MQ=255)

                                                                 ggTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAAcgc                        >  2:1498097/1‑93 (MQ=255)

                                                                 ggTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAAcgc                        >  2:2995166/1‑93 (MQ=255)

                                                                 ggTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAAcgc                        <  2:1308047/93‑1 (MQ=255)

                                                                  gTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAAcgcc                       <  2:3446566/93‑1 (MQ=255)

                                                                    aaCGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAAcgccgc                     >  1:249552/1‑93 (MQ=255)

                                                                     aCGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACgccgcc                    >  1:462464/1‑93 (MQ=255)

                                                                         tGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATcg                >  1:3492701/1‑93 (MQ=255)

                                                                          ggCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATcgc               <  1:2991373/93‑1 (MQ=255)

                                                                            cgcgGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATcgccg             <  1:3646443/93‑1 (MQ=255)

                                                                            cgcgGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATGTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATcgccg             <  2:2694575/93‑1 (MQ=255)

                                                                             gcgGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATcgccgc            >  1:4347150/1‑93 (MQ=255)

                                                                               ggTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAg          >  2:1835147/1‑93 (MQ=255)

                                                                               ggTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAg          >  1:1016680/1‑93 (MQ=255)

                                                                                gTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGt         <  1:693682/93‑1 (MQ=255)

                                                                                 tACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCCCCGCAGTa        <  2:4584395/93‑1 (MQ=255)

                                                                                  aCTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTac       >  1:2947709/1‑93 (MQ=255)

                                                                                  aCTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTac       >  2:750719/1‑93 (MQ=255)

                                                                                    ttCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTacac     >  2:688514/1‑93 (MQ=255)

                                                                                    ttCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTacac     <  2:3368724/93‑1 (MQ=255)

                                                                                    ttCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTacac     <  1:4294709/93‑1 (MQ=255)

                                                                                     tCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACt    <  1:546095/93‑1 (MQ=255)

                                                                                       tgtgCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTgg  >  2:4274790/1‑93 (MQ=255)

TGTTGCTGATCGTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGG  >  NC_000913/36962‑37141

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