BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,873,493	T→C	C124C (TGT→TGC)	yidX →	putative replicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,873,493	0	T	C	99.6%	918.0 / ‑5.1	281	C124C (TGT→TGC)	yidX	putative replicase
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (132/147); total (134/147)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 285 total aligned reads displayed.

CGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGTACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCAT  >  NC_000913/3873403‑3873581

 ggCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCAcc                                                                                       <  1:1216672/93‑1 (MQ=255)

  gCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCa                                                                                      >  2:4795915/1‑93 (MQ=255)

   cacaATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAg                                                                                     <  2:6245476/93‑1 (MQ=255)

     caATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGt                                                                                   >  1:5712947/1‑93 (MQ=255)

       atatCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCa                                                                                 >  2:3902194/1‑93 (MQ=255)

        tatCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCAc                                                                                >  2:3654354/1‑93 (MQ=255)

           cGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTgg                                                                             <  2:277557/93‑1 (MQ=255)

            gCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGa                                                                            >  2:2908243/1‑93 (MQ=255)

             cTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGaa                                                                           >  1:3575674/1‑93 (MQ=255)

              tttACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGTACCAGGTCACTGGaaa                                                                          >  2:5064340/1‑93 (MQ=255)

              tttACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGaaa                                                                          >  1:2460951/1‑93 (MQ=255)

               ttACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGaaaa                                                                         <  2:4280749/93‑1 (MQ=255)

               ttACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGaaaa                                                                         >  1:4506536/1‑93 (MQ=255)

                tactactacTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAAc                                                                        <  2:2297702/93‑1 (MQ=255)

                 actactactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAAcc                                                                       >  1:535534/1‑93 (MQ=255)

                  ctactactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCt                                                                      >  2:151960/1‑93 (MQ=255)

                   tactactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTa                                                                     >  2:3310516/1‑93 (MQ=255)

                   tactactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTa                                                                     <  2:4372987/93‑1 (MQ=255)

                   tactactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTa                                                                     >  2:5821533/1‑93 (MQ=255)

                   tactactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTa                                                                     <  2:909541/93‑1 (MQ=255)

                    actactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaa                                                                    <  2:3137535/93‑1 (MQ=255)

                    actactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaa                                                                    >  1:4841985/1‑93 (MQ=255)

                    actactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaa                                                                    <  1:3015263/93‑1 (MQ=255)

                     ctactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaaa                                                                   <  1:4151909/93‑1 (MQ=255)

                      tactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaaaa                                                                  <  2:3112105/93‑1 (MQ=255)

                      tactGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTaaaa                                                                  <  1:3790880/93‑1 (MQ=255)

                       actGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAg                                                                 <  2:31938/93‑1 (MQ=255)

                       actGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAg                                                                 <  1:4972388/93‑1 (MQ=255)

                         tGGAATAGCGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCGTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGc                                                               <  2:2967860/93‑1 (MQ=255)

                         tGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGc                                                               <  1:5964561/93‑1 (MQ=255)

                         tGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGc                                                               >  2:5453711/1‑93 (MQ=255)

                         tGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGc                                                               >  2:4445564/1‑93 (MQ=255)

                          ggAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGca                                                              >  2:1462753/1‑93 (MQ=255)

                            aaTAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGcaca                                                            <  2:5657541/93‑1 (MQ=255)

                             aTAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACaa                                                           >  2:5386376/1‑93 (MQ=255)

                             aTAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACaa                                                           <  1:2573526/93‑1 (MQ=255)

                             aTAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACaa                                                           >  2:2320103/1‑93 (MQ=255)

                             aTAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACaa                                                           >  2:1987142/1‑93 (MQ=255)

                               aaCGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATc                                                         >  2:3557446/1‑93 (MQ=255)

                                 cGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCa                                                       <  1:1367095/93‑1 (MQ=255)

                                 cGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCa                                                       <  2:553554/93‑1 (MQ=255)

                                 cGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCa                                                       >  1:3250798/1‑93 (MQ=255)

                                  gaagaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAACCCTAAAAGGCACAATCCAt                                                      <  2:6081728/93‑1 (MQ=255)

                                  gaagaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCAt                                                      <  2:3409588/93‑1 (MQ=255)

                                  gaagaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCAt                                                      <  1:1006908/93‑1 (MQ=255)

                                   aagaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATg                                                     <  2:4491087/93‑1 (MQ=255)

                                   aagaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATg                                                     <  2:4331098/93‑1 (MQ=255)

                                    agaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCGCCAGGTCACTGGAAAACCTAAAAGGCACAATCCATga                                                    <  1:1911739/93‑1 (MQ=255)

                                    agaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATga                                                    <  1:3333472/93‑1 (MQ=255)

                                     gaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATgag                                                   <  2:1319946/93‑1 (MQ=255)

                                     gaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATgag                                                   <  1:1173589/93‑1 (MQ=255)

                                     gaagaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATgag                                                   >  2:523107/1‑93 (MQ=255)

                                        gaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGCAAACCTAAAAGGCACAATCCATGAGaaa                                                >  2:2637391/1‑93 (MQ=255)

                                        gaaAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGaaa                                                <  1:1289599/93‑1 (MQ=255)

                                         aaaaGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGaaaa                                               <  1:3663079/93‑1 (MQ=255)

                                             ggCACTGGCAACTCATTATGGTTTTGCCTGTGGTGTTCAACACGGCATCAGGTCACTGGAGAACCTAAAAGGCACAATCCATGAGAAAAATaa                                           <  2:1006385/93‑1 (MQ=255)

                                             ggCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATaa                                           >  1:2915618/1‑93 (MQ=255)

                                              gCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATaaa                                          <  1:2712341/93‑1 (MQ=255)

                                              gCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATaaa                                          >  2:778204/1‑93 (MQ=255)

                                               cACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATaaaa                                         <  1:39940/93‑1 (MQ=255)

                                                aCTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAACAGGCAGAATCCATGAGAAAAATaaaaa                                        >  1:5012861/1‑93 (MQ=255)

                                                 cTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAAc                                       >  1:4054795/1‑93 (MQ=255)

                                                 cTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAAc                                       >  2:4201517/1‑93 (MQ=255)

                                                 cTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAAc                                       >  1:5223576/1‑93 (MQ=255)

                                                 cTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAAc                                       <  2:5091624/93‑1 (MQ=255)

                                                     cAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATgg                                   <  1:3117835/93‑1 (MQ=255)

                                                     cAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATgg                                   <  1:5845844/93‑1 (MQ=255)

                                                     cAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATgg                                   <  1:5861164/93‑1 (MQ=255)

                                                      aaCGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGa                                  <  2:1749272/93‑1 (MQ=255)

                                                      aaCGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGa                                  <  2:5016626/93‑1 (MQ=255)

                                                       aCGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGAt                                 >  2:4903249/1‑93 (MQ=255)

                                                        cGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGAtt                                <  1:5677531/93‑1 (MQ=255)

                                                        cGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGAtt                                <  1:2631129/93‑1 (MQ=255)

                                                        cGCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGAtt                                >  2:2565612/1‑93 (MQ=255)

                                                         gCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTa                               >  1:3760179/1‑93 (MQ=255)

                                                         gCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTa                               >  2:3491438/1‑93 (MQ=255)

                                                         gCATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTa                               >  1:1613218/1‑93 (MQ=255)

                                                          cATTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTAc                              <  2:5712947/93‑1 (MQ=255)

                                                           aTTATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACt                             <  2:4668126/93‑1 (MQ=255)

                                                            ttATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTc                            <  2:5786152/93‑1 (MQ=255)

                                                             tATGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCa                           >  2:2403604/1‑93 (MQ=255)

                                                              aTGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCaa                          >  1:6246164/1‑93 (MQ=255)

                                                              aTGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCaa                          >  2:3972118/1‑93 (MQ=255)

                                                               tGGTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCaaa                         <  2:4088332/93‑1 (MQ=255)

                                                                ggTTTTGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAg                        <  1:2752796/93‑1 (MQ=255)

                                                                  ttttGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGt                      >  1:2264055/1‑93 (MQ=255)

                                                                   tttGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTg                     <  1:4458029/93‑1 (MQ=255)

                                                                    ttGCCTGTGGTGTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGa                    <  1:738372/93‑1 (MQ=255)

                                                                          gtggtgTTCAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGt              >  1:155145/1‑93 (MQ=255)

                                                                                 tCAACACTGCACCAGGTCATTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAcca       >  1:1694684/1‑93 (MQ=255)

                                                                                  cAACACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccat      >  1:6165059/1‑93 (MQ=255)

                                                                                   agcacTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATCACTCAAAGGTGATGGCGTTCTAccatc     <  2:2079115/91‑1 (MQ=255)

                                                                                   aaCACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatc     >  2:2246300/1‑93 (MQ=255)

                                                                                   aaCACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatc     <  1:1504627/93‑1 (MQ=255)

                                                                                   aaCACTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatc     <  2:2669444/93‑1 (MQ=255)

                                                                                     cacTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatcca   <  1:3966741/93‑1 (MQ=255)

                                                                                     cacTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatcca   >  1:4409372/1‑93 (MQ=255)

                                                                                     cacTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatcca   >  2:1117732/1‑93 (MQ=255)

                                                                                     cacTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatcca   <  2:1868444/93‑1 (MQ=255)

                                                                                      acTGCACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTAccatccat  >  1:2116870/1‑93 (MQ=255)

CGGCACAATATCGCTTTACTACTACTGGAATAACGAAGAAGAAAAGGCACTGGCAACGCATTATGGTTTTGCCTGTGGTGTTCAACACTGTACCAGGTCACTGGAAAACCTAAAAGGCACAATCCATGAGAAAAATAAAAACATGGATTACTCAAAGGTGATGGCGTTCTACCATCCAT  >  NC_000913/3873403‑3873581

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