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.7%	909.1 / ‑5.2	290	V258A (GTA→GCA)	caiC	putative crotonobetain
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (144/145); total (145/145)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-01

Only 100 of 304 total aligned reads displayed.

TGTTGCTGATCGTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGG  >  NC_000913/36962‑37141

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

 gtTGCTGATCGTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCt                                                                                        >  2:701278/1‑93 (MQ=255)

  ttGCTGATCGTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTg                                                                                       >  1:1001608/1‑93 (MQ=255)

  ttGCTGAGCGTTCTCGTAAGGCTGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGTTAACGGTGGCGCGGTACTTCTGTGCCTg                                                                                       <  2:3114698/93‑1 (MQ=255)

     cTGATCGTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTcc                                                                                    >  2:4920512/1‑93 (MQ=255)

     cTGATCGTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTcc                                                                                    >  2:5384382/1‑93 (MQ=255)

     cTGATCGTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTcc                                                                                    <  2:6888980/93‑1 (MQ=255)

      tGATCGTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTccc                                                                                   <  1:5092668/93‑1 (MQ=255)

      tGATCGTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTccc                                                                                   <  2:7768784/93‑1 (MQ=255)

       gATCGTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTcccc                                                                                  >  1:2309044/1‑93 (MQ=255)

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

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

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

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

          cGTTCGCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGa                                                                               <  1:94739/93‑1 (MQ=255)

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

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

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

               gCTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACGTCTGTGCCTGTCCCCAGAAGgcg                                                                          <  2:7652562/93‑1 (MQ=255)

                cTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGgcgc                                                                         >  1:2911981/1‑93 (MQ=255)

                cTGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGgcgc                                                                         <  1:5234737/93‑1 (MQ=255)

                 tGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGcgcg                                                                        <  1:549355/93‑1 (MQ=255)

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

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

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

                 tGAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGCGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGcgcg                                                                        >  1:1728048/1‑93 (MQ=255)

                  gAAGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCgg                                                                       >  2:7873617/1‑93 (MQ=255)

                   aaGGCGGCTGTCCCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCggg                                                                      <  2:3591053/93‑1 (MQ=255)

                   aaGGCGGCTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCggg                                                                      >  1:1852081/1‑93 (MQ=255)

                     ggcggcTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGgcg                                                                    >  2:3404126/1‑93 (MQ=255)

                     ggcggcTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGgcg                                                                    <  1:1269536/93‑1 (MQ=255)

                       gggcTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCt                                                                  <  1:5930981/92‑1 (MQ=255)

                       cggcTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCt                                                                  <  1:3535821/93‑1 (MQ=255)

                       cggcTGTACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCt                                                                  <  1:6844624/93‑1 (MQ=255)

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

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

                            gTACCATCAACGTACGGGTCATCATCGGGAAACATTCGGGAAAGGGGGGGCCGGACTTCTTTGGCCGTCCCCCGAAAGCGCGGGGGCCGtttt                                                             >  2:5696529/1‑93 (MQ=255)

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

                             tACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTAttt                                                            >  2:4765959/1‑93 (MQ=255)

                             tACCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTAttt                                                            >  1:2789803/1‑93 (MQ=255)

                              aCCATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTtc                                                           >  1:205204/1‑93 (MQ=255)

                                cATCAACGTACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTtctc                                                         <  1:2559276/93‑1 (MQ=255)

                                  tCAACGTACTGTTCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGa                                                       <  1:5629523/93‑1 (MQ=255)

                                  tCAACGTACGGAGCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGCGCGCTGTATGTCTCGa                                                       >  2:7923581/1‑93 (MQ=255)

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

                                        tACGGATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGca                                                 <  2:2998761/93‑1 (MQ=255)

                                         actgATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGcac                                                <  1:5722109/90‑1 (MQ=255)

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

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

                                            gATCATCATCGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACaaa                                             <  1:6357931/93‑1 (MQ=255)

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

                                                 tcatcGGAATACATTCGGTATCGGTGGCGCGGTACTTCTGTGCCTGTCGCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTgg                                        <  2:3748191/93‑1 (MQ=255)

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

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

                                                  catcGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGc                                       <  2:5426746/93‑1 (MQ=255)

                                                  catcGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGc                                       <  2:5792248/93‑1 (MQ=255)

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

                                                    tcGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGccc                                     >  1:3979301/1‑93 (MQ=255)

                                                     cGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGcccc                                    <  2:1559049/93‑1 (MQ=255)

                                                      ggAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCg                                   >  2:2502803/1‑93 (MQ=255)

                                                      ggAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCg                                   >  2:4382371/1‑93 (MQ=255)

                                                      ggAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCg                                   >  2:3239171/1‑93 (MQ=255)

                                                      ggAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCg                                   >  1:7213605/1‑93 (MQ=255)

                                                       gAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCgg                                  >  2:4991913/1‑93 (MQ=255)

                                                       gAATACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCgc                                  >  2:1484293/1‑92 (MQ=255)

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

                                                         aTACATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCa                                <  1:2186690/93‑1 (MQ=255)

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

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

                                                            cATTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGaa                             >  1:3055489/1‑93 (MQ=255)

                                                             aTTCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGaaa                            <  2:2961560/93‑1 (MQ=255)

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

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

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

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

                                                               tCGGTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAAc                          <  1:6227253/93‑1 (MQ=255)

                                                                  gTAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAAcgcc                       >  2:6265901/1‑93 (MQ=255)

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

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

                                                                  gTAACGGTGGCGCGGTACTTCTGTGCCTGACCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAActcc                       <  1:1353627/93‑1 (MQ=255)

                                                                   tAACGGTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAAcgccc                      >  2:4029764/1‑92 (MQ=255)

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

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

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

                                                                       ggTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCAt                  <  2:2390727/93‑1 (MQ=255)

                                                                        gTGGCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATc                 <  1:2420074/93‑1 (MQ=255)

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

                                                                          ggCGCGGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATcgc               >  2:3283375/1‑93 (MQ=255)

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

                                                                             gcgGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATcgccgc            <  1:2455883/93‑1 (MQ=255)

                                                                             gcgGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATcgccgc            <  1:1658440/93‑1 (MQ=255)

                                                                              cgGTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCa           >  2:2065615/1‑93 (MQ=255)

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

                                                                               ggTACTTCTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAg          <  1:457895/93‑1 (MQ=255)

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

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

                                                                                    ttCTGTGCCGGCCCCCAGAAGGCGCGGGCGCTGTATCTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTacac     <  2:5622536/93‑1 (MQ=255)

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

                                                                                      cTGTGCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTg   <  2:7222772/93‑1 (MQ=255)

                                                                                       tgtgCCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTgg  >  1:4732059/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: `‑`