BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,183,028	T→C	V70A (GTC→GCC)	b3042 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,183,028	0	T	C	99.3%	443.6 / ‑4.7	144	V70A (GTC→GCC)	b3042	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (71/72); total (71/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.36e-01

Only 100 of 150 total aligned reads displayed.

ATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGCCAG  >  NC_000913/3182938‑3183117

aTGTGGAGAAAAAAATCCGCCAACCCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGccc                                                                                         <  1:4460297/93‑1 (MQ=255)

  gtgGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTg                                                                                       >  2:4134545/1‑93 (MQ=255)

   tgGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCTAGCCCTGt                                                                                      <  2:4571626/93‑1 (MQ=255)

   tgGAGAAAAAAACCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTGGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGt                                                                                      <  2:2997399/93‑1 (MQ=255)

     aagaAAAAAACCCGCCAACCCCTACAAGCGCAGCGGCCGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTa                                                                                    <  2:2654008/92‑1 (MQ=255)

      agaAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTAc                                                                                   >  2:3274848/1‑93 (MQ=255)

      agaAAAAAATCCGCCAAACCCGACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTAc                                                                                   <  2:387371/93‑1 (MQ=255)

       gaAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACg                                                                                  >  1:2498179/1‑93 (MQ=255)

       aaaaaaaaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACg                                                                                  <  2:4316901/92‑1 (MQ=255)

        aaaaaaaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACgg                                                                                 >  1:1993557/1‑93 (MQ=255)

          aaaaaCCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGAc                                                                               <  1:226949/93‑1 (MQ=255)

           aaaaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACt                                                                              >  2:3719788/1‑93 (MQ=255)

            aaaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTc                                                                             >  1:4350651/1‑93 (MQ=255)

             aaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCg                                                                            >  2:4195733/1‑93 (MQ=255)

             aaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCg                                                                            >  2:1864467/1‑93 (MQ=255)

             aaTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCg                                                                            >  2:4002684/1‑93 (MQ=255)

              acccGCCAACCCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGt                                                                           <  2:946955/91‑1 (MQ=255)

              aTCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGt                                                                           <  1:3027177/93‑1 (MQ=255)

               tCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTg                                                                          >  1:1553804/1‑93 (MQ=255)

                  gCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGaaa                                                                       >  2:3814706/1‑93 (MQ=255)

                   ccAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGaaaa                                                                      >  1:1889257/1‑93 (MQ=255)

                       ccccTGCAGGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTg                                                                  <  2:3823100/92‑1 (MQ=255)

                       aCCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTg                                                                  >  1:1232838/1‑93 (MQ=255)

                        cccTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTgg                                                                 >  2:3069471/1‑93 (MQ=255)

                        cccAACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTgg                                                                 <  1:56177/93‑1 (MQ=255)

                          cTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGgcg                                                               >  1:667321/1‑93 (MQ=255)

                           tACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGgcgc                                                              >  2:400458/1‑93 (MQ=255)

                           tACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGgcgc                                                              <  2:3687163/93‑1 (MQ=255)

                           tACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGgcgc                                                              >  2:3606176/1‑93 (MQ=255)

                           tACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGgcgc                                                              >  1:4137486/1‑93 (MQ=255)

                           tACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGgcgc                                                              <  1:3747860/93‑1 (MQ=255)

                            aCAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCgct                                                             >  2:1606505/1‑93 (MQ=255)

                            aCAAGCCCAGCGGCCGCGCGTCGACCTTGTAAGCCGGGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCgct                                                             <  1:816972/93‑1 (MQ=255)

                              aaGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCgctgc                                                           <  1:4230609/93‑1 (MQ=255)

                               aGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCgctgct                                                          >  2:1967684/1‑93 (MQ=255)

                                gcgcAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCtt                                                         <  1:465613/93‑1 (MQ=255)

                                gcgcAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCtt                                                         <  1:3812273/93‑1 (MQ=255)

                                gcgcAGCTGACGCGCCTCGACCTTGTAACCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCtt                                                         <  1:1122040/93‑1 (MQ=255)

                                  gcAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGa                                                       >  2:4501349/1‑93 (MQ=255)

                                  gcAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGa                                                       >  2:2721080/1‑93 (MQ=255)

                                   cAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGaa                                                      <  1:2389363/93‑1 (MQ=255)

                                    cgCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAAc                                                     <  1:4179204/92‑1 (MQ=255)

                                    aGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAAc                                                     <  1:1609513/93‑1 (MQ=255)

                                    aGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAAc                                                     <  1:2623214/93‑1 (MQ=255)

                                      cTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACaa                                                   <  2:1709391/93‑1 (MQ=255)

                                       tGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACaaa                                                  >  1:1017413/1‑93 (MQ=255)

                                        gtcgcgCCTCGATCGCTTAAGCCGTGAAGATTACGACGTCGAACCGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAg                                                 <  2:2389788/91‑1 (MQ=255)

                                        cgcgcgCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAg                                                 <  2:2516453/91‑1 (MQ=255)

                                         aCGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGa                                                <  2:450218/93‑1 (MQ=255)

                                         aCGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGa                                                <  2:2295232/93‑1 (MQ=255)

                                            cgcCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATc                                             <  1:2058589/93‑1 (MQ=255)

                                            cgcCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATc                                             <  2:4304002/93‑1 (MQ=255)

                                            cgcCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATc                                             <  1:3961765/93‑1 (MQ=255)

                                            cgcCTCGACCTTGTAAGCCGTGAAGAGTTCCACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATc                                             <  2:787931/93‑1 (MQ=255)

                                             gcCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAGAACTGGCGCTGCTTGAACAAAGAATCa                                            <  1:1827335/93‑1 (MQ=255)

                                             gcCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCa                                            <  2:3040364/93‑1 (MQ=255)

                                             gcCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCa                                            <  1:1556472/93‑1 (MQ=255)

                                             gcCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCa                                            <  2:736203/93‑1 (MQ=255)

                                              ccTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAg                                           >  2:1799235/1‑93 (MQ=255)

                                              ccTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAg                                           <  1:3940785/93‑1 (MQ=255)

                                               cTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGc                                          >  1:1719821/1‑93 (MQ=255)

                                                tCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCg                                         <  1:4458578/93‑1 (MQ=255)

                                                   aCCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAAc                                      <  1:1800969/93‑1 (MQ=255)

                                                   aCCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAAc                                      >  2:3257795/1‑93 (MQ=255)

                                                   aCCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAAc                                      >  2:4044592/1‑93 (MQ=255)

                                                    ccTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACt                                     >  2:1929021/1‑93 (MQ=255)

                                                    ccTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACt                                     >  1:2911093/1‑93 (MQ=255)

                                                        gTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGaa                                 >  2:2973065/1‑93 (MQ=255)

                                                        gTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGaa                                 >  1:1399459/1‑93 (MQ=255)

                                                         tAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGaaa                                <  1:104166/93‑1 (MQ=255)

                                                         tAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGaaa                                <  2:3398766/93‑1 (MQ=255)

                                                         tAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGaaa                                >  1:4556417/1‑93 (MQ=255)

                                                           aGCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAAt                              <  2:93006/93‑1 (MQ=255)

                                                            gCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATc                             <  1:36355/93‑1 (MQ=255)

                                                            gCCGTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATc                             <  2:1908811/93‑1 (MQ=255)

                                                               gTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCa                          <  2:3493684/93‑1 (MQ=255)

                                                               gTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCa                          <  1:2609747/93‑1 (MQ=255)

                                                               gTGAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTAGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCa                          <  2:774556/93‑1 (MQ=255)

                                                                 gAAGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGt                        >  1:1254590/1‑93 (MQ=255)

                                                                  aaGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTa                       >  1:3120133/1‑93 (MQ=255)

                                                                  aaGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTa                       <  1:3274848/93‑1 (MQ=255)

                                                                  aaGAGTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTa                       >  2:243142/1‑93 (MQ=255)

                                                                   agagTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTAc                      >  1:2764545/1‑93 (MQ=255)

                                                                     agTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTg                    <  2:1170332/93‑1 (MQ=255)

                                                                     agTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTg                    >  2:472870/1‑93 (MQ=255)

                                                                     agTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTg                    <  2:2160598/93‑1 (MQ=255)

                                                                      gTTCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTga                   <  2:1358399/93‑1 (MQ=255)

                                                                      attCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTga                   <  2:4178880/92‑1 (MQ=255)

                                                                        tCGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTgaga                 >  2:1714800/1‑93 (MQ=255)

                                                                         cGACGTCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGAt                <  1:3363374/93‑1 (MQ=255)

                                                                          gACGTCCAAACGCAAGCCCTGTTACGGACGCGTGAAAAACTGGCGCTGCGTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATc               <  1:1610645/93‑1 (MQ=255)

                                                                              tCCAAACGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCaaaa           <  2:150349/93‑1 (MQ=255)

                                                                                 aaaCGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAagc        <  2:2712406/93‑1 (MQ=255)

                                                                                  aaCGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAagca       >  1:316068/1‑93 (MQ=255)

                                                                                   aCGCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAagcag      <  2:3745157/93‑1 (MQ=255)

                                                                                     gCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGcc    >  2:4251741/1‑93 (MQ=255)

                                                                                     gCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGcc    >  2:198624/1‑93 (MQ=255)

                                                                                     gCAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGcc    <  2:693449/93‑1 (MQ=255)

                                                                                      cAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGCCa   <  1:2323720/93‑1 (MQ=255)

                                                                                      cAAGCCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGCCa   >  2:1394591/1‑93 (MQ=255)

ATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGCGAACTGGAAAATCGCAGTACTGAGATCAAAAAGCAGCCAG  >  NC_000913/3182938‑3183117

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