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.6%	1007.4 / ‑4.2	275	R204R (CGG→CGA)	yidX	putative replicase
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (139/135); minor base C (1/0); total (140/135)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.73e-01

Only 100 of 283 total aligned reads displayed.

CATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGGGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTCAAGGTTAGAT  >  NC_000913/3873643‑3873822

      tGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAAcgc                                                                                   >  2:168700/1‑93 (MQ=255)

      tGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGAAGTAAACCAATAATCAGTAAGCGATCAAAcgc                                                                                   >  2:3031774/1‑93 (MQ=255)

       gCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAAcgcg                                                                                  >  2:1160783/1‑93 (MQ=255)

        ccAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGt                                                                                 >  1:1472727/1‑93 (MQ=255)

         cAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGtt                                                                                <  1:4006190/93‑1 (MQ=255)

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

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

            tGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTAt                                                                             >  1:2342777/1‑93 (MQ=255)

            tGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTAt                                                                             >  2:3228177/1‑93 (MQ=255)

             gACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATg                                                                            >  2:1103843/1‑93 (MQ=255)

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

               cTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCt                                                                          <  1:3120552/93‑1 (MQ=255)

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

                 gTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGt                                                                        >  2:2741063/1‑93 (MQ=255)

                 gTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGt                                                                        >  1:963560/1‑93 (MQ=255)

                   tACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGttt                                                                      <  2:1087100/93‑1 (MQ=255)

                   tACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGttt                                                                      <  2:1931010/93‑1 (MQ=255)

                    aCGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTg                                                                     >  1:3816475/1‑93 (MQ=255)

                      gCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTgcc                                                                   <  2:2493007/93‑1 (MQ=255)

                      gCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTCTATGCTGTTTgcc                                                                   <  2:4603939/93‑1 (MQ=255)

                       cTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTgccc                                                                  >  1:859975/1‑93 (MQ=255)

                       cTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTgccc                                                                  >  1:4066490/1‑93 (MQ=255)

                       cTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTgccc                                                                  <  2:3100680/93‑1 (MQ=255)

                         ggACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTgcccgc                                                                <  1:2616210/93‑1 (MQ=255)

                           aCATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTgcccgccc                                                              >  2:3288463/1‑93 (MQ=255)

                            catcatTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCca                                                             <  1:1245631/93‑1 (MQ=255)

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

                              tcatTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCcaca                                                           <  2:3885994/93‑1 (MQ=255)

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

                               catTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAg                                                          >  2:642731/1‑93 (MQ=255)

                                 ttACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGAt                                                        >  1:3320544/1‑93 (MQ=255)

                                 ttACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGAt                                                        <  1:1861434/93‑1 (MQ=255)

                                 ttACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGAt                                                        >  1:2715300/1‑93 (MQ=255)

                                 ttACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGAt                                                        >  1:3283230/1‑93 (MQ=255)

                                 ttACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGAt                                                        <  1:3057750/93‑1 (MQ=255)

                                    cTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTaa                                                     <  2:4180356/93‑1 (MQ=255)

                                       cACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTc                                                  <  1:3848965/93‑1 (MQ=255)

                                       cACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTc                                                  <  1:2993787/93‑1 (MQ=255)

                                       cACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTc                                                  <  2:2882734/93‑1 (MQ=255)

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

                                         ccGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAg                                                >  2:4413361/1‑93 (MQ=255)

                                         ccGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAg                                                <  2:437645/93‑1 (MQ=255)

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

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

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

                                              gCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACAt                                           >  1:32022/1‑93 (MQ=255)

                                                aaTTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCACACGCGTTTATGCTGTTTGCCCGACCACAGATTAATTCAGCACATAc                                         >  1:1581384/1‑93 (MQ=255)

                                                aaTTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATAc                                         >  2:4322773/1‑93 (MQ=255)

                                                 aTTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACt                                        >  1:834454/1‑93 (MQ=255)

                                                 aTTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACt                                        <  1:1142907/93‑1 (MQ=255)

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

                                                   ttCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTtc                                      <  2:1270684/93‑1 (MQ=255)

                                                    tCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTtct                                     >  1:4239863/1‑93 (MQ=255)

                                                    tCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTtct                                     >  2:4675396/1‑93 (MQ=255)

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

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

                                                        gTTGTATATGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAAt                                 >  2:185471/1‑93 (MQ=255)

                                                         ttGTATATGCAGTAAACCAATAATCAGTAAGCGATCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATa                                <  2:4331749/93‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                tGCAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAacg                         <  2:1415927/93‑1 (MQ=255)

                                                                  cAGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAacgcc                       <  1:1959835/93‑1 (MQ=255)

                                                                   aGTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAacgcca                      <  1:1414898/93‑1 (MQ=255)

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

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

                                                                    gTAAACCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAacgccac                     <  2:386396/93‑1 (MQ=255)

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

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

                                                                     tAAACCAATAATCAGTAAGAGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAacgccacg                    >  2:4569095/1‑93 (MQ=255)

                                                                       aaCCAATAATCAGTAAGCGAGCAAAGGCGTTTATGCTGTTTGCCCGCCCACAGCTTAATTCAGCACATACTTCTCAATAGCAAacgccacgcc                  >  2:1431904/1‑93 (MQ=255)

                                                                       aaCCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAacgccacgcc                  >  2:300161/1‑93 (MQ=255)

                                                                       aaCCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAacgccacgcc                  <  1:2220368/93‑1 (MQ=255)

                                                                        aCCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAAcgccacgcca                 >  1:1553947/1‑93 (MQ=255)

                                                                        aCCAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAAcgccacgcca                 >  1:4461669/1‑93 (MQ=255)

                                                                         ccAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCACACAGATTAATTCAGGACATACTTCTCAATAGCAAACGCAACGACAt                >  2:796133/1‑93 (MQ=255)

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

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

                                                                          cAATAATCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATc               <  2:3380969/93‑1 (MQ=255)

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

                                                                           aataatCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCt              >  1:1102702/1‑93 (MQ=255)

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

                                                                             taatCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTc            <  1:3671860/93‑1 (MQ=255)

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

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

                                                                             taatCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTc            <  1:4514046/93‑1 (MQ=255)

                                                                               atCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGACCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTCaa          >  1:2829445/1‑93 (MQ=255)

                                                                                tCAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTCAAg         >  1:2781452/1‑93 (MQ=255)

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

                                                                                 cAGTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTCAAgg        <  1:1761674/93‑1 (MQ=255)

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

                                                                                   gTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTCAAGGtt      >  1:1822265/1‑93 (MQ=255)

                                                                                   gTAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTCAAGGtt      >  1:3777135/1‑93 (MQ=255)

                                                                                    tAAGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTCAAGGTTa     <  2:138465/93‑1 (MQ=255)

                                                                                      aGCGAGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTCAAGGTTAGa   <  1:1535546/93‑1 (MQ=255)

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

CATTACTGCCAGTGACTGTTACGCTGGACATCATTACTGCACCGCTGCAATTTCTGGTTGTATATGCAGTAAACCAATAATCAGTAAGCGGGCAAACGCGTTTATGCTGTTTGCCCGCCCACAGATTAATTCAGCACATACTTCTCAATAGCAAACGCCACGCCATCTTCAAGGTTAGAT  >  NC_000913/3873643‑3873822

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