BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,789,986	T→C	D20G (GAC→GGC)	b1707 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,789,986	0	T	C	99.3%	469.5 / ‑4.5	147	D20G (GAC→GGC)	b1707	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (67/78); total (69/78)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

Only 100 of 154 total aligned reads displayed.

GCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTC  >  NC_000913/1789896‑1790074

gcTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGcct                                                                                        >  1:2874924/1‑93 (MQ=255)

 cTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGcctc                                                                                       >  1:2911420/1‑93 (MQ=255)

 cTATAACCCGACTCGTCGGGATACGAACAGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGcctc                                                                                       >  1:4127393/1‑93 (MQ=255)

    taccccGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTa                                                                                    <  2:662176/90‑1 (MQ=255)

    taACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTa                                                                                    <  2:662177/93‑1 (MQ=255)

     aaCCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTAt                                                                                   >  2:1839803/1‑93 (MQ=255)

     aaCCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTAt                                                                                   <  2:3277885/93‑1 (MQ=255)

     aaCCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTAt                                                                                   <  1:4420871/93‑1 (MQ=255)

      aCCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATc                                                                                  <  2:4583702/93‑1 (MQ=255)

          gACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCa                                                                              <  1:1081389/93‑1 (MQ=255)

            cTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAgg                                                                            >  2:1560975/1‑93 (MQ=255)

            cTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAgg                                                                            <  2:2676467/93‑1 (MQ=255)

             tcgtcgGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGa                                                                           >  2:765550/1‑93 (MQ=255)

             tcgtcgGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGa                                                                           >  2:2404427/1‑93 (MQ=255)

                tcgGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGCAAt                                                                        >  2:3950371/1‑93 (MQ=255)

                tcgGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAAt                                                                        <  1:718480/93‑1 (MQ=255)

                tcgGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAAt                                                                        >  1:3785774/1‑93 (MQ=255)

                tcgGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAAt                                                                        >  1:3941083/1‑93 (MQ=255)

                tcgGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAAt                                                                        <  2:413821/93‑1 (MQ=255)

                tcgGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAAt                                                                        >  1:1793690/1‑93 (MQ=255)

                  gggATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAg                                                                      <  2:2053988/93‑1 (MQ=255)

                   ggATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGc                                                                     <  1:2593877/93‑1 (MQ=255)

                   ggATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGc                                                                     <  1:3918235/93‑1 (MQ=255)

                     aTACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCaa                                                                   <  2:2081348/93‑1 (MQ=255)

                      tACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAAt                                                                  >  2:4139214/1‑93 (MQ=255)

                        cGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCa                                                                >  2:2386218/1‑93 (MQ=255)

                        cGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCa                                                                >  2:1969275/1‑93 (MQ=255)

                          aaCCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGa                                                              <  1:1675991/93‑1 (MQ=255)

                           aCCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGaa                                                             <  2:1394382/93‑1 (MQ=255)

                           aCCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGaa                                                             >  2:3738494/1‑93 (MQ=255)

                            ccGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAAt                                                            <  2:2007477/93‑1 (MQ=255)

                            ccGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAAt                                                            <  2:1895552/93‑1 (MQ=255)

                            ccGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAAt                                                            >  1:234415/1‑93 (MQ=255)

                              gTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGa                                                          <  2:2592649/93‑1 (MQ=255)

                               tACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGAt                                                         <  2:365432/93‑1 (MQ=255)

                                aCCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATa                                                        <  1:4657548/93‑1 (MQ=255)

                                 ccATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATaa                                                       >  1:3666519/1‑93 (MQ=255)

                                 ccATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATaa                                                       >  1:4269479/1‑93 (MQ=255)

                                  cATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATaaa                                                      >  2:4508905/1‑93 (MQ=255)

                                   aTCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAg                                                     <  1:447538/93‑1 (MQ=255)

                                    tCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGa                                                    <  2:3751497/93‑1 (MQ=255)

                                     cTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGAt                                                   <  2:2815161/93‑1 (MQ=255)

                                     cTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGAt                                                   <  1:1447098/93‑1 (MQ=255)

                                     cTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGAt                                                   >  1:2009938/1‑93 (MQ=255)

                                      ttCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGAtt                                                  >  2:1265366/1‑93 (MQ=255)

                                      ttCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGAtt                                                  >  2:4259080/1‑93 (MQ=255)

                                       tCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATtc                                                 >  1:4336322/1‑93 (MQ=255)

                                       tCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATtc                                                 <  1:850647/93‑1 (MQ=255)

                                            ggAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCaa                                            <  1:2997869/93‑1 (MQ=255)

                                             gagaAATGGGTAATTAACTCAACTCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCaaa                                           <  1:2290273/93‑1 (MQ=255)

                                              agaAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCaaaa                                          <  1:319540/93‑1 (MQ=255)

                                               gaAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCaaaaa                                         >  2:3272457/1‑93 (MQ=255)

                                               gaAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCaaaaa                                         >  2:4264552/1‑93 (MQ=255)

                                               gaAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCaaaaa                                         >  2:1600447/1‑93 (MQ=255)

                                                 aaTGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCaaaaaaa                                       <  1:1953683/93‑1 (MQ=255)

                                                 aaTGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCaaaaaaa                                       >  2:4091431/1‑93 (MQ=255)

                                                    gggTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCt                                    >  2:3600017/1‑93 (MQ=255)

                                                     ggTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCtt                                   <  1:2492803/93‑1 (MQ=255)

                                                     ggTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCtt                                   <  2:2687559/93‑1 (MQ=255)

                                                      gTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTc                                  <  1:108405/93‑1 (MQ=255)

                                                       tAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCa                                 <  1:569575/93‑1 (MQ=255)

                                                       tAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCa                                 <  2:3963734/93‑1 (MQ=255)

                                                       tAATTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCa                                 <  1:278398/93‑1 (MQ=255)

                                                        aaTTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCAt                                >  2:3580551/1‑93 (MQ=255)

                                                         aTTAACTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATg                               >  1:3220229/1‑93 (MQ=255)

                                                          ttAACTCAACGCCGACAAGGACTTGCTGGTTGACTCTTATCGGCAGGAAATAGGCATCACAATGATCCAAATTCTCCAAAAAAATCTTGATgg                              >  2:3689246/1‑93 (MQ=255)

                                                            aaCTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTa                            >  1:673821/1‑93 (MQ=255)

                                                            aaCTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTa                            <  2:3647878/93‑1 (MQ=255)

                                                             aCTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTaa                           >  1:3338647/1‑93 (MQ=255)

                                                              cTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCCTGGCAAc                          >  2:586453/1‑93 (MQ=255)

                                                              cTCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAAc                          >  2:2315718/1‑93 (MQ=255)

                                                               tCAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAAcc                         >  1:2867505/1‑93 (MQ=255)

                                                                cAACGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCg                        >  1:4161631/1‑93 (MQ=255)

                                                                 aaCGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGt                       <  1:2293554/93‑1 (MQ=255)

                                                                  aCGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTc                      >  2:1295235/1‑93 (MQ=255)

                                                                  aCGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTc                      >  1:1564590/1‑93 (MQ=255)

                                                                   cGCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTcc                     <  2:1719916/93‑1 (MQ=255)

                                                                    gaagACAAGGGCTTGCTGGTTGACTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAAGCTTCATGGTAACCGTccc                    <  2:1637401/90‑1 (MQ=255)

                                                                    gCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTccc                    <  1:4051496/93‑1 (MQ=255)

                                                                    gCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTccc                    <  2:1637397/93‑1 (MQ=255)

                                                                    gCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTccc                    >  1:3545697/1‑93 (MQ=255)

                                                                    gCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTccc                    <  2:3797781/93‑1 (MQ=255)

                                                                    gCCGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTccc                    <  2:473306/93‑1 (MQ=255)

                                                                     ccGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCct                   <  1:1774460/93‑1 (MQ=255)

                                                                      cGACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCctc                  <  1:802232/93‑1 (MQ=255)

                                                                       gACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCctct                 >  1:3392170/1‑93 (MQ=255)

                                                                       gACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCctct                 >  2:2730369/1‑93 (MQ=255)

                                                                       gACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCctct                 <  1:723654/93‑1 (MQ=255)

                                                                       gACAAGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATCATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCctct                 >  2:1447926/1‑93 (MQ=255)

                                                                          aaGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGa              >  1:1384338/1‑93 (MQ=255)

                                                                          aaGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGa              >  2:2161456/1‑93 (MQ=255)

                                                                          aaGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGa              >  2:2575548/1‑93 (MQ=255)

                                                                           aGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGaa             >  1:263160/1‑93 (MQ=255)

                                                                           aGGACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGaa             >  1:2603913/1‑93 (MQ=255)

                                                                            ggACTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAg            <  2:2272720/93‑1 (MQ=255)

                                                                               cTTGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGAt         <  1:377904/93‑1 (MQ=255)

                                                                                 tGCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGt       >  2:2928529/1‑93 (MQ=255)

                                                                                  gCTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGtt      <  1:2500567/93‑1 (MQ=255)

                                                                                   cTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGttt     <  1:3371645/93‑1 (MQ=255)

                                                                                   cTGGTTGCCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGttt     >  1:4062714/1‑93 (MQ=255)

GCTATAACCCGACTCGTCGGGATACGAACCGTACCATCTTCACTGGAGAAATGGGTAATTAACTCAACGCCGACAAGGACTTGCTGGTTGTCTCTTATCGGCAGGAAATAGCAATCAGAATGATAAAGATTCTCCAAAAAAATCTTCATGGTAACCGTCCCTCTCGAAGGATGTTTTTC  >  NC_000913/1789896‑1790074

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