BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,873,733	G→A	R204R (CGG→CGA)	yidX →	putative replicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,873,733	0	G	A	99.8%	1651.6 / NA	450	R204R (CGG→CGA)	yidX	putative replicase
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (223/225); minor base C (1/0); total (225/225)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 460 total aligned reads displayed.

CATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGGGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTCAAGGTTAGATTTG  >  NC_000913/3873643‑3873825

  ttACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCaa                                                                                          >  1:3880553/1‑93 (MQ=255)

  ttACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCaa                                                                                          <  1:7294818/93‑1 (MQ=255)

   tccTGCCGGTGACTGTTCCGCTGGCCATCATTCCTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCaaa                                                                                         <  2:7000861/91‑1 (MQ=255)

    aCTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAAc                                                                                        <  2:7869356/93‑1 (MQ=255)

    aCTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAAc                                                                                        <  1:4483351/93‑1 (MQ=255)

    aCTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAAc                                                                                        <  1:7515858/93‑1 (MQ=255)

      tGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAAcgc                                                                                      <  2:7199029/93‑1 (MQ=255)

       gCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAAcgcg                                                                                     <  2:3306469/93‑1 (MQ=255)

        ccAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGt                                                                                    <  2:4550148/93‑1 (MQ=255)

          aGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGttt                                                                                  >  2:1850319/1‑93 (MQ=255)

           gTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTa                                                                                 >  1:681373/1‑93 (MQ=255)

           gTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTa                                                                                 <  2:4893607/93‑1 (MQ=255)

           gTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTa                                                                                 <  1:4582625/93‑1 (MQ=255)

            tGACTGTTACGCTGGACATCA‑TACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGGTTATg                                                                               >  1:1875769/1‑93 (MQ=255)

              aCTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGc                                                                              >  1:5319309/1‑93 (MQ=255)

              aCTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGc                                                                              >  2:5926615/1‑93 (MQ=255)

              aCTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGc                                                                              <  1:1588456/93‑1 (MQ=255)

                tGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTg                                                                            >  1:2997923/1‑93 (MQ=255)

                tGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTg                                                                            >  1:2626290/1‑93 (MQ=255)

                 gTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGt                                                                           <  2:6149109/93‑1 (MQ=255)

                  ttACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGtt                                                                          <  1:6962987/93‑1 (MQ=255)

                  ttACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGtt                                                                          <  2:5231409/93‑1 (MQ=255)

                  ttACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGtt                                                                          >  1:5620920/1‑93 (MQ=255)

                  ttACGCTGGACATCATTACTGCACCGCTGCAATGTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGtt                                                                          <  1:7607668/93‑1 (MQ=255)

                        tGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTgcccg                                                                    >  1:3188772/1‑93 (MQ=255)

                          gACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTgcccgcc                                                                  <  2:271355/93‑1 (MQ=255)

                          gACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTgcccgcc                                                                  <  1:5897300/93‑1 (MQ=255)

                             atcatTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCcac                                                               >  2:422102/1‑93 (MQ=255)

                              tcatTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCcaca                                                              >  1:3360760/1‑93 (MQ=255)

                              tcatTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGGAAGCGAGGAAACGCGTTTATGCTGGTTGGCCGCCcaca                                                              >  2:3718407/1‑93 (MQ=255)

                                atTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGa                                                            <  2:4463535/93‑1 (MQ=255)

                                atTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGa                                                            >  1:7035187/1‑93 (MQ=255)

                                atTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGa                                                            >  2:7336529/1‑93 (MQ=255)

                                 gtACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCGGTTTGCCCGCCCACAGAt                                                           <  2:7078495/92‑1 (MQ=255)

                                    cTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTaa                                                        >  1:6123346/1‑93 (MQ=255)

                                      gCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAAtt                                                      >  2:1421561/1‑93 (MQ=255)

                                        aCCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCa                                                    <  1:7800319/93‑1 (MQ=255)

                                          cGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGc                                                  <  1:341414/93‑1 (MQ=255)

                                           gCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGca                                                 >  1:6545172/1‑93 (MQ=255)

                                           gCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGca                                                 >  1:7431757/1‑93 (MQ=255)

                                            cTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCTTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGcac                                                <  2:1606627/93‑1 (MQ=255)

                                             tGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGcaca                                               <  1:5179664/93‑1 (MQ=255)

                                             tGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGcaca                                               <  1:4659921/93‑1 (MQ=255)

                                             tGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGcaca                                               >  2:2918242/1‑93 (MQ=255)

                                             tGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGcaca                                               <  1:197164/93‑1 (MQ=255)

                                               cAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATa                                             <  2:3831774/93‑1 (MQ=255)

                                                aaTTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATAc                                            <  2:1211518/93‑1 (MQ=255)

                                                  tttCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACtt                                          <  2:2271604/93‑1 (MQ=255)

                                                  tttCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACtt                                          >  1:1065559/1‑93 (MQ=255)

                                                  tttCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACtt                                          >  2:5314868/1‑93 (MQ=255)

                                                   ttCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTtc                                         >  1:6007484/1‑93 (MQ=255)

                                                     cTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTtctc                                       <  2:5690679/93‑1 (MQ=255)

                                                      tGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCa                                      >  2:6274687/1‑93 (MQ=255)

                                                      tGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCa                                      >  1:3550027/1‑93 (MQ=255)

                                                       ggTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCaa                                     >  2:5028013/1‑93 (MQ=255)

                                                       ggTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCaa                                     >  1:3961760/1‑93 (MQ=255)

                                                         ttGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATa                                   >  1:1038810/1‑93 (MQ=255)

                                                         ttGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATa                                   <  2:3534627/93‑1 (MQ=255)

                                                         ttGTATATGCAGTAAACCAATAATCAGTAAGCAAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATa                                   <  1:5847122/93‑1 (MQ=255)

                                                          tGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAg                                  >  1:721624/1‑93 (MQ=255)

                                                          tGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAg                                  >  2:80242/1‑93 (MQ=255)

                                                          tGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAg                                  >  1:497346/1‑93 (MQ=255)

                                                          tGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAg                                  <  1:3513917/93‑1 (MQ=255)

                                                           gTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGc                                 <  1:4156092/93‑1 (MQ=255)

                                                            tataTGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCa                                >  2:377174/1‑93 (MQ=255)

                                                             atatGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCaa                               >  2:4805771/1‑93 (MQ=255)

                                                              tatGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCaaa                              >  2:6875540/1‑93 (MQ=255)

                                                              tatGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCaaa                              <  1:954074/93‑1 (MQ=255)

                                                              tatGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCaaa                              >  1:2222922/1‑93 (MQ=255)

                                                               atGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAac                             <  2:4152498/93‑1 (MQ=255)

                                                               atGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAac                             <  1:7226602/93‑1 (MQ=255)

                                                                tGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAacg                            >  1:1409243/1‑93 (MQ=255)

                                                                   aGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAacgcca                         >  1:2778420/1‑93 (MQ=255)

                                                                    gTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAacgccac                        >  2:7630877/1‑93 (MQ=255)

                                                                     tAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAacgccacg                       >  1:205279/1‑93 (MQ=255)

                                                                        aCCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAAcgccacgcca                    <  2:6050428/93‑1 (MQ=255)

                                                                        aCCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAAcgccacgcca                    <  2:1066928/93‑1 (MQ=255)

                                                                        aCCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAAcgccacgcca                    <  1:4803891/93‑1 (MQ=255)

                                                                         ccAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCAt                   >  2:407464/1‑93 (MQ=255)

                                                                          cAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATc                  >  1:841933/1‑93 (MQ=255)

                                                                           aataatCAGTAAGCGCGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACAAACTTCTCAATAGCAAACGCCCCGCCATCt                 >  2:2763028/1‑93 (MQ=255)

                                                                            ataatCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCtt                >  1:1804256/1‑93 (MQ=255)

                                                                            ataatCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCtt                >  2:6709868/1‑93 (MQ=255)

                                                                            ataatCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCtt                <  2:115884/93‑1 (MQ=255)

                                                                             taatCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTc               >  2:3126029/1‑93 (MQ=255)

                                                                              aatCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTCa              >  1:6403641/1‑93 (MQ=255)

                                                                               atCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTCaa             <  2:5493605/93‑1 (MQ=255)

                                                                               atCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTCaa             <  1:1247289/93‑1 (MQ=255)

                                                                               atCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAACAGCAAACGCCACGCCATCTTCaa             <  2:4677728/93‑1 (MQ=255)

                                                                                tCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTCAAg            <  1:2435336/93‑1 (MQ=255)

                                                                                 cAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTCAAgg           >  1:4262128/1‑93 (MQ=255)

                                                                                  aGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTCAAGGt          >  2:7510682/1‑93 (MQ=255)

                                                                                  aGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTCAAGGt          >  2:39536/1‑93 (MQ=255)

                                                                                    tAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTCAAGGTTa        >  1:6302806/1‑93 (MQ=255)

                                                                                     aaGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTCAAGGTTAg       <  2:217292/93‑1 (MQ=255)

                                                                                      aGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTCAAGGTTAGa      >  2:580638/1‑93 (MQ=255)

                                                                                      aGCAGGCAAAGGCGTTTTGCCTGTTTGCCCGCCCCCAGATTAATTAAGACCTTACTTTTCATAACCAACCGCCCCGCCATCTTCAAGGttggt      >  2:511529/1‑90 (MQ=255)

                                                                                       gCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTCAAGGTTAGAt     >  2:896326/1‑93 (MQ=255)

                                                                                       gCGAGCAAACGCGTTTATGCTGTGTGCCCGCCCACAGGTTAATTCAGCACATTCTTCTCAATAGCAAACGCCACGCCATCTTCAAGGTTAGAt     <  1:382826/93‑1 (MQ=255)

                                                                                          ggCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTCAAGGTTAGATTTg  <  2:6534686/93‑1 (MQ=255)

CATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGGGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTCAAGGTTAGATTTG  >  NC_000913/3873643‑3873825

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