BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,360,184	G→A	D84N (GAC→AAC)	ycjC →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,360,184	0	G	A	96.9%	560.5 / 1.9	162	D84N (GAC→AAC)	ycjC	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (2/3); new base A (92/65); total (94/68)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.66e-01

Only 100 of 170 total aligned reads displayed.

AAAGCTGCTGAAGGTGTATGGTCTGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACGACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTGGCGATGATCTTTGAAACG  >  NC_000913/1360093‑1360273

aaaGCTGCTGAAGGTGTATGGTCTGTCACTCTCGGAATTCTTTTCCGAGCCGG‑AAAACCTGATGAGCCGCAGGTCGTCATTAATCAGgacaac                                                                                         <  1:1052097/93‑1 (MQ=255)

 aaGCTGCTGAAGGTGTATGGTCTGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGgacaac                                                                                         <  2:1742031/93‑1 (MQ=255)

 aaGCTGCTGAAGGTGTATGGTCTGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGgacaac                                                                                         >  1:116244/1‑93 (MQ=255)

      gctgAAGGTGTATGGTCTGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAAt                                                                                    <  1:3798054/93‑1 (MQ=255)

      gctgAAGGTGTATGGTCTGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAAt                                                                                    >  2:2568706/1‑93 (MQ=255)

      gctgAAGGGGTATGGTCTGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAAt                                                                                    <  2:3271954/93‑1 (MQ=255)

       ctgAAGGTGTATGGTCTGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAAtt                                                                                   >  2:517524/1‑93 (MQ=255)

       ctgAAGGTGTATGGTCTGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAAtt                                                                                   >  2:300629/1‑93 (MQ=255)

        tgAAGGTGTATGGTCTGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTg                                                                                  >  2:2829940/1‑93 (MQ=255)

        tgAAGGTGTATGGTCTGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTg                                                                                  >  2:3774718/1‑93 (MQ=255)

          aaGGTGTATGGTCTGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTgag                                                                                >  1:4732896/1‑93 (MQ=255)

           aGGTGTATGGTCTGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTgaga                                                                               >  2:668075/1‑93 (MQ=255)

           aGGTGTATGGTCTGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTgaga                                                                               >  1:514227/1‑93 (MQ=255)

            ggTGTATGGTCTGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGAt                                                                              >  1:4634369/1‑93 (MQ=255)

            ggTGTATGGTCTGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGAt                                                                              >  2:2835031/1‑93 (MQ=255)

             gtgtATGGTCTGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATg                                                                             >  2:3482176/1‑93 (MQ=255)

              tgtATGGTCTGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATgg                                                                            <  2:536711/93‑1 (MQ=255)

                 aTGGTCTGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTa                                                                         <  1:441553/93‑1 (MQ=255)

                  tGGTCTGTCACTCTCGGAATTCTTTTCCGAGCCGGCAAAACCTGATTATCCGCAGGTCGCCTATATTCAGGATAACTTAATTGAGTTGGGTat                                                                        >  1:2917428/1‑92 (MQ=255)

                  tGGTCTGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACGACTTAATTGAGATGGGTAg                                                                        <  2:4271544/93‑1 (MQ=255)

                  tGGTCTGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAg                                                                        <  1:2481539/93‑1 (MQ=255)

                   ggTCTGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGt                                                                       <  1:2825712/93‑1 (MQ=255)

                    gTCTGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTc                                                                      >  1:3537288/1‑93 (MQ=255)

                      cTGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAg                                                                    >  1:1504402/1‑93 (MQ=255)

                       tGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAgg                                                                   >  2:1174832/1‑93 (MQ=255)

                       tGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAgg                                                                   >  2:2467042/1‑93 (MQ=255)

                           cctctCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGgtgt                                                               <  1:3329856/92‑1 (MQ=255)

                           aCTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGgtgt                                                               >  2:790029/1‑93 (MQ=255)

                            ctctCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGgtgtg                                                              <  2:13936/93‑1 (MQ=255)

                             tctcGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGgtgtgt                                                             >  2:1954238/1‑93 (MQ=255)

                               tcGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCa                                                           <  1:740344/93‑1 (MQ=255)

                                cGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCaa                                                          <  1:391856/93‑1 (MQ=255)

                                 ggAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAAt                                                         <  1:3340563/93‑1 (MQ=255)

                                 ggAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAAt                                                         >  1:2281782/1‑93 (MQ=255)

                                 ggAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAAt                                                         >  2:1697132/1‑93 (MQ=255)

                                 ggAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGGGTGTCAAt                                                         >  1:870779/1‑93 (MQ=255)

                                 ggAATTCTTTACGGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCGAt                                                         <  2:2688910/93‑1 (MQ=255)

                                  gAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATg                                                        >  2:3491987/1‑93 (MQ=255)

                                   aaTTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGa                                                       >  1:3706581/1‑93 (MQ=255)

                                   aaTTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGa                                                       <  1:2726633/93‑1 (MQ=255)

                                       cTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCt                                                   >  2:2959317/1‑93 (MQ=255)

                                         tttCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTgg                                                 >  2:4682622/1‑93 (MQ=255)

                                          ttCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGt                                                <  2:2115123/93‑1 (MQ=255)

                                          ttCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGt                                                >  2:3720530/1‑93 (MQ=255)

                                          ttCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGt                                                <  1:1064227/93‑1 (MQ=255)

                                          ttCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGt                                                <  1:1124310/93‑1 (MQ=255)

                                           tCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGtt                                               <  1:304000/93‑1 (MQ=255)

                                            ccGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTc                                              <  2:2317074/93‑1 (MQ=255)

                                            ccGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTc                                              <  2:4624373/93‑1 (MQ=255)

                                             cGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCa                                             >  1:1112157/1‑93 (MQ=255)

                                               aGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATa                                           >  1:1641916/1‑93 (MQ=255)

                                                gCCGGAAAAGCCTGATGAGCCGCAGGTCGTCATTAAGCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATaa                                          <  2:2548559/93‑1 (MQ=255)

                                                gCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATaa                                          <  2:3408499/93‑1 (MQ=255)

                                                gCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATaa                                          <  1:2795896/93‑1 (MQ=255)

                                                gCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATaa                                          >  1:1254667/1‑93 (MQ=255)

                                                 ccGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAAc                                         <  1:845455/93‑1 (MQ=255)

                                                 ccGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAAc                                         <  2:3927162/93‑1 (MQ=255)

                                                   ggAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACgg                                       >  2:1832209/1‑93 (MQ=255)

                                                    gAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGt                                      >  1:937274/1‑93 (MQ=255)

                                                     aaaaaCCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTa                                     >  1:2702678/1‑93 (MQ=255)

                                                       aaaCCTGATGAGCCGCAGGTCGTCATTAATCAGGACGACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAAc                                   <  2:1588937/93‑1 (MQ=255)

                                                       aaaCCTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAAc                                   >  2:3746262/1‑93 (MQ=255)

                                                        aaCCTGATGAGCCGCAGGTCGTCATTAATCAGGACGACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAAcc                                  >  2:4026790/1‑93 (MQ=255)

                                                          ccTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCg                                <  2:448431/93‑1 (MQ=255)

                                                           cTGATGAGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGa                               <  2:2121124/93‑1 (MQ=255)

                                                            tgatgaGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGaa                              <  2:2116551/93‑1 (MQ=255)

                                                            tgatgaGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGaa                              >  1:3173490/1‑93 (MQ=255)

                                                             gatgaGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAAt                             >  1:4497962/1‑93 (MQ=255)

                                                             gatgaGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAAt                             >  1:2385585/1‑93 (MQ=255)

                                                              atgaGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATc                            <  2:4185855/93‑1 (MQ=255)

                                                                gaGCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGAGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGc                          >  2:3023161/1‑93 (MQ=255)

                                                                  gCCGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCAc                        <  1:942237/93‑1 (MQ=255)

                                                                   ccGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACg                       <  1:3510971/93‑1 (MQ=255)

                                                                    cGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGc                      <  1:1927386/93‑1 (MQ=255)

                                                                    cGCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGc                      <  1:4589120/93‑1 (MQ=255)

                                                                     gCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCt                     <  2:1924228/93‑1 (MQ=255)

                                                                     gCAGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCt                     >  1:414900/1‑93 (MQ=255)

                                                                       aGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTgg                   >  2:4582210/1‑93 (MQ=255)

                                                                       aGGTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTgg                   >  2:277500/1‑93 (MQ=255)

                                                                        ggTCGTCATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTGGc                  >  1:3075659/1‑93 (MQ=255)

                                                                         gtcgtcATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTcgcg                 >  1:3416533/1‑93 (MQ=255)

                                                                          tcgtcATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTGGCga                <  1:2460034/93‑1 (MQ=255)

                                                                          tcgtcATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTGGCga                >  1:3382021/1‑93 (MQ=255)

                                                                           tgttaTTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTGGCgat               <  2:1314885/89‑1 (MQ=255)

                                                                           cgtcATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTGGCgat               <  2:1314884/93‑1 (MQ=255)

                                                                            gtcATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTGGCgatg              >  1:245930/1‑93 (MQ=255)

                                                                            gtcATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTGGCgatg              >  1:1282246/1‑93 (MQ=255)

                                                                            gtcATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTGGCgatg              >  2:4107923/1‑93 (MQ=255)

                                                                              cATTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTGGCgatgat            >  1:651828/1‑93 (MQ=255)

                                                                               aTTAATCAGGACGACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTGGCGATGATc           >  2:1985618/1‑93 (MQ=255)

                                                                               aTTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTGGCGATGATc           >  1:4303806/1‑93 (MQ=255)

                                                                               aTTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTGGCGATGATc           >  2:3113264/1‑93 (MQ=255)

                                                                               aTTAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGCCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTGGCGATGATc           <  1:1951753/93‑1 (MQ=255)

                                                                                ttAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTGGCGATGATCt          >  1:2060867/1‑93 (MQ=255)

                                                                                 tAATCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTGGCGATGATCtt         <  1:1183798/93‑1 (MQ=255)

                                                                                   aTCAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTGGCGATGATCTTTg       >  2:2276761/1‑93 (MQ=255)

                                                                                     cAGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTGGCGATGATCTTTGaa     >  1:4000807/1‑93 (MQ=255)

                                                                                      aGGACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTGGCGATGATCTTTGaaa    >  1:4247169/1‑93 (MQ=255)

                                                                                       ggACAACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTGGCGATGATCTTTGAAAc   <  1:4548344/93‑1 (MQ=255)

                                                                                        gacaacTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTGGCGATGATCTTTGaaaag  >  2:3911357/1‑91 (MQ=255)

AAAGCTGCTGAAGGTGTATGGTCTGTCACTCTCGGAATTCTTTTCCGAGCCGGAAAAACCTGATGAGCCGCAGGTCGTCATTAATCAGGACGACTTAATTGAGATGGGTAGTCAGGGTGTGTCAATGAAGCTGGTTCATAACGGTAACCCGAATCGCACGCTGGCGATGATCTTTGAAACG  >  NC_000913/1360093‑1360273

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