BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	996,684	A→G	S18P (TCC→CCC)	ycbP ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	996,684	0	A	G	99.2%	376.3 / ‑4.8	117	S18P (TCC→CCC)	ycbP	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (56/60); total (56/61)
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.43e-01

Only 100 of 119 total aligned reads displayed.

GTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGAGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCTTTTTTTGTTAACA  >  NC_000913/996594‑996774

 tAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCt                                                                                         <  1:2003581/93‑1 (MQ=255)

  aGATCTTCCGGGGCGAATTTTTGCAGATTCCAGTGATAACCCTCTACACCCAGGCCATTTAGTTTTCCCCGCGCATATTCCAGCAGGGGGCTg                                                                                        <  2:696202/93‑1 (MQ=255)

  aGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTg                                                                                        <  2:617701/93‑1 (MQ=255)

  aGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTg                                                                                        <  2:4304308/93‑1 (MQ=255)

    aTCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGa                                                                                      <  2:2186261/93‑1 (MQ=255)

    aTCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGa                                                                                      <  2:2498674/93‑1 (MQ=255)

    aTCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGa                                                                                      <  2:3150769/93‑1 (MQ=255)

     tCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAg                                                                                     >  2:865524/1‑93 (MQ=255)

      cTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGc                                                                                    >  2:2672953/1‑93 (MQ=255)

       ttCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCg                                                                                   <  1:2489897/93‑1 (MQ=255)

         ccGGGGCGAAGTTTTGCAGATTCCAGTGATAAGCCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCgag                                                                                 <  1:3282161/93‑1 (MQ=255)

         ccGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCgag                                                                                 >  2:1026988/1‑93 (MQ=255)

          cGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCgaga                                                                                <  2:1004699/93‑1 (MQ=255)

            gggCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAg                                                                              <  1:929972/93‑1 (MQ=255)

              gCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGa                                                                            >  2:852639/1‑93 (MQ=255)

              gCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGa                                                                            >  1:2663367/1‑93 (MQ=255)

               cGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGaa                                                                           <  2:4652116/93‑1 (MQ=255)

                gAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGaaa                                                                          <  2:4323116/93‑1 (MQ=255)

                gAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGaaa                                                                          <  1:129438/93‑1 (MQ=255)

                gAAGTGTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGaaa                                                                          <  2:3815587/93‑1 (MQ=255)

                 aaGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAg                                                                         <  2:2691930/93‑1 (MQ=255)

                 aaGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAg                                                                         >  2:252000/1‑93 (MQ=255)

                   gTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCg                                                                       >  1:3476355/1‑93 (MQ=255)

                      ttGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAgg                                                                    <  1:2408204/93‑1 (MQ=255)

                       tGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGa                                                                   >  2:1827465/1‑93 (MQ=255)

                          aGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTa                                                                >  1:2625551/1‑93 (MQ=255)

                             ttCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTAccc                                                             >  1:879349/1‑93 (MQ=255)

                             ttCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTAccc                                                             <  2:2181700/93‑1 (MQ=255)

                                cAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTACCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGcc                                                          <  2:4433851/93‑1 (MQ=255)

                                 aGTGATAAGCCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCa                                                         <  1:1903139/93‑1 (MQ=255)

                                 aGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCa                                                         >  1:1446737/1‑93 (MQ=255)

                                  gTGATAACCCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAg                                                        <  1:3883940/93‑1 (MQ=255)

                                  gTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAg                                                        >  2:3489822/1‑93 (MQ=255)

                                  gTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAg                                                        <  1:3485877/93‑1 (MQ=255)

                                   tGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAgg                                                       >  1:3022399/1‑93 (MQ=255)

                                   tGATAAACCTCTACATCCAGGCCAGTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAATCGAGGACTCCCCGCCAgg                                                       <  2:685048/93‑1 (MQ=255)

                                    gATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAggg                                                      <  1:3779648/93‑1 (MQ=255)

                                      tAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGtg                                                    >  1:721928/1‑93 (MQ=255)

                                      tAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGtg                                                    >  1:1211288/1‑93 (MQ=255)

                                       aaaCCTCTTCCTCCAGGCCATTTAGTTTTTCCCGCGCATTTTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCCAGGACTACCCGCCAGGGtga                                                   <  2:4593739/93‑1 (MQ=255)

                                       aaaCCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGtga                                                   >  1:3872218/1‑93 (MQ=255)

                                       aaaCCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGtga                                                   <  1:545040/93‑1 (MQ=255)

                                         aCCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGtgatg                                                 <  1:2599678/93‑1 (MQ=255)

                                          ccTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGtgatga                                                >  2:947382/1‑93 (MQ=255)

                                            tctACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGaca                                              <  2:3152664/93‑1 (MQ=255)

                                            tctACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGaca                                              <  2:3152663/93‑1 (MQ=255)

                                            tctACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGAGAGGACTACCCGCCAGGGTGATGaca                                              >  1:188886/1‑93 (MQ=255)

                                             ctACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGacac                                             <  1:2101993/93‑1 (MQ=255)

                                              tACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACg                                            >  1:4324186/1‑93 (MQ=255)

                                              tACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACg                                            >  1:2616365/1‑93 (MQ=255)

                                              tACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACg                                            >  1:887801/1‑93 (MQ=255)

                                              tACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACg                                            >  1:4016069/1‑93 (MQ=255)

                                              tACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACg                                            >  2:1144678/1‑93 (MQ=255)

                                               aCATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGc                                           >  2:1952230/1‑93 (MQ=255)

                                                 aTCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCAt                                         >  1:1399211/1‑93 (MQ=255)

                                                 aTCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCAt                                         >  1:4607333/1‑93 (MQ=255)

                                                  tCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATa                                        >  2:4347960/1‑93 (MQ=255)

                                                  tCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATa                                        >  2:2939229/1‑93 (MQ=255)

                                                  tCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATa                                        >  1:4238011/1‑93 (MQ=255)

                                                  tCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATa                                        >  2:3805666/1‑93 (MQ=255)

                                                   ccAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATAc                                       <  2:1492025/93‑1 (MQ=255)

                                                   ccAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATAc                                       <  2:3789804/93‑1 (MQ=255)

                                                   ccAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATAc                                       <  2:3931371/93‑1 (MQ=255)

                                                   ccAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATAc                                       >  1:2942938/1‑93 (MQ=255)

                                                    cAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATAct                                      >  2:3130003/1‑93 (MQ=255)

                                                    cAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATAct                                      >  1:1925254/1‑93 (MQ=255)

                                                    cAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATAct                                      >  2:2870905/1‑93 (MQ=255)

                                                    cAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATAct                                      >  2:3697139/1‑93 (MQ=255)

                                                     aGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATActc                                     >  1:3580485/1‑93 (MQ=255)

                                                     aGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATActc                                     >  1:3580484/1‑93 (MQ=255)

                                                      ggCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATActct                                    >  2:2338788/1‑93 (MQ=255)

                                                      ggCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATActct                                    >  1:3438322/1‑93 (MQ=255)

                                                       gCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATActctc                                   <  1:1917487/93‑1 (MQ=255)

                                                         cATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACtctctc                                 >  1:2436597/1‑93 (MQ=255)

                                                          aTTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTcc                                >  1:2083217/1‑93 (MQ=255)

                                                           tttAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCt                               <  1:514481/93‑1 (MQ=255)

                                                           tttAGTTTTTCCCGCGCATATTACAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCt                               <  1:2046889/93‑1 (MQ=255)

                                                              aGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTtat                            <  1:2721939/93‑1 (MQ=255)

                                                               gtgtttCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTtata                           <  1:1321042/90‑1 (MQ=255)

                                                                tttttCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATaa                          <  1:2961208/93‑1 (MQ=255)

                                                                    tCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCaa                      <  1:2467921/93‑1 (MQ=255)

                                                                      ccGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAAtt                    <  1:2043100/93‑1 (MQ=255)

                                                                        gcgcATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGt                  <  2:1579765/93‑1 (MQ=255)

                                                                        gcgcATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCCCTCCTTATAACCAATTGt                  <  1:2348609/93‑1 (MQ=255)

                                                                         cgcATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGtt                 <  2:2229385/93‑1 (MQ=255)

                                                                          gcATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTc                >  2:2845431/1‑93 (MQ=255)

                                                                          gcATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTc                <  2:3947301/93‑1 (MQ=255)

                                                                          gcATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTc                >  1:1927193/1‑93 (MQ=255)

                                                                          gcATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTc                <  1:4533882/93‑1 (MQ=255)

                                                                            atatTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCtt              <  2:3379084/93‑1 (MQ=255)

                                                                               ttCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCttttt           >  1:1251442/1‑93 (MQ=255)

                                                                                tCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCtttttt          <  1:2979869/93‑1 (MQ=255)

                                                                                   aGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCTTTTTTTGt       <  2:4343877/93‑1 (MQ=255)

                                                                                      aaGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCTTTTTTTGTTaa    <  1:2829718/93‑1 (MQ=255)

                                                                                      aaGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGGGATGACACGCATACTCTCTCCTTATAACCAATTGTTCTTTTTTTGTTac    >  2:16511/1‑92 (MQ=255)

                                                                                      aaGGGGCTGGAGCGAAAAGGAAACCGACGAGTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTTTAACCAATTGTTCTTTTTTTGTTac    <  2:3016774/93‑2 (MQ=255)

                                                                                       aGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCTTTTTTTGTTAAc   >  2:72028/1‑93 (MQ=255)

                                                                                       aGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCTTTTTTTGTTAAc   <  2:1857248/93‑1 (MQ=255)

                                                                                       aGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGCGATGACACGCATACTCTCTCCTTATAACCAATTGTTCTTTTTTTGTTAAc   >  2:72027/1‑93 (MQ=255)

                                                                                        ggAGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCTTTTTTTGTTAACa  <  2:1575130/93‑1 (MQ=255)

GTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGAGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCTTTTTTTGTTAACA  >  NC_000913/996594‑996774

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