BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,604,832	A→G	R146R (CGT→CGC)	yjjT ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,604,832	0	A	G	99.4%	555.9 / ‑4.6	178	R146R (CGT→CGC)	yjjT	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (93/84); total (93/85)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.05e-01

Only 100 of 193 total aligned reads displayed.

AAACACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGACGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCACGCC  >  NC_000913/4604742‑4604921

aaaCACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAggcg                                                                                         >  1:3554393/1‑93 (MQ=255)

 aaCACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAggcga                                                                                        >  1:1305637/1‑93 (MQ=255)

 aaCACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAggcga                                                                                        >  1:1641635/1‑93 (MQ=255)

        cAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAAt                                                                                 <  1:1515990/93‑1 (MQ=255)

        cAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAAt                                                                                 <  2:4144683/93‑1 (MQ=255)

        cAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAAt                                                                                 >  2:190949/1‑93 (MQ=255)

         aGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATa                                                                                <  1:1352098/93‑1 (MQ=255)

          ggCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATag                                                                               >  1:4486249/1‑93 (MQ=255)

            cAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATagag                                                                             >  1:4421018/1‑93 (MQ=255)

             aGCGTTTTGGCCGTCAGGCCATCGACGCTGTATTCACCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAgg                                                                            <  2:350764/93‑1 (MQ=255)

             aGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAgg                                                                            <  1:3674761/93‑1 (MQ=255)

              gCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGc                                                                           >  2:529963/1‑93 (MQ=255)

                gTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCca                                                                         <  1:316975/93‑1 (MQ=255)

                gTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCca                                                                         >  1:3582397/1‑93 (MQ=255)

                 ttttGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCcac                                                                        >  2:285967/1‑93 (MQ=255)

                 ttttGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCcaa                                                                        >  2:1231584/1‑92 (MQ=255)

                  tttGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCcaca                                                                       >  1:3723509/1‑93 (MQ=255)

                    tGACCGTCATGCCATCGACGCTGTATGCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGc                                                                     <  2:1780799/93‑1 (MQ=255)

                    tGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGc                                                                     <  1:2652806/93‑1 (MQ=255)

                      aCCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGcga                                                                   >  1:4308789/1‑93 (MQ=255)

                       ccGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGcgac                                                                  >  1:4567173/1‑93 (MQ=255)

                       ccGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGcgac                                                                  >  2:2816978/1‑93 (MQ=255)

                       ccGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGcgac                                                                  <  1:757501/93‑1 (MQ=255)

                         gTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCTAATACCGGCTGTTTTTCCATGCGACCAAAATAGAGGCCACAGcgacga                                                                >  2:4208748/1‑93 (MQ=255)

                         gTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGcgacga                                                                >  2:243923/1‑93 (MQ=255)

                          tCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAg                                                               <  2:122200/93‑1 (MQ=255)

                           cAGTCCATCCACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAgc                                                              <  1:2643923/93‑1 (MQ=255)

                            aGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAgcg                                                             >  2:1623472/1‑93 (MQ=255)

                             ggCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAgcgc                                                            <  1:2583864/93‑1 (MQ=255)

                             ggCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAgcgc                                                            <  2:4405588/93‑1 (MQ=255)

                               ccATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTg                                                          >  1:4624622/1‑93 (MQ=255)

                                cATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGt                                                         >  2:3523571/1‑93 (MQ=255)

                                 aTCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTc                                                        <  1:31860/93‑1 (MQ=255)

                                 aTCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTc                                                        <  2:2772868/93‑1 (MQ=255)

                                 aTCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTc                                                        >  2:3591564/1‑93 (MQ=255)

                                 aTCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTc                                                        >  2:389395/1‑93 (MQ=255)

                                   cGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGa                                                      >  1:1061053/1‑93 (MQ=255)

                                    gACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGAc                                                     >  2:2340569/1‑93 (MQ=255)

                                    gACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGAc                                                     >  1:3919058/1‑93 (MQ=255)

                                    gACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTACCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGAc                                                     <  2:3639244/93‑1 (MQ=255)

                                     aCGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACt                                                    <  1:648960/93‑1 (MQ=255)

                                     aCGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACt                                                    <  1:4450135/93‑1 (MQ=255)

                                       gCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTTCAGGCGACCAAAATAGAGGCCACAGCGACGACCGCGGCTGACttt                                                  >  2:2993086/1‑93 (MQ=255)

                                       gCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACttt                                                  >  1:3424396/1‑93 (MQ=255)

                                        cTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTa                                                 >  2:895629/1‑93 (MQ=255)

                                        cTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTa                                                 >  2:4105578/1‑93 (MQ=255)

                                        cTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTa                                                 >  1:4473917/1‑93 (MQ=255)

                                        cTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTa                                                 <  1:2607990/93‑1 (MQ=255)

                                          gTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTAtt                                               >  2:3920301/1‑93 (MQ=255)

                                            aTTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCa                                             <  2:4160039/93‑1 (MQ=255)

                                            aTTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCa                                             >  1:2023613/1‑93 (MQ=255)

                                            aTTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCa                                             >  2:4535615/1‑93 (MQ=255)

                                             ttCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCaa                                            <  1:3100909/93‑1 (MQ=255)

                                               cGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACg                                          <  1:3283058/93‑1 (MQ=255)

                                               cGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACg                                          <  2:1365069/93‑1 (MQ=255)

                                                 ccccAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGgc                                        <  2:61742/93‑1 (MQ=255)

                                                 ccccAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGgc                                        <  1:2419036/93‑1 (MQ=255)

                                                     aGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCAt                                    >  2:4634578/1‑93 (MQ=255)

                                                     aGAATTCCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCAt                                    <  1:543244/93‑1 (MQ=255)

                                                       aaTTTCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATaa                                  <  1:3991445/93‑1 (MQ=255)

                                                         tttCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATc                                >  1:841729/1‑93 (MQ=255)

                                                         tttCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATc                                >  1:3133784/1‑93 (MQ=255)

                                                         tttCTCCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATc                                <  2:1059865/93‑1 (MQ=255)

                                                         ttgctcCGCATCAAATACCGGCTGTGTTGCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATc                                <  1:2404994/90‑1 (MQ=255)

                                                             tcCGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGcc                            >  1:2176641/1‑93 (MQ=255)

                                                              ccGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCa                           <  2:1020807/93‑1 (MQ=255)

                                                              ccGCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCa                           <  2:3786660/93‑1 (MQ=255)

                                                                gCATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGc                         <  1:4283035/93‑1 (MQ=255)

                                                                 cATCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCa                        <  1:3404496/93‑1 (MQ=255)

                                                                  aTCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCTACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCAt                       >  2:4547968/1‑93 (MQ=255)

                                                                  aTCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCAt                       >  2:2792588/1‑93 (MQ=255)

                                                                  aTCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCAt                       >  1:1772001/1‑93 (MQ=255)

                                                                   tCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATc                      >  2:2264833/1‑93 (MQ=255)

                                                                   tCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATc                      >  1:1376916/1‑93 (MQ=255)

                                                                   tCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATc                      >  1:987126/1‑93 (MQ=255)

                                                                   tCAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATc                      <  1:4153494/93‑1 (MQ=255)

                                                                    cAAATACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCt                     >  2:1279845/1‑93 (MQ=255)

                                                                     aaaTACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTg                    >  1:1191721/1‑93 (MQ=255)

                                                                     aaaTACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTg                    >  2:741804/1‑93 (MQ=255)

                                                                      aaTACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGc                   <  2:272272/93‑1 (MQ=255)

                                                                      aaTACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGc                   <  1:3783792/93‑1 (MQ=255)

                                                                      aaTACCGGCTGTTTTTCCAGGCGACAAAAATAGAGGCCACAGCGACGAACGCTGTCGACTTTATTCAACGGCGAATAATCTGCCAGCATCTGc                   >  1:41170/1‑93 (MQ=255)

                                                                       aTACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCt                  >  2:4295481/1‑93 (MQ=255)

                                                                        ttctggCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTc                 <  2:1431446/89‑1 (MQ=255)

                                                                        tACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTc                 <  1:1104441/93‑1 (MQ=255)

                                                                        tACCGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTc                 >  1:4666073/1‑93 (MQ=255)

                                                                          ccGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCgg               <  2:4053336/93‑1 (MQ=255)

                                                                          ccGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCgg               >  1:71946/1‑93 (MQ=255)

                                                                          ccGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCgg               >  1:609301/1‑93 (MQ=255)

                                                                           cgtcTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGgc              <  2:3705438/90‑1 (MQ=255)

                                                                           cGGCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGgc              >  2:3842083/1‑93 (MQ=255)

                                                                            ggCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGgcg             >  1:1906526/1‑93 (MQ=255)

                                                                             gCTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGgcgc            >  2:819846/1‑93 (MQ=255)

                                                                              cTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCt           >  1:195686/1‑93 (MQ=255)

                                                                              cTGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCt           >  1:1148991/1‑93 (MQ=255)

                                                                               tGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTg          >  1:3336633/1‑93 (MQ=255)

                                                                               tGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTg          >  2:1887592/1‑93 (MQ=255)

                                                                               tGTTTTTCCAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTg          >  1:3872119/1‑93 (MQ=255)

                                                                                      ccAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCAcgc   >  2:2586968/1‑93 (MQ=255)

                                                                                       cAGGCGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCAcgcc  >  1:3474536/1‑93 (MQ=255)

AAACACGCCAGGCAGCGTTTTGACCGTCAGGCCATCGACGCTGTATTCGCCCCAGAATTTCTCCGCATCAAATACCGGCTGTTTTTCCAGACGACCAAAATAGAGGCCACAGCGACGAGCGCTGTCGACTTTATTCAACGGCGCATAATCTGCCAGCATCTGCTCGGCGCTGCGCACGCC  >  NC_000913/4604742‑4604921

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