BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,606,400	G→A	Q62Q (CAG→CAA)	yjjG →	putative phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,400	0	G	A	99.6%	986.4 / NA	269	Q62Q (CAG→CAA)	yjjG	putative phosphatase
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (148/120); minor base C (0/1); total (148/121)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 274 total aligned reads displayed.

CACCTTTACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGA  >  NC_000913/4606310‑4606489

cACCTTTACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACt                                                                                         <  1:4847123/93‑1 (MQ=255)

cACCTTTACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACt                                                                                         >  2:6482700/1‑93 (MQ=255)

 cccTTTCCCGCTGAAGATTTTAAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACtt                                                                                        <  1:1007100/92‑1 (MQ=255)

 aCCTTTACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACtt                                                                                        >  1:2303324/1‑93 (MQ=255)

  cttttACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTc                                                                                       <  1:1366160/91‑1 (MQ=255)

   cTTTACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCa                                                                                      >  2:2722784/1‑93 (MQ=255)

   cTTTACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCa                                                                                      <  2:1425853/93‑1 (MQ=255)

   cTTTACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCa                                                                                      >  1:4073648/1‑93 (MQ=255)

    tttACCGCTGAAGATTTTCAGGACTATCAGGACGTTAACAAGCCACTTTGGGTGGATTATCAAAACGGCGCGTTCACTTATTTACAACTTCac                                                                                     >  1:3531727/1‑92 (MQ=255)

     ttACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGc                                                                                    <  2:5682434/93‑1 (MQ=255)

        ccGCTGAAGTTTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACg                                                                                 <  2:4877338/93‑1 (MQ=255)

        ccGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACg                                                                                 >  1:5206074/1‑93 (MQ=255)

        ccGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACg                                                                                 >  1:4598705/1‑93 (MQ=255)

           cTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGc                                                                              >  2:4789475/1‑93 (MQ=255)

            tGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCg                                                                             >  1:3889862/1‑93 (MQ=255)

            tGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCg                                                                             >  2:2131374/1‑93 (MQ=255)

             gAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCgg                                                                            >  1:2041630/1‑93 (MQ=255)

               aGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGtt                                                                          >  1:3168280/1‑93 (MQ=255)

                gATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGttt                                                                         >  2:3035190/1‑93 (MQ=255)

                atttttGAGTATTATCAGGCCTTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGttt                                                                         <  1:3113208/91‑1 (MQ=255)

                 aTTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTg                                                                        >  1:6032646/1‑93 (MQ=255)

                  ttttCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTga                                                                       >  2:5415937/1‑93 (MQ=255)

                  ttttCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTga                                                                       <  2:903037/93‑1 (MQ=255)

                  ttttCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTga                                                                       >  2:334027/1‑93 (MQ=255)

                    ttCAGGATTACCAGTCCTTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTgaga                                                                     <  1:5421963/93‑1 (MQ=255)

                    ttCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTgaga                                                                     >  2:4005367/1‑93 (MQ=255)

                    ttCAGGACTATAAGCCGGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTgaga                                                                     <  2:1487817/93‑1 (MQ=255)

                     tCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGagag                                                                    >  1:816121/1‑93 (MQ=255)

                      cAGGACTATGAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGc                                                                   <  2:3505595/93‑1 (MQ=255)

                       aGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCt                                                                  >  2:1569311/1‑93 (MQ=255)

                        ggACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTg                                                                 >  2:3245703/1‑93 (MQ=255)

                        ggACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTg                                                                 >  2:119531/1‑93 (MQ=255)

                        ggACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTg                                                                 <  2:2859678/93‑1 (MQ=255)

                         gACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTgg                                                                <  2:6580479/93‑1 (MQ=255)

                           cTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAGACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGc                                                              <  2:5884916/93‑1 (MQ=255)

                           cTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGc                                                              <  1:6203295/93‑1 (MQ=255)

                           cTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGAGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGc                                                              >  2:4252582/1‑93 (MQ=255)

                            tATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTgggtc                                                             <  1:4200192/93‑3 (MQ=255)

                            tATCAGGACGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGcc                                                             >  2:1928341/1‑93 (MQ=255)

                                aGGCCTTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAAc                                                         <  2:5731376/93‑1 (MQ=255)

                                aGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAAc                                                         >  1:4141174/1‑93 (MQ=255)

                                 ggCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACg                                                        >  1:2437849/1‑93 (MQ=255)

                                 agCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACg                                                        >  1:3257670/2‑93 (MQ=255)

                                   ccGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGc                                                      <  1:6080612/93‑1 (MQ=255)

                                    cGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCt                                                     <  1:3041466/93‑1 (MQ=255)

                                      ttAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGa                                                   >  2:4525389/1‑93 (MQ=255)

                                       tAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGaa                                                  <  2:2326774/93‑1 (MQ=255)

                                        aaCAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAAc                                                 >  2:3389459/1‑93 (MQ=255)

                                         aCAAGCCGCTGTGGGTGGATTATCAAAACGGGGCGATCGCTTCATTACAACTTCACCACGGGCGGTTTGAGGGCTGGGCCGAACGGCTGAACg                                                >  1:4001464/1‑93 (MQ=255)

                                         aCAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACg                                                >  2:6712743/1‑93 (MQ=255)

                                           aaGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTc                                              <  2:3126937/93‑1 (MQ=255)

                                             gCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGa                                            <  1:6586844/93‑1 (MQ=255)

                                             gCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGa                                            <  1:645059/93‑1 (MQ=255)

                                                gcTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGcc                                         <  2:1787423/92‑1 (MQ=255)

                                                aCTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGcc                                         <  2:1301412/93‑1 (MQ=255)

                                                aCTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGcc                                         >  2:3301591/1‑93 (MQ=255)

                                                aCTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGcc                                         >  1:2776289/1‑93 (MQ=255)

                                                 cTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGGCGAGccg                                        >  1:2751314/1‑92 (MQ=255)

                                                  tgtgGGTGGATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAg                                       <  1:3833336/93‑1 (MQ=255)

                                                   gtgggtggATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAgg                                      <  2:3921287/93‑1 (MQ=255)

                                                   gtgggtggATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAgg                                      <  1:4376409/93‑1 (MQ=255)

                                                     gggtggATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTa                                    >  2:4627493/1‑93 (MQ=255)

                                                      ggtggATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTaa                                   >  1:1184577/1‑93 (MQ=255)

                                                      ggtggATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTaa                                   >  1:5767565/1‑93 (MQ=255)

                                                        tggATTATCAAAACGGCGCGATCGCTTCATTACACCTGCAGCACGGGCGGTTGGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAAc                                 <  1:464999/93‑1 (MQ=255)

                                                        tggATTATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAAc                                 >  1:2900937/1‑93 (MQ=255)

                                                            ttATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCaa                             >  2:3719921/1‑93 (MQ=255)

                                                             tATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGGGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAAc                            >  2:3277862/1‑93 (MQ=255)

                                                             tATCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAAc                            <  1:1730663/93‑1 (MQ=255)

                                                              aTCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACg                           >  2:6236232/1‑93 (MQ=255)

                                                               tCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGa                          <  2:2791009/93‑1 (MQ=255)

                                                               tCAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGa                          >  1:5710574/1‑93 (MQ=255)

                                                                cAAAACGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGaa                         <  1:1250514/93‑1 (MQ=255)

                                                                 aaaaCGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAg                        >  1:5989803/1‑93 (MQ=255)

                                                                 aaaaCGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAg                        <  2:3884150/93‑1 (MQ=255)

                                                                 aaaaCGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAg                        >  2:1694091/1‑93 (MQ=255)

                                                                 aaaaCGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAg                        >  1:2698481/1‑93 (MQ=255)

                                                                  aaaCGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGc                       <  1:2143585/93‑1 (MQ=255)

                                                                  aaaCGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGc                       >  2:4757956/1‑93 (MQ=255)

                                                                   aaCGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGcc                      >  1:1892277/1‑93 (MQ=255)

                                                                    aCGGCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCt                     <  2:4090020/93‑1 (MQ=255)

                                                                      ggCGCGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCttt                   <  1:700788/93‑1 (MQ=255)

                                                                       gcgcGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTtta                  <  1:1623565/93‑1 (MQ=255)

                                                                       gcgcGATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTtta                  >  2:5071450/1‑93 (MQ=255)

                                                                         gcgATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTttatt                >  1:6112521/1‑93 (MQ=255)

                                                                           gATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTaa              <  1:6523924/93‑1 (MQ=255)

                                                                           gATCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGACGAGCCAGGTACACTCAACGAAGCCTTTATTaa              >  2:267094/1‑93 (MQ=255)

                                                                            aTCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAAt             <  1:2809269/93‑1 (MQ=255)

                                                                            aTCACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAAt             >  1:1313338/1‑93 (MQ=255)

                                                                              cACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGc           <  1:5089034/93‑1 (MQ=255)

                                                                              cACTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGc           <  2:4181514/93‑1 (MQ=255)

                                                                                cTTCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGa         <  2:3004627/93‑1 (MQ=255)

                                                                                 ttCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGAt        >  1:5617899/1‑93 (MQ=255)

                                                                                 ttCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGAt        <  2:356849/93‑1 (MQ=255)

                                                                                 ttCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGAt        >  2:106144/1‑93 (MQ=255)

                                                                                  tCATTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATg       >  1:2179829/1‑93 (MQ=255)

                                                                                    aTTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGc     >  2:5225212/1‑93 (MQ=255)

                                                                                    aTTACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGc     <  1:5027867/93‑1 (MQ=255)

                                                                                      tACAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCgg   >  2:5787607/1‑93 (MQ=255)

                                                                                       aCAACTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGa  <  2:6006968/93‑1 (MQ=255)

CACCTTTACCGCTGAAGATTTTCAGGACTATCAGGCCGTTAACAAGCCACTGTGGGTGGATTATCAAAACGGCGCGATCACTTCATTACAGCTTCAGCACGGGCGGTTTGAGAGCTGGGCCGAACGGCTGAACGTCGAGCCAGGTAAACTCAACGAAGCCTTTATTAATGCGATGGCGGA  >  NC_000913/4606310‑4606489

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